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Vorwort des Herausgebers 


Der schnelle technische Fortschritt im Turbomaschinenbau, der durch extreme technische Forde- 
rungen und starken internationalen Wettbewerb geprägt ist, verlangt einen effizienten Austausch 
und die Diskussion von Fachwissen und Erfahrung zwischen Universitäten und industriellen 
Partnern. Mit der vorliegenden Reihe haben wir versucht, ein Forum zu schaffen, das neben 
unseren Publikationen in Fachzeitschriften die aktuellen Forschungsergebnisse des Instituts für 
Thermische Strömungsmaschinen am Karlsruher Institut für Technologie (KIT) einem möglichst 
großen Kreis von Fachkollegen aus der Wissenschaft und vor allem auch der Praxis zugänglich 
macht und den Wissenstransfer intensiviert und beschleunigt. 


Flugtriebwerke, stationäre Gasturbinen, Turbolader und Verdichter sind im Verbund mit den zuge- 
hörigen Anlagen faszinierende Anwendungsbereiche. Es ist nur natürlich, dass die methodischen 
Lösungsansätze, die neuen Messtechniken, die Laboranlagen auch zur Lösung von Problemstel- 
lungen in anderen Gebieten - hier denke ich an Otto- und Dieselmotoren, elektrische Antriebe 
und zahlreiche weitere Anwendungen - genutzt werden. Die effiziente, umweltfreundliche und 
zuverlässige Umsetzung von Energie führt zu Fragen der ein- und mehrphasigen Strömung, der 
Verbrennung und der Schadstoffbildung, des Wärmeübergangs sowie des Verhaltens metallischer 
und keramischer Materialien und Verbundwerkstoffe. Sie stehen im Mittelpunkt ausgedehn- 
ter theoretischer und experimenteller Arbeiten, die im Rahmen nationaler und internationaler 
Forschungsprogramme in Kooperation mit Partnern aus Industrie, Universitäten und anderen 
Forschungseinrichtungen durchgeführt werden. 


Es sollte nicht unerwähnt bleiben, dass alle Arbeiten durch enge Kooperation innerhalb des Insti- 
tuts geprägt sind. Nicht ohne Grund ist der Beitrag der Werkstätten, der Technik-, der Rechner- 
und Verwaltungsabteilungen besonders hervorzuheben. Diplomanden und Hilfsassistenten tra- 
gen mit ihren Ideen Wesentliches bei, und natürlich ist es der stets freundschaftlich fordernde 
wissenschaftliche Austausch zwischen den Forschergruppen des Instituts, der zur gleichbleibend 
hohen Qualität der Arbeiten entscheidend beiträgt. Dabei sind wir für die Unterstützung unserer 
Förderer außerordentlich dankbar. 


Die in jüngster Zeit in ersten Turbofantriebwerken vollzogene Änderung der Triebwerksarchi- 
tektur, die Integration eines Planetengetriebes zwischen dem Gebläse (Fan) des Triebwerks und 
der Niederdruckwelle, führt zu einer drastischen Reduktion des Brennstoffverbrauchs in der 
Größenordnung von 10 - 15% im Vergleich zur vorigen Triebwerksgeneration. Eine wesentliche 
Herausforderung bei der Auslegung des Untersetzungsgetriebes von Triebwerken für Großraum- 
Langstreckenflugzeuge besteht in der Beherrschung seines Thermalhaushalts. Aufgrund der 
immensen übertragenen Leistung führen selbst geringe relative mechanische Verluste des Ge- 
triebes zur Dissipation einer Leistung von mehreren hundert Kilowatt. Die dissipierte Leistung 
ist durch eine effektive und effiziente Kühlung der Zahnräder aus dem Getriebe abzuleiten. 
Dabei steht die Kühlung, die durch Besprühen der Zahnflanken mit Ölstrahlen erfolgt, vor der 
Herausforderung, eine ausreichende Temperaturabsenkung der Zahnräder sicherzustellen und 
gleichzeitig aber die eingesetzte Öl€menge zu begrenzen, um keine zusätzlichen Planschverluste 
zu generieren. Um dieses Ziel erreichen zu können, ist die genaue Kenntnis des Wärmeüber- 


gangs zwischen Olstrahl und Zahnflanke unbedingte Voraussetzung. Im vorliegenden Band 
der Schriftenreihe befasst sich der Autor folgerichtig mit der experimentellen Untersuchung 
des Wärmeübergangs zwischen Ölstrahlen und den Flanken eines repräsentativen Zahnrads. 
Die Ergebnisse der ortsaufgelösten Messungen des Wärmeübergangs ermöglichen, die Wirk- 
zusammenhänge zwischen den Randbedingungen Olvolumenstrom, Drehzahl, Diisenwinkel, 
Strahleindringtiefe und Aufprallwinkel und dem mittleren Wärmeübergang zu identifizieren und 
physikalisch verständlich zu erklären. Im Gegensatz zu bisher verwendeten Korrelationen zur 
Abschätzung des Wärmeübergangs wird erstmals der dominierende Effekt der Prallkühlung des 
auf die Zahnflanke auftreffenden Ölstrahls erfasst. Die hohe Qualität der aus den Messungen 
abgeleiteten Korrelation für den mittleren Wärmeübergang ermöglicht nun eine erheblich siche- 
rere wärmetechnische Auslegung von Hochleistungsgetrieben für zukünftige Flugtriebwerke mit 
gesteigerten Nebenstromverhältnissen und höchsten Vortriebswirkungsgraden. 


Karlsruhe, im September 2021 Hans-Jörg Bauer 


Vorwort des Autors 


Die vorliegende Arbeit ist während meiner Mitarbeit am Institut für Thermische Strömungsma- 
schinen (ITS) am Karlsruher Institut für Technologie (KIT) entstanden. Das zugrunde liegende 
Forschungsprojekt zur Untersuchung der Einspritzkühlung von Leistungsgetrieben in Flug- 
triebwerken wurde in enger Zusammenarbeit mit Rolls-Royce Deutschland im Rahmen des 
Luftfahrtforschungsprogramms des Bundesministeriums für Wirtschaft und Energie (BMWi) 
durchgeführt. 


Im Gegensatz zu der linearen, deterministischen und kühlen Form, die das Endprodukt Disser- 
tation naturgemäß annimmt, ist ihre Entstehungsgeschichte verwinkelt, zunächst uneindeutig 
und bisweilen auch hitzig. Dabei traten einige Hürden auf, die in den Kapiteln dieser Arbeit 
unerwähnt bleiben. Die Überwindung dieser Hürden wurde erst durch meine Unterstützer und 
Mitstreiter möglich, denen ich an dieser Stelle danken möchte. 


Herrn Professor Hans-Jörg Bauer danke ich dafür, am ITS arbeiten zu dürfen und für die Betreu- 
ung meiner Arbeit. Wie auch zahlreichen Vorworten dieser Schriftenreihe zu entnehmen ist, stellt 
die Einbindung in die vielfältigen akademischen Aktivitäten des Instituts eine hervorragende 
Möglichkeit dar, sich fachlich weiterzuentwickeln. Eine aus meiner Sicht besonders wertvolle 
Lehre entstammt jedoch der Arbeit an industrienahen Forschungsprojekten des ITS. In diesen 
Projekten wurden mir große Freiheiten zuteil, sowohl in der Organisation als auch in der wis- 
senschaftlichen und technischen Herangehensweise. Die daraus hervorgehende Verantwortung 
bringt eine professionelle Weiterentwicklung mit sich, von der ich heute besonders profitiere. 
In Anbetracht dessen, dass es sich um Industrieprojekte mit entsprechenden Ansprüchen der 
Auftraggeber handelt, bin ich Herrn Professor Bauer für das mir entgegengebrachte Vertrauen 
besonders dankbar. 


In ähnlicher Weise möchte ich Frau Dr. Corina Schwitzke, in deren Forschungsgruppe ich 
arbeiten durfte, danken. Neben ausgezeichneter fachlicher Kompetenz besitzt Frau Dr. Schwitzke 
selbst in stressigen, schwierigen Situationen die Gelassenheit einer Jazzband und gleichzeitig 
die Organisiertheit eines Raumfahrtprogramms. Diese einzigartige Qualität strahlt auf alle 
Beteiligten aus, wodurch die Arbeit am Institut nicht nur gut funktioniert, sondern auch große 
Freude bereitet. 


Frau Professorin Andrea Luke danke ich für die Übernahme des Korreferats und das entgegen- 
gebrachte Interesse an meiner Arbeit. Herrn OStD Peter Labouvie danke ich für das Lektorat 
meiner Arbeit und die interessanten Diskussionen zu wissenschaftlicher Sprache. 


Experimentelle Ergebnisse stellen einen wesentlichen Teil dieser Arbeit dar. Die Versuche 
hinter diesen Ergebnissen waren stark von Murphys Gesetz betroffen - was schief gehen kann, 
wird schief gehen. Im wiederholten Aufbau und Betrieb, der Fehlersuche und eben häufig 
auch Reparatur der Prüfstände konnte ich glücklicherweise auf das breit aufgestellte Team 
der technischen Mitarbeiter des ITS bauen. Ihre große Fähigkeit und Erfahrung haben sich 
in Verbindung mit der badischen Herzlichkeit dieses Teams als unschlagbares Gegenmittel zu 
Murphys Gesetz erwiesen. Für diese Unterstützung gilt ihnen mein besonderer Dank. 


Der wissenschaftliche Austausch mit den Doktoranden des Instituts führte regelmäßig zu neuen 
Ideen und verbesserten Lösungen für Problemstellungen. Aber auch für mein Leben in Karlsruhe 
war dieses Umfeld ein großer Zugewinn, wofür ihnen mein herzlicher Dank gebührt. Hervorheben 
möchte ich einerseits Herrn Dr. Maximilian Elfner, mit dem ich über lange Zeit ein Büro 
teilte. Seine kreativen Ideen und sein herausragendes Umsetzungsvermögen haben meiner 
Arbeit zahlreiche wichtige Impulse gegeben. Andererseits Herrn Dr. Franz Pütz, der unter den 
Doktoranden stets neue (gerne ökonomisch-wirtschaftliche) Blickwinkel eingebracht und mein 
Leben in mehrfacher Hinsicht bereichert hat. 


Besonders hervorheben möchte ich meine Familie, die meinen Weg weit über das hier Darstell- 
bare hinaus geprägt hat. Meine Lebensgefährtin Constanze, für unzählige fruchtvolle Diskurse, 
Impulse und Erdungen in einem bunten Spektrum an Themen. Meine Eltern, für die vielen 
richtigen Anregungen und die beständige Unterstützung in meinen Entscheidungen. Meine 
Schwester Caroline, die in der gesamten Bandbreite von entscheidendsten Problemen bis zu 
amüsanten Nichtigkeiten eine Wellenlänge mit mir teilt. 


München, im Dezember 2021 Felix von Plehwe 
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Symbolverzeichnis 


1 Einführung 


Eine der stärksten Triebfedern der Forschung im Bereich der Luftfahrt im Allgemeinen und der 
Turbomaschinen im Speziellen ist zweifellos das Streben nach besserer Umweltverträglichkeit 
des Luftverkehrs. Aus Gründen des Wettbewerbs und durch politisch gesetzte Anreize, wie Emis- 
sionsgrenzwerte, herrscht eine hohe Nachfrage nach effizienten, emissionsarmen Triebwerken. 
Die Beherrschung von Technologien, die Treibstoffeinsparungen ermöglichen, bedeutet daher 
einen direkten Wettbewerbsvorteil für einen Triebwerkshersteller. 


Steigerungen des Gesamtwirkungsgrads von Flugzeugtriebwerken wurden in der Vergangen- 
heit durch höhere Wirkungsgrade des Kerntriebwerks erreicht, die wiederum durch erhöhte 
Turbineneintrittstemperaturen und Druckverhältnisse ermöglicht wurden. Diese Fortschritte 
wurden zum Teil durch verbesserte Materialien für Turbinenschaufeln möglich. Seit den 1960er 
Jahren leistet die Kühlung von Turbinenschaufeln durch interne Kühlung und Filmkühlung einen 
erheblichen Beitrag zur Steigerung des thermischen Wirkungsgrads von Triebwerken. Aufgrund 
des hohen Wissensstands in diesen Forschungsbereichen sind weitere Fortschritte nur durch 
enormen Aufwand zu erreichen. 


Neben dem thermischen Wirkungsgrad ist auch der Vortriebswirkungsgrad seit Einführung der 
ersten zivilen Strahltriebwerke kontinuierlich gestiegen. Turbofantriebwerke des 21. Jahrhun- 
derts weisen Nebenstromverhältnisse von bis zu 10 auf und sind dadurch wesentlich effizienter 
als die ersten in der Luftfahrt eingesetzten Strahltriebwerke (Jane’s Information Group (2005)). 


Bei weiterer Erhöhung des Nebenstromverhältnisses gerät das Konzept des konventionellen 
Mantelstromtriebwerks an seine Grenzen. Höhere Nebenstromverhältnisse erfordern größere 
Fandurchmesser. Um transsonische Verluste an den Fanspitzen zu vermeiden und die mechani- 
sche Belastbarkeit des Fans zu gewährleisten, ist die Drehzahl der Niederdruckwelle begrenzt. 
Dadurch erfolgt der Enthalpieumsatz in der Turbine bei geringen Umfangsgeschwindigkeiten 
und geringer Stufenarbeit. Zum Antrieb des Fans können bis zu siebenstufige Niederdrucktur- 
binen nötig sein. Höheres Gewicht, eine größere Anzahl an Bauteilen und höhere Kosten des 
Triebwerks sind die Folge hoher Stufenzahlen. 


Durch ein Untersetzungsgetriebe zwischen Niederdruckturbine und Fan können beide Kom- 
ponenten bei jeweils optimalen Drehzahlen betrieben werden. Bei einer weiteren Steigerung 
des Nebenstromverhältnisses können so höhere Vortriebswirkungsgrade bei einer geringeren 
Stufenanzahl der Niederdruckturbine realisiert werden. Das im Jahr 2016 eingeführte PW 1000G, 
das unter anderem im Airbus A320neo eingesetzt wird, weist Nebenstromverhältnisse von bis zu 
12,5 bei einer dreistufigen Niederdruckturbine auf. Dies führt zu Treibstoffeinsparungen von bis 
zu 16% im Vergleich zu konventionellen Turbofantriebwerken (Pratt & Whitney (2018)). 


Dem Gewichtsvorteil durch eine geringere Anzahl Turbinenstufen steht das Gewicht des Getrie- 
bes gegenüber. Die übertragene Leistung in Triebwerken für sogenannte Single-Aisle-Flugzeuge 
liegt im Bereich von 20 MW, die Wellenleistung von Langstreckenflugzeugen kann über 50 MW 
betragen. Die Übertragung der Leistung durch ein Getriebe bei derart hohen Leistungsdich- 
ten stellt eine enorme technische Herausforderung dar. Zusätzlich ist der verfügbare Bauraum 
innerhalb des Triebwerks begrenzt. 


2 Einfiihrung 


Neben festigkeits- und schwingungsmechanischen Fragestellungen tritt mit der gesteigerten 
Leistungsdichte die Kühlung der Zahnräder und des gesamten Getriebes in den Vordergrund. 
Eine effiziente und effektive Kühlung mittels Öl ist erforderlich, um die thermische Belastung 
des Getriebes zu limitieren. Gleichzeitig wird aus Gründen der Wirtschaftlichkeit angestrebt, die 
Verluste im Getriebe zu minimieren, was einer Forderung nach möglichst geringen Ölvolumen- 
strömen gleichkommt. Ein Optimum zwischen Schmierung, Kühlung und minimalen Verlusten 
kann nur dann erreicht werden, wenn der Wärmeübergang am Zahnrad durch Einspritzkühlung 
genau bekannt und verstanden ist. Dadurch leistet ein verbessertes Verständnis der Öleinspritz- 
kühlung in Getrieben in mehrfacher Hinsicht einen Beitrag zur Senkung der Emissionen und 
Steigerung der Umweltverträglichkeit von Flugzeugtriebwerken. 


In Untersuchungen zum Thermalhaushalt von Getrieben wird häufig die analytische Korrelation 
von DeWinter und Blok (1974) aufgegriffen, um den Wärmeübergang durch Einspritzschmie- 
rung zu bestimmen. Spätere experimentelle Untersuchungen haben jedoch gezeigt, dass die 
gemessenen Temperaturen und abgeschätzten Wärmeübergänge dadurch nicht immer zu erklären 
sind. Experimentelle Untersuchungen, in denen der Wärmeübergang an rotierenden Zahnrädern 
getrennt von Reibwärmeströmen und ortsaufgelöst bestimmt wird, existieren derzeit nicht. 


Diese Lücke soll durch die vorliegende Arbeit geschlossen werden. Es wird eine experimentelle 
Methode und ihre Umsetzung in einem Prüfstand vorgestellt, mit der erstmalig der Wärmeüber- 
gang an einem einspritzgekühlten Zahnrad bestimmt werden kann. 


Folgendes Vorgehen wurde gewählt, um einen Forschungsbeitrag zum Getriebefantriebwerk zu 
leisten: Zunächst wird in Kapitel 2 der derzeitige Stand der Wissenschaft zum Thermalhaushalt 
von hochbelasteten Getrieben geschildert. Darin wird deutlich, weshalb die Bestimmung des 
Öl-Zahnrad-Wärmeübergangs zur Entwicklung effizienter, hochbelasteter Getriebe für Flugtrieb- 
werke unerlässlich ist. Aus dieser Erkenntnis wird die Zielsetzung dieser Arbeit abgeleitet. Da 
gegenwärtig keine Arbeiten bekannt sind, die das System der Einspritzkühlung klassifizieren, 
wird diese in Kapitel 3 zunächst beschrieben. Die zur Beschreibung des Einspritzvorgangs 
verwendeten Größen werden durch die Analyse des Systems abgeleitet. In Kapitel 4 wird die 
entwickelte Methode zur Bestimmung des Öl-Zahnrad-Wärmeübergangs vorgestellt und durch 
die Monte-Carlo-Methode hinsichtlich ihrer Messunsicherheit untersucht. Der zur experimen- 
tellen Umsetzung der Messmethode konzipierte Prüfstand wird in Kapitel 5 beschrieben. Die 
Diskussion der experimentellen Ergebnisse zum Öl-Zahnrad-Wärmeübergang, die durch diese 
Arbeit erstmalig quantitativ bestimmt werden, erfolgt in Kapitel 6. Die Einflüsse der Betriebspa- 
rameter werden diskutiert, anschließend wird eine Korrelation abgeleitet. In Kapitel 7 werden 
die Ergebnisse dieser Arbeit zusammengefasst. 


2 Grundlagen und wissenschaftlicher Kenntnisstand 


Getriebe und Zahnräder sind ein elementarer Bestandteil der industrialisierten Welt. Dement- 
sprechend umfangreich ist die Literatur zur Getriebetechnik. Die Arbeiten zum Thermalhaushalt 
von hochbelasteten Getrieben sind jedoch wesentlich weniger zahlreich. Besonders der Ol- 
Zahnrad-Wärmeübergang durch Einspritzkühlung bedarf eines genaueren Verständnisses. Diese 
Notwendigkeit soll in diesem Kapitel herausgearbeitet werden. 


Die Vorgänge und Phänomene, die sich in einem Getriebe abspielen, sind vielfältig, komplex und 
umspannen ein breites naturwissenschaftliches Feld. Für einfache Konstruktionen, wie sie die 
Mehrzahl der global eingesetzten Getriebe darstellen, ist die genaue Kenntnis dieser Vorgänge 
erlässlich und lässt sich durch Auslegung hin zu großen Sicherheitsvielfachen ersetzen. Mit 
steigenden Anforderungen an das Getriebe steigt die Notwendigkeit, die beteiligten Vorgänge 
genau zu verstehen. Wie für andere Maschinenelemente auch stellen die Luftfahrt und der 
Turbomaschinenbau im Besonderen extreme Anforderungen an Getriebe. Daraus entstehen 
erhöhte Anforderungen an das Detailverständnis der beteiligten physikalischen Vorgänge. 


Die Herausforderung, ein Höchstleistungsgetriebe für ein Luftfahrttriebwerk zu entwickeln, 
zeichnet sich aus durch die Forderungen nach 


e sehr hoher übertragener Leistung 
e geringem Gewicht und geringem Bauraum 
e hoher Zuverlässigkeit 


e hohem Wirkungsgrad. 


Die Erfüllung dieser Forderungen ist dadurch erschwert, dass sie stark gegensätzlich sind. Die 
Forderung nach hoher Leistungsdichte bei geringem Gewicht und Bauraum bedeutet unweiger- 
lich hohe spezifische Wärmeströme durch Reibung und damit hohe lokale Materialtemperaturen, 
wie in Kapitel 2.1.1 erläutert wird. Über die Elastohydrodynamik, die in Kapitel 2.1.2 umrissen 
wird, lässt sich zeigen, dass hohe Materialtemperaturen für eine gegebene Getriebekonfiguration 
zu geringeren Schmierfilmdicken führen. Abnehmende Schmierfilmdicken haben Mischreibung 
im Eingriff zur Folge, bei der der Schmierfilm örtlich unterbrochen und die erwartete Lebens- 
dauer der Zahnräder herabgesetzt wird. Dies steht der Forderung nach hoher Zuverlässigkeit 
entgegen. Die Begegnung dieses Problems durch den verstärkten Einsatz von Ol, dessen Haupt- 
funktion die Kühlung der Zahnräder ist, bringt jedoch Einbußen im Wirkungsgrad mit sich, da 
Ventilations- und Planschverluste in Getrieben mit der verwendeten Ölmenge ansteigen. Der 
derzeitige wissenschaftliche Kenntnisstand hält keine Hinweise darauf bereit, wie die Kühlwir- 
kung der eingesetzten Olmenge maximiert werden kann. Die in Abschnitt 2.2 beschriebenen 
Ansätze zur Thermalmodellierung erlauben zwar die Modellierung der sich einstellenden Zahn- 
flankentemperaturen. Die benötigten Randbedingungen zum Zahnflankenwärmeübergang bei 
Einspritzkühlung sind jedoch nur als grobe Abschätzung bekannt, wie in Abschnitt 2.3 verdeut- 
licht wird. Die Bestimmung der Randbedingungen aus allgemeineren Arbeiten zur Prallkühlung 


4 Grundlagen und wissenschaftlicher Kenntnisstand 


ist nur überschlägig möglich, wie in Abschnitt 2.4 dargelegt wird. Wird also ein Optimum 
zwischen den gegensätzlichen Anforderungen gesucht, müssen zunächst die gekoppelten und 
komplexen Vorgänge im Getriebe verstanden werden. Erst dann können die oben genannten 
Forderungen bestmöglich erfüllt werden. 


2.1 Getriebemechanik 


Das Zahnrad als Element zur Übertragung von Leistung und zur Übersetzung von Drehzahl 
und Drehmoment gilt als eines der ältesten Maschinenelemente der Welt. Die erste eingehende 
mathematische Beschreibung von Zahnformen geht auf das Jahr 1694 und Philippe de la Hire 
zurück, der darin die epizykloide Zahnform beschreibt. Die heute gebräuchlichste Form, die 
Evolventenverzahnung (siehe Abbildung 2.1), wurde von Leonard Euler 1760 erstmals beschrie- 
ben (Klocke und Brecher (2017)). Ihr Vorteil besteht in der gleichmäßigen und damit weitgehend 
vibrationsfreien Übertragung der Drehbewegung. Diese Eigenschaft ist für den verschleiß- 
und geräuscharmen Betrieb von großer Bedeutung. Sie lässt sich mathematisch als Forderung 
nach = = const. für die Rotationsgeschwindigkeit des treibenden und angetriebenen Zahnrads, 
@, und @», zu jedem Zeitpunkt während des Eingriffs ausdrücken. Sie ist dann erfüllt, wenn 
die Berührnormale des Eingriffspunkts den Wälzpunkt schneidet. Dieser Zusammenhang ist 
als Verzahnungsgesetz bekannt, das außer bei der erwähnten Evolventenverzahnung auch bei 
Zykloiden-, Triebstock-, Kreisbogen- und Wildhaber-Novikov-Verzahnungen erfüllt ist (Klocke 
und Brecher (2017)). Aufgrund der herausragenden Bedeutung der Evolventenverzahnung im 
Luftfahrtumfeld und für den Kontext des Höchstleistungsgetriebes soll diese in dieser Arbeit 
ausschließlich betrachtet werden. 


2.1.1 Reibung und Verluste 


In Abbildung 2.1 sind zwei im Eingriff befindliche Zahnräder mit Evolventenverzahnung schema- 
tisch dargestellt. Der Punkt A stellt dabei den Beginn, Punkt E das Ende eines Zahnradeingriffs 
dar. Die Summe aller Berührpunkte zwischen den Punkten A und E wird Eingriffstrecke genannt. 
Bei der dargestellten Evolventenverzahnung ist diese Strecke eine Gerade, die sogenannte Ein- 
griffsgerade AE. Im Berührpunkt des treibenden und des angetriebenen Zahnrads liegt dabei, mit 
Ausnahme des Wälzpunktes, stets eine Relativgeschwindigkeit der Oberflächen vor. Dies wird 
durch die Darstellung der Gleitgeschwindigkeiten des treibenden und des getriebenen Zahnrads, 
Vg1 und vg>, deutlich. Am Wälzpunkt C erfahren sie einen Vorzeichenwechsel. Der Verlauf der 
Geschwindigkeiten ist linear über die Eingriffstrecke. Diese Gleitbewegung ist ursächlich für die 
entstehende Reibwärme im Getriebe, die einen wesentlichen Anteil an den Gesamtverlusten hat. 


Die Verlustleistung kann bei gering belasteten Getrieben allein durch Wärmeleitung über die 
Zahnradkörper abgeführt werden, wodurch die Materialtemperaturen in einem unkritischen 
Bereich bleiben. Mit steigender Last können jedoch kritische Temperaturwerte erreicht werden, 
die zu einer Lebensdauerreduzierung oder gar zum unmittelbaren Versagen des Getriebes füh- 
ren. Um für erhöhte Lasten Dauerfestigkeit sicherzustellen, werden Getriebe daher mit einem 
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Abbildung 2.1: Schematische Darstellung des Eingriffs zweier evolventischer Zahnräder und der 
auftretenden Geschwindigkeiten nach Klocke und Brecher (2017) 


Schmiermittel versehen. Zum einen dient dieses dazu, den Reibwert u und damit die Verlust- 
bzw. Reibleistung Py, herabzusetzen. Zum anderen kann durch das Schmiermittel auch Wärme 
abgeführt werden. In einfachen Konstruktionen gibt das Schmiermittel die aufgenommene 
Wärme an das Gehäuse ab (Sumpfschmierung), in aufwendigeren Konstruktionen können dafür 
Ölkühler eingesetzt werden. Mit steigender Leistungsdichte tritt der Aspekt der Kühlung durch 
das Öl zunehmend in den Vordergrund. Das gilt auch für andere tribologische Maschinenele- 
mente wie beispielsweise Wälzlager, deren Schmierung und Kühlung durch Öl am Institut für 
Thermische Strömungsmaschinen seit mehreren Jahrzehnten Gegenstand intensiver Forschung 
ist (Glahn et al. (1997), Gorse et al. (2003)). 


Die Kühlung und Schmierung von Wälzlagern in Triebwerken erfolgt, wie im vorliegenden Fall, 
häufig durch Einspritzkühlung (Krug (2018)). Zur Kühlung der Wälzlager ist eine ausreichende 
Zufuhr von Öl sicherzustellen. Das Öl interagiert, nachdem es die Einspritzdüse verlässt, mit 
der von der Welle mitbewegten Luft in der Lagerkammer und kann dadurch zerstäubt werden. 
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Dadurch erreicht nur ein Anteil des eingespritzten Olvolumenstroms das Lager, der durch den 
Auffangwirkungsgrad beschrieben wird. (Krug et al. (2014)) zeigen, dass die Wellendrehzahl 
und der eingespritzte Olvolumenstrom entscheidend fiir den Auffangwirkungsgrad sind. Das 
eingespritzte Ol nimmt anschließend an den Strömungsvorgängen der Lagerkammer teil und 
beeinflusst diese. Kurz und Bauer (2014) identifizieren zwei charakteristische Strömungsregime, 
die vor allem vom Lagerkammerdurchmesser und dem Ölvolumenstrom durch die Lagerkammer 
abhängen. Die Strömung in der Lagerkammer, die sich durch die Geometrie und die Betriebs- 
bedingungen einstellt, ist entscheidend für den Wärmeübergang an den Lagerkammerwänden. 
Sowohl Glahn (1995) als auch Kurz (2018) charakterisieren diesen Zusammenhang experimen- 
tell. Die Kenntnis der komplexen Vorgänge dient im Kontext der Lagerkammer zur Steigerung 
der umsetzbaren Leistungsdichte. 


Durch den Einsatz von Öl als Schmier- und Kühlmittel können auch in einem Getriebe höhere 
Leistungsdichten realisiert werden. Gleichzeitig treten dadurch zusätzliche Verluste auf, da das 
Öl zwangsläufig an den komplexen Strömungsvorgängen beteiligt ist. Eine häufig getroffene 
Aufteilung der Getriebeverluste ist die von Niemann und Winter (2003): 


Zahnräder Lager 
i aa 
Py = Pyz + Py20 + Pvp + Pvpo + Pvp + Pvx (2.1) 


Darin sind Pyz und Pyzo die Verzahnungsverluste, Pyg und Pygo die Lagerverluste, Pyp die 
Dichtungsverluste und Pyx sonstige Verluste. In Gleichung 2.1 kennzeichnet der Index 0 diejeni- 
gen Verluste, die lastunabhängig sind. Diese lastunabhängigen Verluste stellen sich unabhängig 
von der übertragenen Leistung bzw. dem übertragenen Moment ein und beinhalten Plansch-, 
Ventilations- und Verdrängungsverluste. Während für langsam drehende Getriebe und hohe Dreh- 
momente die lastabhängigen Verluste dominieren, überwiegen bei Hochdrehzahlanwendungen 
die lastunabhängigen Verluste (Diab et al. (2006a)). 


Bei der Sumpfschmierung, die auch Tauchschmierung genannt wird, tauchen die Zahnräder in 
ein Olbad ein und werden allein durch ihre Eigenbewegung mit Öl versorgt und gekühlt. Die 
Sumpfschmierung stellt die einfachste konstruktive Methode der Ölversorgung dar. Diesem 
Vorteil stehen Verluste entgegen, die durch das vergleichsweise hohe Ölvolumen entstehen, das 
an den Strömungsvorgängen im Getriebe beteiligt ist. Durch die Reduzierung der Füllstandshöhe 
des Öls können diese Verluste gesenkt werden (Seetharaman et al. (2009), Michaelis et al. (2011)). 
Bei zu geringen Füllstandshöhen ist jedoch sowohl mit mangelnder Benetzung der Zahnräder mit 
Öl als auch mit unzureichender Kühlung zu rechnen, was zur Beschädigung der Zahnräder führen 
kann (Höhn et al. (2008)). Schober (1983b) nennt eine Grenze der Umfangsgeschwindigkeit von 
Zahnrädern von etwa 60 m/s, bis zu der Sumpfschmierung sinnvoll eingesetzt werden kann. Bei 
Planetengetrieben, wie sie für den Einsatz in Höchstleistungsgetrieben für Getriebefantriebwerke 
in Frage kommen, scheidet der Einsatz von Sumpfschmierung aus. 


Durch Einspritzen des Öls können Schmierung und Kühlung der Zahnräder bei höheren Um- 
fangsgeschwindigkeiten ermöglicht werden, als es mit Sumpfschmierung möglich ist. Dem 
größeren konstruktiven Aufwand stehen dabei geringere Verluste gegenüber (Schober (1983a), 
Andersson et al. (2017)). In einspritzgeschmierten Getrieben sind die lastunabhängigen Verluste 
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vor allem von der Drehzahl, der Geometrie der Zahnräder und der unmittelbaren Umgebung 
sowie der Dichte des umgebenden Fluids abhängig (Eastwick und Johnson (2008)). 


Der Zusammenhang zwischen der Dichte des Fluids im Getriebe und lastunabhängigem Verlust 
bedeutet, dass eine zunehmende Konzentration von Öl in Form von feinen Tropfen eine Erhöhung 
der Verluste mit sich bringt. Erhöhte Konzentrationen sind dann zu erwarten, wenn mehr 
Öl am Schmier- und Kühlprozess teilnimmt. Zusätzlich hat der Impulsaustausch zwischen 
eingespritztem Öl und Zahnrad unmittelbar Verluste zur Folge, die mit zunehmender Öl€menge 
steigen (Massini et al. (2017)). Eine geringere Menge eingespritzten Öls führt daher zu geringeren 
Verlusten. Darüber hinaus können so die benötigten Pumpen und Übertragungssysteme kleiner 
dimensioniert werden. Dies führt zu einem geringeren benötigten Bauraum der Systeme und zu 
Gewichtseinsparungen. 


Für fortschrittliche, hochbelastete Getriebe mit hohen Drehzahlen sind Wirkungsgrade zwischen 
98% und 99% erreichbar (Eastwick und Johnson (2008)). Im Fall des Hilfsgerätegetriebes von 
Luftfahrttriebwerken können die Strömungsverluste dabei mehr als 50% des Gesamtverlustes 
ausmachen (Arisawa et al. (2014)). 


Fiir den Thermalhaushalt des Zahnrads, der in Abschnitt 2.2 eingehend betrachtet wird, ist 
besonders die Verlustleistung im Eingriff, die durch Py, erfasst wird, von Bedeutung. Diese 
wird unmittelbar an den belasteten Zahnflanken und somit auf kleiner Fläche freigesetzt, was 
zu starker Erwärmung der Zahnflanken führt. Die lastunabhängigen Verluste entstehen durch 
Dissipationsvorgänge, die sich im gesamten vom Getriebe eingeschlossenen Volumen abspielen. 
Sie führen zu einer Erwärmung des Kühlmittels, sind jedoch für die Zahnradtemperaturen nur 
von mittelbarer Bedeutung. 


Die durch die Gleitbewegung entstehende Verlustleistung Py, kann ausgedrückt werden als: 


Pyz = Frvg = (1/pe) [Fre dg (2.2) 
& 


mit der Reibkraft Fr, der Gleitgeschwindigkeit vg, der Eingriffsteilung pe, der Normalkraft Fn und 
der Eingriffstrecke g. Der Reibungsbeiwert u ist eine örtliche Größe, die unter anderem von der 
Materialpaarung und der Art der Schmierung abhängig ist. In Gleichung 2.2 ist nur der Reibungs- 
beiwert u nicht von der Konfiguration des Getriebes abhängig. Neben ihrer Bedeutung für den 
Thermalhaushalt des Getriebes ist die Verlustleistung Py, auch für die Wirkungsgradbetrachtung 
in Getrieben wichtig. 


Die Bestimmung des Reibungsbeiwerts u kann experimentell erfolgen. Mittlere Werte für u 
können jedoch auch durch empirische Näherungsformeln wie der von Winter und Michaelis 
(1984) bestimmt werden: 


Kıkn/b 


VEmPm 


0,2 
Um = 0,045 ( ) 02. (2.3) 


Darin ist Ka ein Anwendungsfaktor, der von der Art der Belastung abhängt, b die Breite des 
Zahnrads, vy die mittlere Summengeschwindigkeit und pm der mittlere Ersatzkrümmungs- 
radius. Die Summengeschwindigkeit ist die Summe der Tangentialgeschwindigkeiten (siehe 
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Abbildung 2.1) und beeinflusst den Aufbau des Schmierfilms maßgeblich. no ist die Viskosität 
des Öls bei Massentemperatur, also der Temperatur der Zahnradoberfläche. Xp ist ein Faktor für 
die Rauigkeit, mit 

Xp = 3,8(Ra/dw)?. (2.4) 


Zur Berechnung des instantanen Reibungsbeiwerts an einem beliebigen Punkt der Eingriffsgerade 
werden in der Literatur häufig die Korrelation von Kelley und Lemanski (1967) 


1 50, 8F,/b | 
= 0,0099 —— — 10 25 
H IRB es: G 
oder die Korrelation von Benedict und Kelley (1961) 
3,17: 10°F, /b 
u = 0,0127 -10810 (e) (2.6) 
NOV VS 


verwendet. Gleichungen 2.5 und 2.6 liefern Reibungsbeiwerte, die im Allgemeinen mit den Wer- 
ten experimenteller Untersuchungen übereinstimmen. Allen genannten Bestimmungsansätzen 
für den Reibwert ist gemein, dass der Reibwert bei höheren Oberflächen- und Gleitgeschwindig- 
keiten abnimmt, während höhere bezogene Normalkräfte zu höheren Reibwerten führen. 


Diab et al. (2006b) zeigen, dass die genannten Ansätze den Reibwert aufgrund der geringen 
relativen Gleitgeschwindigkeiten an Zahnrädern tendenziell überschätzen. Effekte wie die Ani- 
sotropie und die Orientierung der Oberflächenrauigkeit gefertigter Zahnräder führen zusätzlich 
zu Abweichungen des tatsächlichen Reibwerts im Vergleich zu Korrelationen. 


Wesentlich genauere Vorhersagen des Reibwerts könnten getroffen werden, indem die elastohy- 
drodynamischen Vorgänge vollständig modelliert werden. Britton et al. (2000) zeigen anhand 
derartiger Berechnungen und experimenteller Untersuchungen unter identischen Bedingungen 
sehr geringe Abweichungen der Werte. 


Festzuhalten ist, dass der Reibwert u von der Konfiguration des Getriebes, dem Lastzustand und 
dem thermischen Zustand der Zahnräder abhängt. Hohe Normalkräfte im Eingriff können dabei 
den Reibwert erhöhen, jedoch ist die Abhängigkeit nach Gleichung 2.3 mit 


Um x FO? (2.7) 
oder nach Gleichungen 2.5 und 2.6 mit 


H œ logio (Fn) (2.8) 


schwach. Wird die Leistung in einem Getriebe auf kleinerem Bauraum umgesetzt, ist nur mit 
geringen Steigerungen des Reibwerts und somit der lastabhängigen Verluste zu rechnen. Die 
Fläche, auf der die Reibleistung in Form von Wärme freigesetzt wird, sinkt jedoch. Dadurch 
ist in Höchstleistungsgetrieben, wie sie im Rahmen dieser Arbeit betrachtet werden, mit ho- 
hen Wärmestromdichten durch Reibung zu rechnen. Der Kühlung der Zahnräder bei hohen 
Umfangsgeschwindigkeiten kommt somit besondere Bedeutung zu. 
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2.1.2 Tribologie 


Schmierstoffe wurden bereits eingesetzt, um Verschleiß und Verluste in Reibkontakten zu redu- 
zieren, bevor die zu Grunde liegenden Mechanismen der Schmierfilmbildung und der Kontakt- 
mechanik bekannt oder verstanden waren (Dowson et al. (1977)). Die erstmalige Beschreibung 
dieser tribologischen Mechanismen erfolgt durch Reynolds (1886). Reynolds schlägt einen ma- 
thematischen Lösungsansatz für die von Tower (1883) beobachteten, äußerst geringen Reibungen 
zwischen vollständig in Öl eingetauchten Oberflächen vor. Dieser beruht auf der Annahme eines 
durchgehenden Schmierfilms zwischen den Oberflächen, einem wesentlichen Gedanken der 
Tribologie. Aufbauend darauf vereinfacht Reynolds die hydrodynamischen Grundgleichungen, 
um zu einem Ausdruck für die sich einstellende Schmierfilmdicke zu gelangen. Die Auswer- 
tung dieses Ausdrucks ergibt für die real umgesetzten hohen Normallasten stets sehr kleine 
Schmierfilmdicken, die wesentlich kleiner als die herstellbaren Oberflächenrauigkeiten sind. Die 
klassische Hydrodynamik scheitert demnach bei der Erklärung eines Schmierfilms für Zahnräder 
unter den technisch bereits umgesetzten Lasten und Kontaktdrücken. 


Erst die Berücksichtigung zweier weiterer wesentlicher Aspekte liefert plausible Schmierfilmdi- 
cken: 


e Die Viskosität typischer Schmieröle steigt unter hohem Druck erheblich an. Der Faktor 
liegt üblicherweise bei 10*. Dies ist als Piezoviskoser Effekt bekannt. 


e Die Kontaktflächen der Zahnräder verformen sich im Eingriff elastisch. 


Die Betrachtung des Schmierfilms unter Berücksichtigung der beiden genannten Effekte wird 
in der Literatur als Elastohydrodynamik (EHD) bezeichnet (Olver (2002)) und wurde in den 
1950er und -60er Jahren etabliert (Grubin und Vinogradova (1949), Barwell (1958), Dowson et al. 
(1977)). Erst durch die Theorie der Elastohydrodynamik lassen sich Filmdicken erklären, die 
mit experimentell ermittelten Werten in Einklang gebracht werden können. Die Erkenntnisse der 
1960er Jahre, die von Dowson und Higginson zusammengefasst werden (Dowson et al. (1977)), 
werden im Bereich der Getriebemechanik aufgegriffen und erweitert, um auch Schmierfilme 
zwischen gleitenden Oberflächen zu beschreiben. Eine solche Gleitbewegung liegt bei der 
Evolventenverzahnung zwangsläufig vor. 


In Abbildung 2.2 ist der qualitative Druck verlauf des Schmierfilms im Eingriff zweier Wälz- 
körper sowohl für die Hertz’sche Theorie als auch für die EHD-Theorie dargestellt. Über die 
Hertz’sche Theorie kann die Druckverteilung sich berührender Festkörper ohne den Einfluss 
eines Schmierfilms bestimmt werden (Hertz (1881)). In weiten Bereichen stimmt die Druck- 
verteilung beider Theorien überein. Ein wichtiges Merkmal des Verlaufs der Schmierfilmdicke, 
die durch die EHD-Theorie bestimmt wird, ist die Einschnürung am Eingriffsauslass. An der 
Stelle der Einschnürung ist ein Druckmaximum zu erkennen, das charakteristisch für EHD- 
Schmierfilme ist und Petrusevich’sche Spitze genannt wird. Diese Einschnürung wird durch die 
Theorie vollständig erklärt (Dowson et al. (1977)) und ist experimentell bestätigt. Ihr kommt 
besondere Bedeutung zu, da sie die Stelle dünnster Schmierfilmdicke darstellt und somit für 
Verschleißbetrachtungen essentiell ist. 


10 Grundlagen und wissenschaftlicher Kenntnisstand 


Pe Petrusevich’ sche 
Spitze 


Druckverlauf 
nach Hertz Se 
LA 


+‘ 


EHD- 
Druckverlauf 


EHD- 
Schmierfilm 


Abbildung 2.2: Druckverteilung im Schmierspalt nach Hertzscher und elastohydrodynamischer 
Theorie und charakteristischer EHD-Schmierfilm nach Theißen (2003) 


Die exakte Vorhersage von Schmierfilmdicken für unterschiedliche tribologische Systeme wie 
Gleitlager oder Zahnräder ist aufgrund der Vielzahl beteiligter Vorgänge schwierig. Die De- 
formation der Festkörper und die Fluidmechanik des Schmierfilms sind auch für einfachste 
Betrachtungen zu berücksichtigen, um korrekte Ergebnisse zu erhalten. Je nach Anwendungsfall 
können weitere physikalische Vorgänge, wie beispielsweise Wärmetransportvorgänge (Pinkus 
(1990)) in oder aus dem Schmierfilm, einen Einfluss auf die Schmierfilmdicke haben. Betrach- 
tung und Bewertung des Einzelfalls sind umfangreich und sollen im Rahmen dieser Arbeit nicht 
angestellt werden. 


Dowson et al. (1977) stellen jedoch fest, dass die Betrachtung der Schmierfilmdicke anhand 
dreier dimensionsloser Kennzahlen G, W und U für alle bekannten theoretischen Lösungen durch 
eine Näherungsgleichung bereits gute Ergebnisse liefert. Dazu wird der Materialparameter 


G= aE’ (2.9) 


definiert. Dieser ist über den Druck-Viskositätskoeffizienten & vom Schmiermittel und den 
Ausdruck E’ = i b vom Wälzkörpermaterial abhängig. Über & wird der eingangs erwähnte, 
starke Viskositätsanstieg der meisten Öle für sehr hohe Drücke berücksichtigt. E’ ist eine 
Funktion des E-Moduls E und der Querkontraktionszahl v. Für eine definierte Kombination aus 


Schmiermittel und Zahnradmaterial ist G eine Konstante. 


Der Lastparameter 
w 


ER 
ist abhängig von der längenbezogenen Normalkraft w = A. F, ist die Kraft, die normal zur 
Berührfläche auf den tribologischen Kontakt bzw. den Eingriff wirkt. / ist die Eingriffslänge, 
über die die Kraft Fp übertragen wird. Für einen gegebenen Lastfall ist auch W eine Konstante. 


(2.10) 
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Der Geschwindigkeitsparameter 
_ Nou 
ER 
ist abhängig vom äquivalenten Wälzkörperradius R, der Oberflächengeschwindigkeit u sowie no. 
Die Viskosität No ist die Viskosität des Schmiermittels bei der Materialtemperatur zu Beginn des 
Eingriffs. Bei gegebenen Betriebszuständen und für ein ausgeführtes Getriebe ist nur No eine 


(2.11) 


Variable. Die minimale Schmierfilmdicke Amin Kann aus den obenstehenden dimensionslosen 
Kennzahlen nach Dowson et al. (1977) durch den Ausdruck 


GPU" 


Amin = 1, ORs 


@Q12) 


bestimmt werden. Durch Gleichung 2.12 ist die fiir die Auslegung von Verzahnungen und Wälz- 
körpereingriffen im Allgemeinen wichtige minimale Schmierfilmdicke, Amin, durch eine sehr 
einfache mathematische Beschreibung zugänglich. Sie ist im Bereich der elastohydrodynami- 
schen Schmierung gültig, die für Bereiche von U und W vorliegt, die nahezu alle real auftretende 
Anwendungsfälle abdecken. Die Kennzahlen der im Rahmen dieser Arbeit betrachteten, hochbe- 
lasteten Getriebe liegen stets im Gültigkeitsbereich von Gleichung 2.12. Durch die Abhängigkeit 
von Amin von W mit einem Exponenten von —0,13 ist die minimale Schmierfilmdicke nur äußerst 
schwach von der Zahnradlast abhängig. 


Die Anwendbarkeit von Gleichung 2.12 ist durch zahlreiche Veröffentlichungen eindeutig 
belegt. Auch in der Normungsliteratur wird sie als Auslegungskriterium zur Sicherheit gegen 
Warmfressen genutzt. Die internationale Norm zur Auslegung von hochbelasteten Getrieben mit 
Sicherheit gegen Graufleckenbildung, ISO TR 15144, wurde im Jahr 2014 veröffentlicht und 
verwendet Gleichung 2.12 als Auslegungskriterium (ISO (2014)). 


Für ein gegebenes Getriebe mit gegebener Lastanforderung und Geometrie sind der Lastparameter 
W und der Materialparameter G konstant. Der Geschwindigkeitsparameter U ist, bei konstantem 
u, ausschließlich eine Funktion von no, der Viskosität des Schmiermittels bei Umgebungsdruck. 
Für die Schmierfilmdicke hmin gilt daher: 


Amin = n9” (To). (2.13) 


Somit ist die temperaturabhängige Viskosität no (To) des Ols unmittelbar vor dem Eingriff der 
Zahnräder entscheidend für die Dicke des Schmierfilms. Dabei ist Tọ nicht etwa die Einspritz- 
temperatur des Öls oder die Temperatur im Ölsumpf, sondern die Temperatur der Zahnräder 
vor dem Eingriff. In der englischsprachigen Literatur wird für Tọ üblicherweise der Begriff bulk 
temperature verwendet. 


Zur Bestimmung der Schmierfilmdicke sind durch die genannten Zusammenhänge die Visko- 
sität des Schmiermittels und ihr Temperaturverhalten entscheidend. Typische Schmieröle im 
Luftfahrtbereich besitzen einen Viskositäts-Temperaturverlauf, der durch die Arrhenius-Andrade- 
Beziehung angenähert werden kann. Grundsätzlich sinkt die Viskosität mit steigender Temperatur. 
Dies gilt auch für das häufig verwendete Mobil Jet Oil II, dessen Temperaturverlauf von Glahn 
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(1995) beschrieben wird. Die EHD-Schmierfilmdicke für Jet Oil II ist in Abbildung 2.3, bezo- 
gen auf die Schmierfilmdicke bei 100°C und fiir einen generischen Lastfall, dargestellt. Eine 
Erhöhung der Temperatur von 100 °C auf 148 °C führt zu einer Halbierung der Schmierfilmdicke. 
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Abbildung 2.3: Einfluss der Oberflächentemperatur vor dem Eingriff auf die minimale relative 
Schmierfilmdicke 


Die Schmierfilmdicke ist entscheidend für mehrere Schadensmechanismen, die die Leistungs- 
dichte und die Lebensdauer eines Getriebes begrenzen. In Abbildung 2.4 sind die möglichen 
Zustände des dünner werdenden Schmierfilms dargestellt. Bei ausreichend großem Amin ist stets 
ein durchgängiger Schmierfilm zwischen den Zahnflanken vorhanden, die Gleitbewegung führt 
lediglich zu einer Scherung des Fluids im Schmierspalt unter Reibwärmeentwicklung (engl.: 
Full-film EHD, Flüssigkeitsreibung). Bei reduzierter Schmierfilmdicke nähern sich lokale Rauig- 
keitsspitzen an, es kommt zunächst zur Ausbildung lokaler, kleinstskaliger Schmierfilme (engl.: 
Micro-EHD) und anschließend zu örtlich begrenzten Unterbrechungen des Schmierfilms (engl.: 
Mixed lubrication, Mischreibung). Ab einer gewissen Konzentration solcher Unterbrechungen 
wird diese Mischreibung kritisch. Durch die Festkörperreibung der beteiligten Flanken kommt 
es zu kurzzeitigen Temperaturspitzen im Metall. Im Falle des sogenannten Fressens erreicht 
das kurzzeitige Temperaturniveau den Schmelzpunkt des Metalls, es kommt zum Verschweißen 
der Berührflächen und unmittelbar darauf folgendem Aufreißen der Verbindung (Niemann und 
Winter (2003)). Ist der Schmierfilm vollständig unterbrochen, spielt das Schmiermittel nur noch 
in kleinstskaligen Vorgängen eine Rolle (engl.: Boundary lubrication). 


Die Oberflächenrauigkeit in Verbindung mit der minimalen Schmierfilmdicke ist entscheidend 
für den Zustand im Schmierspalt. Zur Bewertung der Sicherheit gegen Warmfressen wird die 
spezifische Schmierfilmdicke 


h ; 
Acr = — 2.14 
GF = p, (2.14) 
unter Verwendung der mittleren Rautiefe 
Ra +R 
R, = S tia (2.15) 


2 
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Abbildung 2.4: Schmierfilmzustände im Zahnradeingriff nach Spikes (1997) 


mit den Rauigkeiten der beiden Flanken, R,; und R,>, verwendet. Agr wird typischerweise in 
die in Tabelle 2.1 aufgefiihrten Bereiche aufgeteilt. 


Tabelle 2.1: Bewertung der spezifischen Schmierfilmdicke Agr 


AGE Schmierfilmzustand 
>2 Überwiegend hydrodynamische Schmierung 
0,7—2 Grenzschmierung 
< 0,7 Festkörperreibung 


Der Übergang von reiner Flüssigkeitsreibung zur Festkörperreibung über den Bereich der Misch- 
reibung verläuft kontinuierlich und hängt von einer Vielzahl weiterer Faktoren ab. So spielt auch 
die Gestalt der beteiligten Rauigkeiten eine Rolle (Guegan et al. (2016)). Es wird beobachtet, 
dass bereits ab einem Wert Agr von 2-3 erste Rauigkeitsspitzen in Berührung kommen und bis 
zu einem Wert von 1 der Großteil der Last vom Schmierfilm getragen wird (Spikes (1997)). 


Krantz (2014) untersucht den Einfluss der spezifischen Schmierfilmdicke auf die Lebensdauer 
der beteiligten Zahnräder experimentell. Dazu werden Zahnradpaare in einem Prüfstand unter 
kontinuierlicher Leistungsübertragung so lange betrieben, bis es zu Schäden an den Flanken 
unter Entstehung von Vibrationen und Ausfall der Verzahnung kommt. Durch unterschiedliche 
Kombinationen aus Schmiermittel, Rauigkeit des Zahnrads und Zahnradmaterial werden unter- 
schiedliche spezifische Schmierfilmdicken erzeugt. Die Verteilung der durchlaufenen Zyklen bis 
zum Versagen in Abhängigkeit der spezifischen Schmierfilmdicke ist in Abbildung 2.5 dargestellt. 
Während die Lebensdauer bei Agr > 2 unabhängig von der Schmierfilmdicke ist, nimmt sie bei 
Unterschreiten eines Wertes von etwa 2 stark ab. Bei Schmierfilmen, die mit Agr < 1 dünner 
als die mittlere Rautiefe der Zahnräder sind, werden deutlich reduzierte minimale Werte für die 
Lebensdauer gemessen. 


Grundsätzlich wird das Auftreten von Festkörperreibung bei abnehmenden Schmierfilmdicken 
mit erhöhter Gefahr des Verschleißes in Verbindung gebracht. Bei additivierten Schmierstoffen 
kann bei erhöhten Temperaturen jedoch die chemische Aktivität des Additivs erhöht sein. Die 
dadurch zunehmende Schutzwirkung Kann den negativen Einfluss abnehmender Viskosität und 
zunehmenden Festkörperkontakts hinsichtlich des Verschleißes unter Umständen kompensieren 
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Abbildung 2.5: Einfluss der spezifischen Schmierfilmdicke Agr auf die Lebensdauer nach Krantz 
(2014) 


(Höhn und Michaelis (2004)). In diesem Fall kann jedoch ein Temperaturbereich existieren, in 
dem die Schmierfilmdicke bereits kritisch ist, die chemische Aktivität des Additivs jedoch noch 
nicht ausreichend hoch ist. Zusätzlich kann der additivierte Schmierstoff nicht bei beliebig hohen 
Temperaturen eingesetzt werden. Daher ist gerade bei Einsatz solcher Schmierstoffe die Kenntnis 
der Zahnradtemperaturen wichtig, um schadensanfällige Betriebsbereiche zu vermeiden. 


Neben der Bewertung der Filmdicke hinsichtlich Verschleiß und Lebensdauer ist Agp auch 
für die Entstehung der Reibwärme von Bedeutung. Im Allgemeinen ist der Reibwert u für 
Festkörperreibung wesentlich höher als für Misch- und Flüssigkeitsreibung. 


Die Erkenntnisse der EHD können dazu verwendet werden, die Reibwärmeentwicklung im 
Zahnradeingriff zu berechnen. Die Temperaturen sowohl des Ölfilms als auch des Zahnrads und 
die Reibleistung sind gekoppelt zu betrachten, wie in Abschnitt 2.2 näher beschrieben wird. Es 
ist festzuhalten, dass nach der elastohydrodynamischen Theorie die Temperaturen zu Beginn 
des Eingriffs für die Schmierfilmdicke entscheidend sind (Olver (2002)). Erhöhte Temperaturen 
von Öl und Zahnrad führen zu geringeren Schmierfilmdicken, was bei gleichbleibender Last und 
konstanter Oberflächenrauigkeit der Zahnräder die Wahrscheinlichkeit von Mischreibung erhöht. 
Diese Erkenntnis spiegelt sich auch darin wider, dass die meisten in der Praxis beobachteten Fälle 
des Zahnradfressens durch zu hohe Temperaturen und nicht etwa durch zu geringe Schmierung 
auftreten (Townsend und Akin (1983)). Eine genaue Kenntnis des thermischen Zustands zu 
Beginn des Eingriffs ist demnach unerlässlich, um Schadensprozesse wie Micropitting und 
Warmfressen ausschließen zu können. Im Umkehrschluss können durch die genaue Kenntnis der 
Zahnradtemperaturen höhere Leistungsdichten in Getrieben realisiert werden. 
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2.2 Thermalmodellierung von Zahnrädern 


Die im vorangegangenen Abschnitt geschilderten Arbeiten und Erkenntnisse erlauben die Erklä- 
rung, Vorhersage und Bewertung der auftretenden Schmierfilme. Die Bestimmung der Reibwär- 
me, die durch die tribologischen Vorgänge im Eingriff entsteht, ist durch die in Abschnitt 2.1.1 
erwähnten Arbeiten möglich. Außerdem konnte festgestellt werden, dass die Temperaturen von 
Zahnrad und Öl kurz vor Eingriff des Zahnrads entscheidend für eine ausreichende Schmierung 
sind. Im Kontext der Elastohydrodynamik wird dieser Zustand als Einlasszustand bezeichnet. 
Die Kenntnis der zu erwartenden Temperaturen ist demnach zwingend notwendig, um den 
tribologischen Zustand im Eingriff bestimmen zu Können. Hierfür muss der Wärmehaushalt 
mit allen Wärmequellen und -senken betrachtet werden. Durch die Kopplung von Reibwär- 
meentwicklung und gesamtem Temperaturniveau ist der Zusammenhang komplex. Auch in 
vereinfachten Betrachtungen muss zur Lösung des Wärmeleitungsproblems der Wärmeübergang 
durch Kühlung berücksichtigt werden. In der Literatur sind einige Ansätze zur Modellierung des 
Thermalhaushalts von hochbelasteten Zahnrädern bekannt und sollen im Folgenden geschildert 
werden. Eine Übersicht ist in Tabelle 2.2 dargestellt. 


Tabelle 2.2: Literaturstellen zur Thermalmodellierung von Zahnrädern 


Autor Art der Untersuchung Wärmeübergang abgebildet durch 
Patir und Cheng (1979) Numerisch DeWinter und Blok (1974) 
Handschuh (1992) Experimentell - 
Handschuh und Kicher (1994) Numerisch verschiedene! 
Long et al. (2003) Exp. und Numerisch De Winter und Blok (1974) 
Jie et al. (2013) Numerisch De Winter und Blok (1974) 
Xing und Shaojun (2013) Numerisch DeWinter und Blok (1974) 
Wang et al. (2016) Numerisch CFD 
Li und Tian (2017) Exp. und Numerisch verschiedene 


Patir und Cheng (1979) untersuchen in einer der frühesten Arbeiten zu dieser Fragestellung 
den Thermalhaushalt von Zahnrädern. Die Autoren modellieren das stationäre Temperaturfeld 
eines Zahnrads unter Einfluss von Reibwärme, Kühlung durch Einspritzkühlung und unter 
Berücksichtigung des stirnseitigen Wärmeübergangs. In der Arbeit wird ein numerisches Modell 
zur Bestimmung des Temperaturfeldes beschrieben, dessen Randbedingungen auf Korrelationen 
aus der Literatur beruhen. Das Modell ist angelehnt an eine Arbeit von Wang (1976) und erweitert 
das darin beschriebene Modell um Wärmeübergangskoeffizienten an der Zahnradoberfläche. Das 
Zahnrad wird als periodische Anreihung von Einzelzähnen betrachtet, um den Rechenaufwand 
gering zu halten. Die Randbedingungen des Finite-Elemente-Modells eines Einzelzahns werden 
in Zonen der Einspritzkühlung, der ungekühlten Zahnflanken und der Seitenflächen unterteilt. 
Das schematische Modell des Einzelzahns und die Aufteilung der Randbedingungen, die die 
Autoren beschreiben, werden in zahlreichen späteren Arbeiten übernommen. Das Modell des 
Einzelzahns ist in Abbildung 2.6 dargestellt. 


Der Wärmeeintrag in die belastete Zahnflanke wird in der Arbeit von Patir und Cheng (1979) 
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variiert (Bereich m). Der Wärmeübergangskoeffizient durch die Einspritzkühlung wird für den 
gesamten Zahnflankenbereich als konstant angenommen und durch die analytische Korrelation 
von DeWinter und Blok (1974) bestimmt (Bereiche t und m). Der Wärmeübergang an den 
Stirnflächen des Zahnrads wird durch eine Korrelation für rotierende Scheiben angenähert 
(Bereich s). Die periodischen Ränder (Bereiche g und p) werden als adiabate Ränder betrachtet. 
Die Analyse von Patir und Cheng (1979) erfolgt in dimensionsloser Form. 


Die Autoren stellen fest, dass die qualitative Temperaturverteilung stark vom Verhältnis der 
Wärmeübergangskoeffizienten an den Seitenflächen und den Zahnflanken abhängt. Eine überwie- 
gende Kühlung an den Zahnflanken verursacht demnach starke lokale Gradienten im Zahnbereich, 
während eine hauptsächlich stirnseitige Kühlung Gradienten über den gesamten Zahn hervorruft, 
wie in Abbildung 2.7 dargestellt ist. 


Die Arbeit von Patir und Cheng (1979) ist dahingehend eingeschränkt, dass der Wärmeübergangs- 
koeffizient als über die gesamte Zahnflanke konstant angenommen wird. Diese Einschränkung 
begründen die Autoren mit dem Mangel an Daten oder Korrelationen in der Literatur, die den 
Wärmeübergang ortsaufgelöst beschreiben. Gerade bei hochbelasteten Zahnrädern kann dies 
dazu führen, dass die auftretenden Temperaturen unterschätzt werden. Die Autoren betonen, 
dass keine gesicherten Daten für den Wärmeübergang durch Einspritzkühlung verfügbar sind 
und die Korrelation von DeWinter und Blok (1974) approximativer Natur ist. Die Erkenntnisse 
der Arbeit beschränken sich daher auf qualitative Beobachtungen. Die Autoren schlussfolgern, 
dass die quantitative Bestimmung der Temperaturniveaus und -verteilungen nur bei vorhandenen 
Daten zum Wärmeübergang erfolgen kann. 


Long et al. (2003) modellieren die Temperaturverteilung an einem belasteten, geradverzahnten 
Zahnrad durch das Aufstellen eines FE-Modells und greifen dabei das Modell von Patir und 
Cheng (1979) auf. Es wird angenommen, dass die Wärmezufuhr durch Reibung und die Kühlung 
durch die beteiligten Kühlmechanismen simultan ablaufen. Dies ermöglicht es, die Wärmelei- 


Abbildung 2.6: Periodisches Zahnelement und Wärmeübergangsgebiete nach Patir und Cheng 
(1979) 
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tungsgleichung stationär zu lösen. Während der Seitenwandwärmeübergang auf den gleichen 
Grundlagen wie in der Arbeit von Patir und Cheng (1979) beruht (Hartnett und Deland (1961), 
Popiel und Boguslawski (1975) und Dorfman (1963)), wird eine Korrektur eingeführt, um den 
Einfluss des Zweiphasengemisches aus Luft und Öl in Zahnradumgebung zu berücksichtigen. 
Dafür wird ein Stoffgemisch definiert, dessen Eigenschaften mit einem Quotienten d zwischen 
denen von Luft und denen des verwendeten Öls interpoliert werden. Um diesen Quotienten zu 
bestimmen, verwenden die Autoren Ergebnisse aus einer experimentellen Versuchskampagne. 
Analog zu Patir und Cheng (1979) wird der Wärmeübergang an den Zahnflanken durch die 
Korrelation von DeWinter und Blok (1974) abgeschätzt, die zu einem Wärmeübergangskoeffizi- 
enten umgeformt wird. Die Reibungswärme, die im Eingriff freigesetzt wird, wird anhand der 
Getriebekinematik und Reibungsbeiwerten aus experimentellen Arbeiten von Castro und Seabra 
(1998) und Winter und Michaelis (1984) (Gleichung 2.3) bestimmt und entlang der Zahnkontur 
aufgelöst ausgewertet. Wie in Kapitel 2.1 beschrieben sind die Schergeschwindigkeit vg und 
somit auch der Reibwärmeeintrag ortsabhängig, was erstmalig in der Arbeit von Long et al. 
(2003) berücksichtigt wird. 


Die Autoren vereinfachen das Modell zusätzlich, indem sie ein periodisches Element definieren. 
Dieses besteht aus einem Einzelzahn eines Zahnrads, das an seinen Schnittflächen adiabate Rän- 
der aufweist. Es werden sieben Oberflächengebiete definiert, an denen jeweils unterschiedliche 
thermische Randbedingungen angenommen werden. Der Zahn wird aufgeteilt in die belastete 
und unbelastete Flanke, die Zahnoberseite und den Zahnfuß. Neben den adiabaten Schnittflä- 
chen des periodischen Elements werden die Berührfläche mit der Welle und die Stirnseiten des 
Zahnrads als Ränder definiert. 


Long et al. (2003) variieren die Eingangsparameter des Modells und identifizieren so den 
Einfluss von Zahnlast und Drehzahl auf die Temperaturverteilung im Zahnrad. Dabei zeigt die 
Zahnradlast den größten Einfluss. Sowohl die mittleren als auch die maximalen Temperaturen 


(a) (d) 


Abbildung 2.7: Isolinien der Temperatur bei dominierendem Seitenwandwärmeübergang (links) 
und dominierendem Zahnflankenwärmeübergang (rechts) nach Patir und Cheng (1979) 
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auf dem Zahnrad steigen mit höheren übertragenen Lasten. 


Dadurch, dass der Quotient d iterativ durch die Annäherung der Simulationsergebnisse an die 
experimentellen Messwerte gewonnen wird, erreichen die Autoren eine Übereinstimmung der 
simulierten und gemessenen Temperaturen. Dabei liegt der resultierende Wärmeübergangskoeffi- 
zient für die Seitenflächen des Zahnrads in einer Größenordnung, die sonst nur in Staupunktnähe 
bei flüssiger Prallkühlung erreicht wird. Der Wärmeübergangskoeffizient der gekühlten Zahn- 
flanken ist um eine Größenordnung kleiner als der der Seitenflächen. 


Jie et al. (2013) beschreiben ein den obenstehenden Ansätzen ähnliches Thermalmodell eines 
Zahnrads. Die Bestimmung der Randbedingungen erfolgt analog zu den Arbeiten von Patir 
und Cheng (1979) und Long et al. (2003). Auch Xing und Shaojun (2013) bestimmen die 
Temperaturverteilung eines belasteten Zahnrads durch die Modellierung eines Einzelzahns. 
Dabei verwenden sie Korrelationen für rotierende Scheiben für die Seitenflächen, etablierten 
Zahnreibungskorrelationen für den Reibwärmeeintrag sowie der Wärmeübergangskorrelation 
für die einspritzgekühlten Flächen von DeWinter und Blok (1974). In beiden Arbeiten wird ein 
Temperaturmaximum am Ort des Wärmeeintrags beobachtet. Ebenfalls wird in beiden Arbeiten 
festgestellt, dass die Bestimmung genauer Randbedingungen unerlässlich ist, um verwertbare 
Abschätzungen der zu erwartenden Temperaturverteilungen treffen zu können. 


Wang et al. (2016) modellieren die instationäre Temperaturverteilung an einem belasteten, ein- 
spritzgekühlten Zahnrad. Dazu werden zunächst Wärmeübergangskoeffizienten bei vereinfach- 
tem Aufprallen eines Strahls auf die gekühlte Flanke durch numerische Strömungssimulationen 
berechnet. Der Wärmeübergang an der strömungsabgewandten Flanke wird ebenfalls numerisch 
bestimmt. Die so gewonnenen Werte werden als Randbedingungen eines zeitlich aufgelösten 
FE-Modells verwendet, um die instationäre Temperaturverteilung zu berechnen. Dabei werden 
Temperaturschwankungen von bis zu 17 K während eines Eingriffszyklus festgestellt. 


In einigen Arbeiten wird der Thermalhaushalt von Getrieben und Zahnrädern stark vereinfacht 
und als thermisches Netzwerk betrachtet (Manin und Play (1999), Koffel et al. (2009), Durand 
de Gevigney et al. (2012)). Häufig wird in diesen Arbeiten eine Getriebekonfiguration expe- 
rimentell untersucht, indem charakteristische Temperaturen an den Zahnrädern gemessen und 
mit den Berechnungen durch die thermischen Netzwerke verglichen werden. Auch hier stellt 
die Verfügbarkeit von Daten zum Wärmeübergang an Zahnrädern durch Einspritzkühlung eine 
Lücke dar. 


Im Vorangegangenen wurde eine Übersicht über die Literatur zur Thermalmodellierung gegeben. 
Verschiedene Ansätze und Erkenntnisse zur Bestimmung der Temperaturverteilung an Zahnrä- 
dern wurden vorgestellt. Die Arbeit von Patir und Cheng (1979) wird dabei häufig aufgegriffen, 
indem die Aufteilung der Gebiete des Wärmeübergangs und die verwendeten Korrelationen über- 
nommen werden. Die numerischen Arbeiten, in denen einspritzgekühlte Zahnräder betrachtet 
werden, verwenden dabei fast ausschließlich die Korrelation von DeWinter und Blok (1974), um 
den Wärmeübergang an der belasteten Zahnflanke abzubilden. In mehreren Arbeiten wird auf 
einen Mangel an gesicherten Daten zum Öl-Zahnrad-Wärmeübergang hingewiesen.” 


unter anderem durch Patir und Cheng (1979), Townsend und Akin (1981), Townsend und Akin (1983) und 
Xing und Shaojun (2013) 
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2.3 Prallkühlung an Zahnrädern 


Zur Implementierung der in Abschnitt 2.2 beschriebenen Modelle sind Daten zum Wärmeüber- 
gang durch Kühlung erforderlich. Die verfügbaren Daten in der Literatur sind äußerst begrenzt. 
Die in diesem Zusammenhang bekannten Arbeiten werden im Folgenden zusammengefasst. 


In den Thermalmodellen, die im vorangegangenen Abschnitt beschrieben werden, wird häufig 
auf die Korrelation von DeWinter und Blok (1974) zurückgegriffen. Aufgrund der prominenten 
Bedeutung dieser Korrelation soll der ihr zu Grunde liegende Ansatz beleuchtet werden. Die 
Autoren betrachten den Wärmeübergang an einem Zahnrad analytisch, indem sie ein stark 
vereinfachtes Modell des fluidmechanischen Vorgangs aufstellen. Dazu wird folgendermaßen 
vorgegangen: 


Die Zahnlücke wird als rechteckige Geometrie betrachtet (siehe Abbildung 2.8). 


Zu Beginn des Vorgangs füllt das Fluid den gesamten Zahnquadranten und ist in Ruhe. 


Die Temperaturdifferenz zwischen Oberfläche und Fluid wird als gering betrachtet. Es 
wird angenommen, dass daher auch die Viskositätsunterschiede vernachlässigt und die 
Impuls- und Energiegleichung entkoppelt werden können. 


Die zeitliche Entwicklung der Filmdicke und Geschwindigkeitsverteilung des Fluids 
werden durch eine analytische Lösung der Impulsgleichung bestimmt. Durch die Annahme 
konstanter Viskosität ist diese unabhängig vom Wärmeeintrag in das Fluid. 


Um den Wärmeübergang von der Oberfläche in das Fluid zu bestimmen, wird die Energie- 
gleichung numerisch gelöst. 


Durch Entdimensionieren und numerisches Lösen der Impulsgleichung und der Energieglei- 
chung bestimmen DeWinter und Blok (1974) die Wärme, die zu einem beliebigen Zeitpunkt 
durch das Fluid aufgenommen wurde. Daraus lassen sich unmittelbar der zeitliche Verlauf des 
Wärmeübergangs und seine obere Grenze bestimmen. Die über den gesamten Vorgang des 
Abschleuderns des Öls maximal abführbare Wärme nach DeWinter und Blok (1974) beträgt 


AT 
Otot,max = 5,6m ADE Tes (2.16) 


oder als Wärmeübergangskoeffizient: 


@ Otot,max vo 
ST AT DIR 2,5454 /Apc (2.17) 
Aus der Betrachtung von Gleichung 2.17 wird bereits die approximative Natur der Korrelation 
von DeWinter und Blok (1974) deutlich. Durch die getroffene Annahme, dass zu Beginn des 
Kühlvorgangs die gesamte Zahnlücke stets mit Öl gefüllt ist, findet der eingespritzte Ölvolumen- 
strom keine Berücksichtigung bei der Bestimmung des Wärmeübergangs. Dadurch, dass das 
Fluid zu diesem Zeitpunkt als auf der Flanke ruhend angenommen wird, Kann auch der Impuls 
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Abbildung 2.8: Auszentrifugieren eines Fluidfilms in einem Zahnquadranten nach DeWinter und 
Blok (1974) 


des auf die Zahnradoberfläche auftreffenden Ölstrahls in Gleichung 2.17 nicht berücksichtigt 
werden. 


Townsend und Akin (1983) untersuchen das Einsetzen von Warmfressen an einem belasteten 
Zahnradpaar experimentell und analytisch. Die Autoren beschreiben einen Prüfstand, in dem 
ein Zahnradpaar im Eingriff steht und durch eine Verspannmechanik unter einem einstellbaren 
Drehmoment bei unterschiedlichen Drehzahlen betrieben werden kann. Zusätzlich werden so- 
wohl der Einspritzdruck des Öls als auch der Düsendurchmesser variiert. In den durchgeführten 
Experimenten wird die Temperatur der belasteten Zahnflanke unter Verwendung eines Pyrome- 
ters gemessen. Dadurch, dass das verwendete Pyrometer zeitlich hochaufgelöste Messungen 
zulässt und dadurch, dass der Schnittpunkt der Sichtachse des Pyrometers und der Zahnflanke 
während des Durchlaufs eines Einzelzahns entlang der Flanke wandert, können ortsaufgelöste 
Temperaturen erfasst werden. Zusätzlich wird ein Thermalmodell vorgestellt, das das Modell 
von Patir und Cheng (1979) um die Eindringtiefe des Ölstrahls erweitert. Die Eindringtiefe wird 
dabei durch das vektorielle Modell von Akin et al. (1975) bestimmt. Der Wärmeübergangsko- 
effizient für Einspritzschmierung wird von der Zahnspitze bis zur Eindringtiefe angenommen. 
Den Wärmeübergangskoeffizienten, der durch die Korrelation von DeWinter und Blok (1974) 
bestimmt wird, bewerten die Autoren als zu gering und verwenden daher einen nicht näher 
beschriebenen eigenen Schätzwert. 


Die experimentellen Ergebnisse von Townsend und Akin (1983) zeigen eine Abnahme der 
Zahnflankentemperaturen mit steigendem Düsenvordruck für alle untersuchten Betriebspunkte. 
Der Effekt ist dann besonders ausgeprägt, wenn das übertragene Moment oder die Drehzahl hoch 
sind. Die geringsten Temperaturen am Zahnrad werden dann erreicht, wenn der eingespritzte 
Ölstrahl den Zahnfuß erreicht. Der beobachtete Zusammenhang zwischen Düsenvordruck und 
Zahnradtemperatur ist in Abbildung 2.9 für einige gemessene Betriebspunkte von Townsend 
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und Akin (1983) dargestellt. Dabei wird ein Moment von 3736 N/cm übertragen. Durchgezo- 
gene Linien stellen mittlere Zahnflankentemperaturen, gestrichelte Linien die Minimal- und 
Maximaltemperaturen dar. 
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Abbildung 2.9: Einfluss des Düsenvordrucks und des übertragenen Moments auf die Zahnrad- 
temperatur nach Townsend und Akin (1983) 


Der Vergleich der experimentellen und numerischen Ergebnisse liefert allgemein tibereinstim- 
mende Ergebnisse. Die getroffene Abschätzung des Wärmeübergangs wird jedoch nicht erläutert. 
Dadurch, dass im Experiment die unbelastete Zahnflanke gekühlt wird, während die Messungen 
an der belasteten Zahnflanke aufgezeichnet werden, kann die örtliche Verteilung des Wärme- 
übergangskoeffizienten nicht bestimmt werden. 


Schober (1983a) untersucht den Einfluss der Düsenanordnung und der verwendeten Ölmenge 
auf das Einsetzen von Fressen an einem Zahnradpaar in einem Verspannungsprüfstand. Düsenan- 
ordnungen, bei denen Öl in den beginnenden Eingriff, in den auseinandergehenden Eingriff und 
in radiale Richtung auf den Zahnradumfang gespritzt wird, werden untersucht. Die notwendigen 
Ölmengen, um Fressen zu vermeiden, sind dabei wesentlich geringer, als die im Vorfeld durch 
gängige Auslegungsmethoden bestimmten Ölmengen. 


Für eine Versuchsreihe mit Einspritzung in den beginnenden Eingriff werden Zahnradtempe- 
raturen unter Last und unter Variation des Ölvolumenstroms gemessen. Dabei führen höhere 
Ölvolumenströme zu sinkenden Zahnradtemperaturen. Durch Berechnung und Bilanzierung der 
Verlustleistungen und der abgeführten Wärmeströme bestimmt Schober (1983a) die anteilige 
Kühlwirkung des Öls. Diese steigt bei höheren Ölvolumenströmen an, während bei geringen 
Ölvolumenströmen die Kühlung durch Konvektion am Prüfstandsgehäuse und Wärmeleitung 
durch das Prüftstandsfundament dominiert. 


Schober führt die beobachteten geringen Ölmengen zur Vermeidung von Fressen auf die geringen 
zum vollständigen Aufbau eines tragfähigen EHD-Schmierfilms erforderlichen Ölmengen zu- 
rück. Das Einsetzen von Fressen bei geringen Ölvolumenströmen wird von Schober auf kritische 
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Zahnradtemperaturen durch geringere Kühlung und nicht auf mangelhafte Ölbenetzung zurück- 
geführt. Es wird hervorgehoben, dass im Fall sehr geringer Schmiermengen auch geringere 
Verluste zu erwarten sind. 


Handschuh (1992) untersucht die transiente Temperaturverteilung an schrägverzahnten Kegelrä- 
dern. Dabei wird ein Fächerstrahl eingesetzt, der die gesamte Zahnbreite abdeckt. Die Zahnräder 
werden durch eine Verspannmechanik, ähnlich der von Townsend und Akin (1983), vorgespannt 
und bei unterschiedlichen Drehmomenten und Drehzahlen getestet. Der Einfluss der Position 
des Ölstrahls wird für insgesamt fünf Positionen um den Umfang des Zahnrads untersucht. Die 
Oberflächentemperaturen werden ebenfalls durch Pyrometrie gemessen. Die Temperatur wird 
zusätzlich durch in das Zahnrad eingelassene Thermoelemente erfasst. Auch hier führt erhöhte 
Ölzufuhr zu geringeren Temperaturen für alle untersuchten Betriebspunkte und Düsenpositionen. 
Der Temperaturunterschied bei den untersuchten Ölvolumenströmen ist jedoch wesentlich gerin- 
ger als der, der von Townsend und Akin (1983) beobachtet wird. Die übertragene Leistung hat 
den stärksten Einfluss auf die resultierenden Temperaturen. Die Temperaturen sind stets an den 
Zahnspitzen am höchsten und nehmen zum Zahnfuß hin ab. Die Düsenpositionen unmittelbar 
vor und nach dem Zahneingriff führen zu den niedrigsten Temperaturen und weisen somit die 
höchsten Wärmeübergänge auf. Es wird berichtet, dass die Einspritzung unmittelbar vor dem 
Eingriff jedoch mit erhöhten Verlusten durch die Verdrängung des Öls verbunden ist. 


Mizutani et al. (1989) bestimmen die Temperaturverteilung eines schnelllaufenden Zahnradpaares 
und berechnen die anteiligen Beiträge von Reibungs-, Plansch- und Ventilationsverlusten. Sie 
stellen fest, dass bei hohen Drehzahlen die Plansch- und Ventilationsverluste dominieren. Eine 
Bestimmung des Wärmeübergangskoeffizienten erfolgt nicht. 


Leoni (1991) untersucht die Eignung getrennter Schmierung und Kühlung zur Reduktion der 
eingespritzten Ölmenge. Durch die Messung der Temperaturen an zwei Stellen eines Zahnrads 
kann der Wärmestrom durch das Zahnrad näherungsweise bestimmt werden. Aus Bilanzbe- 
trachtungen wird der Wärmeübergangskoeffizient berechnet. Höhere Ölvolumenströme sowie 
höhere Umfangsgeschwindigkeiten führen in der Untersuchung zu höheren Wärmeübergangs- 
koeffizienten. Die höchsten untersuchten Ölvolumenströme bilden hier eine Ausnahme und 
führen zu geringeren ermittelten Wärmeübergangskoeffizienten. Aufgrund des Versuchsaufbaus 
können die Verlustwärme, die im Eingriff entsteht und der Wärmeübergang durch Kühlung 
nicht voneinander unterschieden werden. Leoni (1991) führt den unerwarteten Rückgang des 
Wärmeübergangs bei hohen Ölvolumenströmen auf diese Einschränkung zurück. 


2.4 Wärmeübergang durch Prallkühlung 


Zahlreiche technische Anwendungen benötigen Kühlmethoden mit hohen Wärmeübergangskoef- 
fizienten. Der Einsatz von Prallkühlung zur Erzeugung sehr hoher Wärmeübergänge ist effektiv, 
robust und konstruktiv meist einfach, was zu einer weiten Verbreitung dieser Methode geführt 
hat. Beim Walzen von Stählen werden sowohl die Walzen als auch das Material durch Prall- und 
Sprühkühlung gekühlt, um thermische Überlastung zu verhindern und die Prozesstemperaturen 
zu optimieren (Chen und Tseng (1992)). In der jüngeren Vergangenheit fand durch den steigen- 
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den Kühlungsbedarf von Computerchips auch hier Prallkühlung Anwendung (Ebadian und Lin 
(2011)). Durch zunehmende Miniaturisierung steigen die Anforderungen an die erreichbaren 
Wärmestromdichten zusätzlich an. In bestimmten militärischen Mikroelektronikanwendungen 
treten Wärmestromdichten in der Größenordnung von 1 MW /m? auf. In zukünftigen Anwendun- 
gen werden 2,5 MW /m? erwartet, die durch effiziente Kühlsysteme abgeführt werden müssen. 
Auch diesen Anforderungen kann durch Prallkühlung begegnet werden. 


Im Bereich der Turbomaschinen wird ebenfalls auf Prallkühlung zurückgegriffen, um thermisch 
hoch belastete Oberflächen und Strukturen zu kühlen. So wird zur Kühlung der Vorderkanten von 
Laufschaufeln Prallkühlung eingesetzt, da nur im Bereich der Vorderkante durch die typische, 
geringe Dickenrücklage genügend Bauraum verfügbar ist, ohne die strukturelle Integrität zu 
gefährden. Bei Leitschaufeln, die keiner Fliehkraftlast standhalten müssen, wird die Prallkühlung 
auch in der Mitte der Sehnenlänge eingesetzt (Han et al. (2013)). 


Wie im vorangegangenen Abschnitt verdeutlicht wird, ist Prallkühlung mittels eines Ölstrahls in 
Getrieben ab gewissen Leistungsdichten notwendig, um ausreichend hohe Wärmestromdichten 
abführen zu können. Im Folgenden sollen die Grundlagen der Prallkühlung geschildert werden, 
um das bereits verfügbare Wissen und die Analogien zum Anwendungsfall aufzuzeigen. 


In der Literatur wird das Feld der Prallkühlung häufig unterteilt in Prallkühlung durch 


e Gasstrahlen (gas jet impingement) 


e Flüssigkeitsstrahlen 
in gasförmiger Umgebung (free surface jet, unsubmerged jet) 


in flüssiger Umgebung (submerged jet) 


e Flüssigkeitsstrahlen mit Phasenübergang. 


Stellenweise wird von einphasiger (single phase) bzw. zweiphasiger (two phase) Prallkühlung 
gesprochen, womit Strömungen ohne bzw. mit Phasenübergang gemeint sind. Die Literatur zur 
Beschreibung und Vorhersage der Strömung und des Wärmeübergangs durch Gasstrahlen ist 
umfangreich und wird von Downs und James (1987) sowie Han und Goldstein (2001) zusam- 
mengefasst. Im vorliegenden Fall ist ausschließlich die Kühlung durch Flüssigkeitsstrahlen in 
gasförmiger Umgebung relevant. Einige wichtige Aspekte dieser Form der Prallkühlung sollen 
im Folgenden umrissen werden. Für eine ausführliche Darstellung sei auf die Übersichtsarbeit 
von Webb und Ma (1995) verwiesen. 


In Arbeiten zur Prallkühlung werden die Strömung und der Wärmeübergang meist durch dimen- 
sionslose Kennzahlen beschrieben. Bei Rundstrahlen ist die Definition der Reynoldszahl mit 
der mittleren Strahlgeschwindigkeit us, dem Austrittsdurchmesser des Rundstrahls d und der 
kinematischen Viskosität v üblich, als: 


Ree (2.18) 
v 
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Die Viskosität, die Wärmeleitfähigkeit, die Dichte und die spezifische Wärmekapazität des Fluids 
werden über die Prandtlzahl ausgedrückt, mit: 


v 
—, (2.19) 
a 

und der sogenannten Temperaturleitfähigkeit a = 2, Der Wärmeübergangskoeffizient h kann 
durch die Nußeltzahl dimensionslos ausgedrückt werden: 


hd 


Nu 7 


(2.20) 
Darin wird im Zusammenhang von Rundstrahlen d als charakteristische Länge eingesetzt, À ist 
die Wärmeleitfähigkeit des Fluids. 


2.4.1 Prallkühlung durch Flüssigkeitsstrahlen 


Die hohen Wärmeübergangskoeffizienten der Prallkühlung werden erreicht, indem Konvektion 
erzwungen wird. Durch den Strahlaufprall an der zu kühlenden Oberfläche wird die Strömung 
umgelenkt, es werden hohe Geschwindigkeitsgradienten in Wandnähe hervorgerufen. Die 
Eigenschaften der sich durch Wandreibung ausbildenden Strömungsgrenzschicht und der damit 
gekoppelten thermischen Grenzschicht sind entscheidend für den erreichbaren Wärmeübergang. 


Free Surface — 


Boundary 
Layer stag 


Approximate Stagnation 
Zone 
Abbildung 2.10: Ausbildung einer Grenzschicht in der Stauströmung durch Prallkühlung mit 
einem unbewegten, flüssigen Rundstrahl in Gasumgebung nach Lienhard (2006) 


Durch Reibung an der Wand und den dadurch entstehenden Impulsverlust in der Strömung 
nimmt die Dicke der Strömungsgrenzschicht eines Prallstrahls in Strömungsrichtung zu. Im 
Gegensatz zur Umströmung einer ebenen Platte oder von Tragprofilen kann die Grenzschicht 
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jedoch nicht ungehindert anwachsen. Bei achsensymmetrischen Prallstrahlen nimmt die Dicke 
des Flüssigkeitsfilms aus Kontinuitätsgründen stetig ab (siehe Abbildung 2.10). An einem 
bestimmten Radius erreicht die Dicke der Stromungsgrenzschicht die Filmhöhe. Ab diesem 
Punkt ist die gesamte Strömung als reibungsbehaftet zu betrachten. Zusätzlich kann die Strömung 
an einem weiter vom Staupunkt entfernten Radius einen hydraulischen Sprung erfahren, an 
dem sie vom schießenden Zustand durch eine schlagartige Dickenzunahme in den strömenden 
Zustand übergeht. Durch die genannten Phänomene sind die Beschreibung und Vorhersage 
der Strömung und des Wärmeübergangs komplex. In der Literatur wird der Wärmeübergang 
durch Prallkühlung daher in Bereiche unterteilt, die sich durch den Charakter ihrer Strömung 
unterscheiden. Die Unterteilung ist in Abbildung 2.11 dargestellt. 


Pin To Ti 
o i poe oe k Fully Turbulent 
Developing Turbulence 


Stagnation Zone 


Abbildung 2.11: Unterteilung eines auftreffenden Prallstrahls in charakteristische Strömungsge- 
biete nach Liu et al. (1991) 


Verschiedene Arbeiten befassen sich mit der Strömung und dem Wärmeübergang in einzelnen der 
in Abbildung 2.11 dargestellten Bereiche. Liu et al. (1991) bestimmen den Wärmeübergang durch 
Prallkühlung in der Stauzone von Flüssigkeitsstrahlen gleichförmigen Geschwindigkeitsprofils 
anhand analytischer Überlegungen zu 


0,715Rey Prt, 0,15<Pr <3 

Nu = peel (2.21) 
0,797 Reg” Pr3, 3<Pr 

Dies steht in guter Übereinstimmung mit den Ergebnissen von Ma et al. (1996), die den Wär- 


meübergang an Oberflächen unter Berücksichtigung inhomogener Randbedingungen analytisch 
bestimmen. Der Einfluss der Randbedingungen wird dabei über einen Faktor € berücksichtigt. 
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Unter homogenen Wärmestromrandbedingungen gilt € = 1 und somit nach Ma et al. (1996) 


0,7212Rey? Pr'*, 0,7<Pr <3 
Nu = ¢ 0,7212Re% Pr? 3<Pr <10 . (2.22) 
0,8597Re) > Pr3,  10<Pr 


Durch Gleichungen 2.21 und 2.22 ist der Wärmeübergang im Staupunktbereich durch einfache 
Ausdrücke und für die Prandtlzahlen der meisten Fluide, die in technischen Anwendungen 
eingesetzt werden, bestimmbar. Der Wärmeübergang in den stromab anknüpfenden Bereichen 
hängt in großem Maße vom Wärmeübergang im Staupunktbereich ab. Durch die außerhalb 
des Staupunktbereichs anwachsenden Grenzschichten nimmt der Wärmeübergang stetig ab. 
Lediglich durch den Umschlag zur turbulenten Strömung kann ein positiver Gradient des lokalen 
Wärmeübergangs in Strömungsrichtung beobachtet werden. Analog zu Gleichungen 2.21 und 
2.22 sind für die Bereiche stromab des Staupunkts analytische Ausdrücke zur Bestimmung 
des Wärmeübergangs bekannt. Da für laminare und turbulente Strömungen und verschiedene 
Bereiche der Prandtlzahl unterschiedliche, teils sehr Komplexe Ausdrücke zur Anwendung 
kommen, ist die Darstellung der Theorie umfangreich und soll im Rahmen dieser Arbeit nicht 
erfolgen. Ein ausführlicher Überblick wird von Webb und Ma (1995) gegeben. Nahezu allen 
Korrelationen für die unterschiedlichen Bereiche der Strömung ist jedoch der Einfluss der 
Reynoldszahl auf den Wärmeübergang gemeinsam, je nach Fall entweder zu 


Nu œ Re! (2.23) 


oder zu 
Nu œ Re®333, (2.24) 


Aufgrund des hohen Dichteverhältnisses zwischen Flüssigkeit und Luft in den meisten Anwen- 
dungsfällen ist die Beeinflussung des Strahls durch die ihn umgebende Luft gering. Dadurch 
ist der Einfluss des Abstands zwischen Düse und Oberfläche gering. Stevens und Webb (1991) 
berichten eine Abhängigkeit des Wärmeübergangs in Abhängigkeit des Düsenabstands von der 
Oberfläche, z, mit 

z ) —0,052 


Nu « © (2.25) 


2.4.2 Prallkühlung an bewegten Oberflächen 


Der Wärmeübergang durch Prallkühlung an bewegten Oberflächen wird in einigen Arbeiten 
diskutiert. Carper et al. (1986) beschreiben den gemittelten Wärmeübergang an einer rotierenden 
Scheibe bei konzentrischen Rotations- und Strahlachsen durch den Ausdruck 


Nup = 0,097 Reg Reo? pr0:448, (2.26) 


mit der Rotationsreynoldszahl 


Reg = —. (2.27) 
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Aus Gleichung 2.26 ist ersichtlich, dass eine Steigerung der Oberflächengeschwindigkeit ortho- 
gonal zur Strahlachse den flächen- und zeitgemittelten Wärmeübergang erhöht. Dies liegt vor 
allem an den dünneren Grenzschichten, die sich durch die laterale Bewegung der Oberfläche 
relativ zum Strahl ausbilden. Eine Erhöhung der Strahlreynoldszahl Reg erhöht ebenfalls den 
mittleren Wärmeübergang an der Scheibe mit Nu « Rey”. Somit liegt trotz grundsätzlich ande- 
rer Zustände ein vergleichbarer Einfluss der Strahlgeschwindigkeit im Vergleich zur Prallkühlung 
unbewegter Oberflächen vor. Lavalle et al. (2007) bestätigen die Korrelation von Carper et al. 
(1986) durch experimentelle Daten mit Abweichungen im Bereich von etwa 20%. Gleiches gilt 
für das Auftreffen eines Strahls außerhalb des Mittelpunkts der Scheibe mit 


Nup = 0,57Reg Rey“ Pr?8! (1 —27/D)°>. (2.28) 


2.4.3 Instationäre Prallkühlung 


In den vorangegangenen Betrachtungen konnte stets von stationären Strömungszuständen und 
thermischen Randbedingungen und somit zeitlich Konstanten Wärmeübergängen ausgegangen 
werden. Im Kontext dieser Arbeit sind jedoch periodisch instationäre Zustände zu erwarten. 
Zumbrunnen und Aziz (1993) untersuchen die instationäre Prallkühlung mit dem Ziel, den 
erreichbaren Wärmeübergang zu erhöhen. Die Autoren beschreiben einen Versuchsaufbau, in 
dem ein flüssiger, ebener Freistrahl durch rotierende Schneiden periodisch unterbrochen wird. 
Als eben gelten dabei Freistrahlen, wenn ihre Tiefenausdehnung wesentlich größer als ihre Brei- 
tenausdehnung ist. Durch den Vergleich der Wärmeübergangskoeffizienten im stationären und 
im intermittierenden Fall wird eine Erhöhung des mittleren Wärmeübergangs festgestellt. Dabei 
stellen die Autoren kein asymptotisches Verhalten des Wärmeübergangs mit steigender Frequenz 
der Intermittenz fest, wie es bei zunehmender Turbulenz des Strahls der Fall ist. Außerdem 
wird im gesamten betrachteten Bereich der Oberfläche eine Erhöhung des Wärmeübergangs 
beobachtet, während bei Wärmeübergangssteigerungen durch Turbulenzerhöhung des Strahls 
nur lokale Erhöhungen beobachtet werden. 


Bei vollständig intermittierenden Strahlen, bei denen die Strömung kurzzeitig unterbrochen wird, 
kann eine Erhöhung des mittleren Wärmeübergangs beobachtet werden. Maximale Steigerungen 
des Wärmeübergangs von 33% (Zumbrunnen und Aziz (1993)) bei Flüssigkeitsstrahlen sowie 
30% (Herwig und Middelberg (2008)) und 15% (Herwig et al. (2004)) bei Gasstrahlen werden 
beobachtet. Dabei wird der Wärmeübergang durch instationäre Strahlen stets im Verhältnis zur 
stationären Anströmung mit der maximalen Geschwindigkeit, die im instationären Fall auftritt, 
verglichen. Jedoch können zeitlich veränderliche Strahlen den mittleren Wärmeübergang auch 
reduzieren. Im Fall von Strahlen, deren Geschwindigkeitsverlauf einen Sinusverlauf annimmt, 
werden Rückgänge des Wärmeübergangs von 20% (Sheriff und Zumbrunnen (1994)) beobachtet. 
In diesem Fall überwiegt der ungewünschte Effekt der Unterbrechung der Fluidzufuhr und damit 
des Wärmeübergangs den günstigen Effekt der geringeren Grenzschichtdicke. 


Das Auftreten von transienten Zuständen der Oberflächentemperaturen der thermischen Grenz- 
schicht führt ebenfalls zu starken Änderungen des Wärmeübergangs im Vergleich zu stationären 
Zuständen (Zumbrunnen (1992)). Dieser Effekt wird vor allem bei kleinen Prandtlzahlen beob- 
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Abbildung 2.12: Verteilung des Wärmeübergangs durch intermittierende Prallkühlung bei unter- 
schiedlichen Unterbrechungsfrequenzen nach Zumbrunnen und Aziz (1993) 


achtet, da bei Pr < 1 auch die freie Außenströmung maßgeblich am Wärmetransport beteiligt 
ist. 


2.4.4 Zusammenfassung 


Aus dem Überblick über die verfügbare Literatur ist ersichtlich, dass der Wärmeübergang durch 
stationäre, achsensymmetrische, orthogonale Prallstrahlen gut beschrieben ist. Der isolierte 
Einfluss verschiedener Abweichungen von diesem grundlegenden Vergleichsfall wird ebenfalls 
in einigen Arbeiten behandelt. So sind der Wärmeübergang durch pulsierende Strahlen und die 
Kühlung bewegter Oberflächen in gewissen Grenzen mathematisch beschrieben. Auch die hier 
nicht dargelegte Anströmung unter einem nicht orthogonalen Winkel wird durch einige Arbeiten 
beschrieben. Die folgenden Einzeleffekte sind bekannt: 


* Der Wärmeübergang durch Prallkühlung verhält sich in weiten Bereichen proportional zu 
Rej mit 0,333 < j < 0,5. 


e Zeitliche Schwankungen des Strahls können den mittleren Wärmeübergang erhöhen, sofern 
es sich um eine intermittierende Strömung handelt. Sinusmodulierte Geschwindigkeits- 
schwankungen führen zu geringeren mittleren Wärmeübergängen. 


e Der Einfluss des Anstellwinkels des Strahls auf den Wärmeübergang ist gering. Prallstrah- 
len, die unter einem Winkel von unter 90° auftreffen, führen zu geringeren Wärmeüber- 
gangskoeffizienten. 


Bei der Prallkühlung von Zahnrädern sind alle genannten Effekte zu erwarten. Die Bestim- 
mung des zu erwartenden Wärmeübergangs durch die Öl-Zahnrad-Interaktion ist nach dem 
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derzeitigen wissenschaftlichen Kenntnisstand nicht möglich, da keine Informationen über die 
Wechselwirkungen zwischen den jeweiligen Effekten vorliegen. 


2.5 Zielsetzung dieser Arbeit 


In Abschnitt 2.1 wird dargestellt, weshalb die Schmierung und Kühlung von Getrieben unerläss- 
lich ist, gleichzeitig aber besonders für schnelldrehende Getriebe zu erheblichen Verlusten führen 
kann. Daher sind ein sparsamer Umgang und eine möglichst effiziente Zufuhr des Kühlmittels 
unerlässlich. Die Elastohydrodynamik, die in Abschnitt 2.1.2 beschrieben wird, zeigt jedoch, 
dass nur durch die Einhaltung auslegungsspezifischer Temperaturmaxima auf der Zahnradober- 
fläche Schäden an den Zahnrädern vermieden und somit ausreichende Lebensdauer sichergestellt 
werden kann. Dazu muss der Thermalhaushalt der Zahnräder verstanden und dementsprechend 
konzipiert werden. Abschnitte 2.2 und 2.3 zeigen auf, dass der gegenwärtige wissenschaftliche 
Kenntnisstand dazu nicht ausreicht. Die meisten bekannten Modellierungsansätze greifen auf 
die Korrelation von DeWinter und Blok (1974) zurück, um den durch die Einspritzkühlung her- 
vorgerufenen Wärmeübergang abzuschätzen. Die Ergebnisse experimenteller Untersuchungen 
belasteter Zahnradpaare können stets nur mäßig mit den aufgestellten Modellen in Einklang 
gebracht werden. Dies liegt vor allem darin begründet, dass die Korrelation von DeWinter 
und Blok (1974) analytischer Natur ist und auf äußerst weitreichenden Vereinfachungen der 
Strömungs- und Wärmetransportvorgänge beruht. 


Eine experimentelle Datenbasis zum Zahnradwärmeübergang existiert nicht. Die bekannten 
experimentellen Untersuchungen erlauben in keinem Fall die Trennung zwischen der Reibwärme, 
die im Zahnradeingriff erzeugt wird und dem Wärmeübergang durch Einspritzkühlung. Die Be- 
rechnung der Wärmeübergangskoeffizienten für die veröffentlichten experimentellen Daten durch 
das Berücksichtigen des Reibwärmeeintrags ist nicht möglich. Durch die Elastohydrodynamik 
kann der Reibwärmeeintrag nur näherungsweise bestimmt werden. Die weiteren thermischen 
Randbedingungen der Untersuchungen sind in keinem der Fälle genau genug bekannt. Örtlich 
aufgelöste Wärmeübergangskoeffizienten sind ebenfalls nicht aus den Daten bestimmbar. 


Aus der Zusammenfassung der Arbeiten zur Prallkühlung in Abschnitt 2.4 wird deutlich, dass 
der Wärmeübergang für einige generische Fälle eingehend beschrieben wird. Auch einzelne 
Einflüsse, wie Prallkühlung durch Flüssigkeitsstrahlen unter einem Winkel, sind innerhalb 
gewisser Grenzen beschrieben. Dadurch, dass im Fall der Einspritzkühlung in Getrieben jedoch 
mehrere Effekte simultan auftreten, reicht der Stand der Wissenschaft nicht zur Bestimmung des 
Öl-Zahnrad-Wärmeübergangs aus. 


Um die Entwicklung hochbelasteter Getriebe zu ermöglichen, muss dieser Mangel behoben und 
der Wärmeübergang zwischen Kühlmittel und Zahnrad experimentell ermittelt werden. Dies ist 
einerseits notwendig, um eine effiziente Auslegung derartiger Getriebe zu gewährleisten. Ande- 
rerseits sind experimentelle Daten unerlässlich, um CFD-Modelle der Ol-Zahnrad-Interaktion 
und des Öl-Zahnrad-Wärmeübergangs zu validieren. 


Die vorliegende Arbeit schließt diese Lücke. Dazu werden eine experimentelle Methode und eine 
Versuchsanordnung vorgestellt, mit welchen der Wärmeübergang an einem einspritzgekühlten 
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Zahnrad bestimmt werden kann. Die Bestimmung des Wärmeübergangs erfolgt örtlich aufgelöst. 
Thermographische, örtlich hochaufgelöste Messungen der Oberflächentemperatur, die häufig 
zur Bestimmung des Wärmeübergangs genutzt werden, sind im vorliegenden Fall nicht möglich. 
Daher wird in dieser Arbeit eine Methode vorgestellt, durch die der Wärmeübergang mit einer 
örtlich begrenzten Anzahl von Messpunkten bestimmt werden kann. Die Messunsicherheit 
der vorgestellten Methode wird statistisch untersucht. Eine Versuchsanordnung, mit der der 
Wärmeübergang am Zahnrad unter Laborbedingungen erzeugt werden kann, wird vorgestellt. 


Zur Erfüllung der Zielsetzung wird folgendermaßen vorgegangen: 


1. Identifikation der relevanten Parameter für den Öl-Zahnrad-Wärmeübergang 


2. Entwicklung einer experimentellen Methode zur Messung des Wärmeübergangs und 
Bestimmung ihrer Messunsicherheit 


3. Umsetzung der Methode in einem Versuchsaufbau 
4. Charakterisierung des Öl-Zahnrad-Wärmeübergangs unter Variation der Parameter 


5. Ableitung einer Korrelation zur Abbildung des mittleren Wärmeübergangs 


3 Untersuchte Öleinspritzkühlung 


Die Analyse des Stands der Wissenschaft in Kapitel 2 zeigt, dass der Kenntnisstand zum Wär- 
meübergang durch Einspritzkühlung lückenhaft ist. Die wenigen existierenden Arbeiten zum 
Öl-Zahnrad-Wärmeübergang befassen sich zudem mit unterschiedlichen konstruktiven Aus- 
gestaltungen der Öleinspritzung, weshalb ein Vergleich nur qualitativ möglich ist. So ist die 
Arbeit von Handschuh (1992) die einzige Arbeit zum Wärmeübergang durch aufgefächerte 
Ölstrahlen. Andere Arbeiten (Townsend und Akin (1981), Townsend und Akin (1983), Akin und 
Townsend (1985)) beschreiben den Wärmeübergang durch zylindrische Düsen bzw. Ölstrahlen. 
Die Korrelation von DeWinter und Blok (1974) wird mit der Modellvereinfachung abgeleitet, 
dass der Raum zwischen zwei Zahnflanken nach Einspritzen durch die Düse mit Öl aufgefüllt ist, 
wodurch die Form der Ölzufuhr außerhalb der Betrachtung liegt. Darüber hinaus wird in den 
bekannten Quellen die geometrische Ausgestaltung meist nur schemenhaft geschildert. 


Um die Zielsetzung der vorliegenden Arbeit zu erfüllen, wird der Wärmeübergang an einem 
einspritzgekühlten Zahnrad untersucht. Um die Ergebnisse in den Kontext der wissenschaftlichen 
Literatur einordnen zu Können und nachfolgenden Arbeiten die Anknüpfung zu ermöglichen, 
soll in diesem Kapitel die untersuchte Öleinspritzkühlung charakterisiert werden. Insbesondere 
für numerische Untersuchungen soll so eine Validierungsbasis geschaffen werden, die bisher 
nicht existiert. 


Zum derzeitigen Punkt sind keine Arbeiten bekannt, in denen der Öl-Zahnrad-Wärmeübergang 
isoliert von anderen Wärmeströmen bestimmt wird. Es ist keine Systematik bekannt, nach der die 
geometrischen und strömungsmechanischen Parameter der Kühlung klassifiziert werden Können. 
Zur Erfüllung der Zielsetzung dieser Arbeit müssen daher in einem ersten Schritt die für den 
Öl-Zahnrad-Wärmeübergang relevanten Parameter identifiziert werden. Dies soll im Folgenden 
anhand der untersuchten Öleinspritzung geschehen. Anschließend werden die physikalischen 
Parameter betrachtet, die einen Betriebspunkt definieren. Danach werden die zur Beschreibung 
des Ölaufpralls verwendeten abgeleiteten Größen eingeführt. 


3.1 Öleinspritzkühlung 


Geometrische Parameter 


Die betrachtete Kühlmethode nutzt die hohen erreichbaren Wärmeübergangskoeffizienten der 
Prallkühlung, indem Öl in Form eines Ölstrahls direkt auf die zu kühlende Zahnradoberfläche 
gespritzt wird. Dazu wird eine Ölbrause verwendet, die aus mehreren identischen Öldüsen und 
einer Zuführbohrung besteht. Diese ist in geringem Abstand zum Zahnrad angeordnet. Die 
Anordnung ist in Abbildung 3.1 dargestellt. 


Die untersuchte Ölbrause besitzt vier Düsen, die gleichmäßig über die Zahnbreite W = 90mm 
verteilt sind. Das Öl wird über die Zuführbohrung zu den Düsen geführt. Diese ist wesentlich 
größer als die Düsenbohrungen, um die Strömungsgeschwindigkeit gering zu halten. Dies sorgt 
dafür, dass an allen Düseneinlässen nahezu identische Anströmbedingungen herrschen und 
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Ölbrause Zuführbohrung 


/ Schwenkachse 


Abbildung 3.1: Charakteristische geometrische Parameter der untersuchten Oleinspritzkiihlung 


die sich ausbildenden Ölstrahlen nur äußerst geringe Unterschiede in ihren Strömungsprofilen 
aufweisen. Der Abstand der Düsen zueinander beträgt 0,25 W, der Abstand der äußeren Düsen zu 
den Stirnseiten des Zahnrads beträgt 0,125W. Die Düsen sind als Bohrungen mit Durchmesser 
d = 1,15mm und einer Länge von / = 7mm ausgeführt. 


Die Strömung des Ölstrahls in der Düse ist entscheidend für den Aufprallvorgang auf dem Zahn- 
rad und somit für den Wärmeübergang. Sie wird vor allem durch die Düsengeometrie beeinflusst, 
wie in Abschnitt 3.2 erläutert wird. Maßgeblichen Einfluss haben der Bohrungsdurchmesser d 
sowie das mit der Länge der Bohrung, /, gebildete Verhältnis L, das im vorliegenden Fall 6,09 
beträgt. 


Aufgrund der gleichmäßigen Aufteilung der Düsen sind für jedes Zahnviertel identische Vor- 
gänge bezüglich des Ölaufpralls und des Wärmeübergangs zu erwarten. Dies ermöglicht die 
Übertragbarkeit auf Fälle, in denen eine periodische Anordnung von Einzeldüsen betrachtet wird. 
Unter der Voraussetzung, dass das Aufeinandertreffen zweier Ölfronten keine stromauf gerichtete 
Rückwirkung hat, gilt dies auch für das Auftreffen eines einzelnen Ölstrahls. Zur Bewertung 
der Übertragbarkeit kann das Verhältnis aus Düsendurchmesser und Breite des periodischen 
Elements herangezogen werden, das in der betrachteten Anordnung 


_ 0,25W 


rg 


= 19,6 (3.1) 
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beträgt. Zur Umsetzung der in Kapitel 4 beschriebenen Methode wird die Periodizität der 
Anordnung ebenfalls ausgenutzt. 


Die Strömung der Luft, die das Zahnrad umgibt, interagiert mit dem Ölstrahl. Neben der Ober- 
flächengeschwindigkeit des Zahnrads, das diese Strömung anregt, ist die Länge der Strecke, 
die der Ölstrahl zwischen Düsenaustritt und Zahnrad zurücklegt, entscheidend für die Interak- 
tion zwischen Luft und Öl (Akin et al. (1975), Krug (2018)). Sie wird bestimmt durch den 
Abstand R zwischen der Schwenkachse der Düse und der Zahnradachse, der Exzentrizität E 
des Düsenaustritts, dem Teilkreisdurchmesser dy des Zahnrads und dem Diisenwinkel a. Der 
Düsenwinkel ist der Winkel zwischen der Bohrungsachse einer Öldüse und der Ebene, die durch 
die Schwenkachse der Ölbrause und die Zahnradachse aufgespannt wird (siehe Abbildung 3.1). 


Eine Übersicht der geometrischen Parameter, die die radiale Öleinspritzkühlung definieren, ist in 
Tabelle 3.1 dargestellt. 


Tabelle 3.1: Geometrische Parameter der untersuchten Einspritzkühlung 


Parameter Symbol Wert 
Düsenradius R 153 mm 
Düsendurchmesser d 1,15 mm 
Länge der Düsenbohrung l 7mm 
Exzentrizität des Diisenauslasses E 10mm 
Zahnbreite W 90 mm 
Teilkreisdurchmesser dw 260 mm 


Die Abmessungen des untersuchten geradverzahnten Zahnrads sind in Tabelle 3.2 zusammenge- 
fasst. Dabei wird mit 1 = 0° ein geradverzahntes Zahnrad eingesetzt, um Gegenlichtaufnahmen 
des Strahlaufpralls zu ermöglichen. 


Tabelle 3.2: Geometrie des untersuchten Zahnrads 


Parameter Symbol Wert 
Zähnezahl Z 65 
Eingriffswinkel Y 25° 
Schrägungswinkel l 0° 
Modul m 4mm 
Profilverschiebung +1 mm 
FuBradius 0,8 mm 
Fußkreisdurchmesser 252 mm 
Wälzkreisdurchmesser dw 260 mm 
Kopfkreisdurchmesser 270 mm 


Das untersuchte Zahnrad wird durch Vorverzahnen und anschlieBendes Schleifen hergestellt. 
Aus Gründen, die sich aus der in Kapitel 4 beschriebenen Methode zur Bestimmung des Wärme- 
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übergangs ergeben, wird die Titanlegierung TiAl6V4 als Zahnradmaterial eingesetzt. Bei dieser 
Kombination aus Herstellverfahren und Material sind mittlere Rauheiten der Flanken von etwa 
Ra = 0,4 — 0,6 um zu erwarten. 


Variationsparameter des Öl-Zahnrad-Wärmeübergangs 


In den meisten Anwendungsfällen in der Luftfahrt sind für Getriebe variable Drehzahlen zu 
erwarten. Dies gilt sowohl für Hilfsgerätegetriebe in Turbofantriebwerken als auch für Leis- 
tungsgetriebe, die in Turboprop- oder Getriebefantriebwerken eingesetzt werden. Auch in der 
vorliegenden Arbeit wird die Drehzahl n des gekühlten Zahnrads daher als Parameter betrachtet 
und variiert. 


Der Ölvolumenstrom V hatte in den Arbeiten von Townsend und Akin (1983), Handschuh (1992) 
und Schober (1983a) einen deutlichen Einfluss auf die gemessenen Zahnradtemperaturen. Zudem 
zeigen die kinematische Betrachtung von Akin et al. (1975) sowie die in Abschnitt 3.3 abgeleitete 
Beschreibung die Abhängigkeit des Strömungsverhaltens vom Ölvolumenstrom auf. Massini 
et al. (2016) bestätigen die Bedeutung des Geschwindigkeitsverhältnisses von Ölstrahl und 
Zahnflanke für das Auftreffen des Ölstrahls. Dies legt es nahe, den Ölvolumenstrom ebenfalls 
als Variationsparameter zu betrachten. 


Bisherige Untersuchungen können unterteilt werden in solche, bei denen Öl in radiale Richtung 
eingespritzt wird, und solche, bei denen das Öl tangential in den Zahnradeingriff eingespritzt 
wird. Bei tangentialer Einspritzung kann die Düse entweder in Drehrichtung vor oder nach 
dem Zahnradeingriff angeordnet sein (engl.: into-mesh und out-of-mesh). Die Anordnung der 
Düse kann durch ihren Einstellwinkel &, der im Vorangegangenen definiert wurde, beschrieben 
werden. Die radiale und tangentiale Anordnung entsprechen dabei den Einstellwinkeln 0° 
und 90°. Fondelli et al. (2015), Keller et al. (2017) und Massini et al. (2017) betrachten den 
Einfluss unterschiedlicher Einstellwinkel numerisch und beobachten deutliche Einflüsse auf 
das Strömungsverhalten. Dies wird durch die kinematische Betrachtung des Einspritzvorgangs 
bestätigt. & wird im Rahmen dieser Arbeit als variabler Parameter betrachtet. 


Aufgrund der starken Temperaturabhängigkeit der Viskosität der typischerweise in der Luftfahrt 
verwendeten Öle ist davon auszugehen, dass auch die Temperatur den Öl-Zahnrad-Wärmeübergang 
beeinflusst. Das Auftreten von laminaren oder turbulenten Ölstrahlen hängt wesentlich von der 
Reynoldszahl ab, die mit der Viskosität gebildet wird. Die Viskosität bzw. die Temperatur des 
Öls 9, wird im Kontext dieser Arbeit ebenfalls als Parameter untersucht und variiert. 


Ein Betriebspunkt des betrachteten Systems wird durch vier variable Parameter beschrieben: 
e Drehzahl n 
e Ölvolumenstrom V 
e Düsenwinkel & 


e Oltemperatur w 


Der untersuchte Wertebereich der obenstehenden Parameter in dieser Arbeit wird in Tabelle 3.3 
zusammengefasst. 


3.2 Düsenströmung 35 


Tabelle 3.3: Wertebereich der Untersuchungen 


Parameter Bereich 


n bis 6000 /min 
V bis 10,41/min 
a -45-45° 

V, 60,80°C 


3.2 Düsenströmung 


Die Gestaltung der vorangegangenen Abschnitt beschriebenen Düse und der gewählte Parame- 
terbereich der Versuche führen zum Auftreten verschiedener Strömungsformen des Ölstrahls. 
Sowohl anliegende als auch abgelöste Ölstrahlen werden beobachtet. Diese unterscheiden sich in 
ihrem Durchmesser, was zu unterschiedlichen Geschwindigkeiten führt, als auch durch den Tur- 
bulenzgrad der Strömung. Um die gemessenen Wärmeübergangsverteilungen sinnvoll auswerten 
zu können, ist die Unterscheidung zwischen diesen Erscheinungsformen notwendig. Im Fol- 
genden werden die zu erwartenden Strömungsformen und ihre charakteristischen Eigenschaften 
wiedergegeben. 


Die Erzeugung eines Ölstrahls mit ausreichender Geschwindigkeit ist einerseits für die Ein- 
spritzkühlung notwendig, um hohe Staudrücke, geringe Grenzschichtdicken und somit hohe 
Wärmeübergänge sicherzustellen. Andererseits setzt auch die Kinematik des Ölstrahls hohe 
Strahlgeschwindigkeiten voraus, um ein tiefes Eindringen des Öls in die Zahnlücke zu ermög- 
lichen. Sowohl das Geschwindigkeitsprofil des Strahls, das aus der Düsengeometrie und der 
eingespritzten Öl€menge resultiert, als auch die Queranströmung durch die Luft, die vom Zahnrad 
mitbewegt wird, sind dafür entscheidend. Die strömungsmechanischen Phänomene, die hierbei 
eine Rolle spielen, werden von verschiedenen Autoren eingehend betrachtet und sollen in dieser 
Arbeit nicht detailliert beleuchtet werden. Eine Übersicht über die Phänomene der Düsenströ- 
mung wird unter anderem von Birouk und Lekic (2009) und von Kuensberg Sarre et al. (1999) 
gegeben. Der Zerfall von querangeströmten Flüssigkeitsstrahlen wird durch Faeth et al. (1995) 
zusammengefasst. 


Schweitzer (1937) beobachtet, dass Flüssigkeitsstrahlen in ruhender Gasumgebung, ähnlich zur 
Rohrströmung, sowohl in laminarer als auch in turbulenter Form auftreten können. Entscheidend 
für die Ausbildung der jeweiligen Strömungsform ist dabei die Reynoldszahl des Strahls, definiert 
als: 


Bas. (3.2) 


Reg 
In Analogie zur Rohrströmung ist unterhalb einer kritischen Reynoldszahl mit laminarer Strö- 
mung des Ölstrahls zu rechnen. Ab der kritischen Reynoldszahl können Störungen in der 
Düsenströmung ungedämpft anwachsen und zu einer vollturbulenten Düsenströmung führen. Im 
Fall von Rohrströmungen liegt die kritische Reynoldszahl im Bereich von etwa Re ~ 2300 (Rotta 
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(1956)). Liegen, durch geeignete Gestaltung der Düse, nur sehr geringe Störungen vor, kann 
sich auch oberhalb der kritischen Reynoldszahl ein laminarer Strahl ausbilden. Entscheidende 
Geometrieeinflüsse sind unter anderem der Einlauf der Düse, die Rauigkeit der Bohrung und die 
Anwesenheit von Krümmungen und Querschnittsänderungen (Giffen und Muraszew (1953)). 
Diese beeinflussen die Bedingungen am Auslass der Düse und somit den anschließenden Zerfall 
des Strahls. Wu et al. (1994) beobachten laminare Strahlen bei Re ~ 2- 10* und einer Düsen- 
geometrie mit L = 7. Bei noch kleineren Verhältnissen l < 4—6 wird auch bei Rex 1 - 10° ein 
laminarer Strahl festgestellt. 

Das Strömungsregime in der Düse hat Auswirkungen auf das Geschwindigkeitsprofil an ihrem 
Austritt und den Durchflussbeiwert, der den notwendigen Öldruck bestimmt. Nach dem Austritt 
aus der Düse interagiert der Ölstrahl mit der ihn umgebenden Luft. Bei einer ausreichend 
langen freien Wegstrecke zerfällt der Strahl zu Tropfen. Dabei treten für laminare und turbulente 
Strahlen unterschiedliche Zerfallsmechanismen auf. Der Zusammenhang zwischen Zerfallslänge 
L des Strahls und der Austrittsgeschwindigkeit u ist qualitativ in Abbildung 3.2 dargestellt. Bei 
geringen Austrittsgeschwindigkeiten und laminarer Strömung dominiert der Rayleigh-Zerfall, 
bei dem die Oberflächenspannung der Flüssigkeit im umgebenden Gas zur Bildung einzelner 
Tropfen führt (Bereich A und B). Mit zunehmender Austrittsgeschwindigkeit steigt der Abstand 
zwischen dem Düsenaustritt und dem Beginn des Zerfalls zunächst an, fällt dann jedoch wieder 
ab. 


Abbildung 3.2: Verlauf der Zerfallslänge eines Ölstrahls in Abhängigkeit der Strömungsge- 
schwindigkeit nach Lin und Reitz (1998) 


Bei weiterer Erhöhung der Austrittsgeschwindigkeit kann die Zerfallslänge weiter absinken, wenn 
der Strahl turbulent wird. In diesem Fall führt die Turbulenz zum Auffächern des Strahls unter 
Verzögerung der mittleren Strahlgeschwindigkeit (Bereich D). Die dabei entstehenden Tropfen 
sind wesentlich kleiner als der Düsendurchmesser. Ist die Strömung in der Düse zusätzlich 
abgelöst, etwa durch eine scharfkantige Einschnürung im Düseneinlass (siehe Abbildung 3.3), 
kann die turbulente Zerfallslänge jedoch auch deutlich ansteigen (Bereich C). 
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Anliegende Strömung Abgelöste Strömung 


--- vena contracta vena contracta 


Panel VASEN DEREN Wiederanlegepunkt 


Abbildung 3.3: Anliegende und vollständig abgelöste Düsenströmungen nach von Kuensberg 
Sarre et al. (1999) 
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3.3 Kinematik des Ölstrahls 


Die Bedeutung der zugeführten Ölmenge für den Wärmehaushalt des Getriebes wird von einigen 
Autoren diskutiert. In der entwicklungsnahen Fachliteratur werden häufig bestimmte, auf Erfah- 
rungswerten beruhende Ölvolumenströme empfohlen, um ausreichende Kühlung sicherzustellen. 
So nennen Winter und Michaelis (1984) einen zur Kühlung notwendigen, auf die Zahnbreite 
bezogenen Mindestölvolumenstrom: 


. ~ 0,08 —0,11/min (3.3) 


Handschuh (1992) bestimmt die Oberflächentemperaturen an Kegelrädern unter Last und Ein- 
spritzschmierung durch einen Fächerstrahl. Die beobachteten Oberflächentemperaturen nehmen 
bei Erhöhung des Ölvolumenstroms ab. Es zeichnet sich ein asymptotischer Verlauf ab, der 
Effekt des Ölvolumenstroms ist mäßig. Diese qualitative Beobachtung wird durch Townsend 
und Akin (1983) für einspritzgekühlte Stirnräder bestätigt. 


Eine systematische Methode zur Bestimmung des Wärmeübergangs in Abhängigkeit der geo- 
metrischen Parameter und Betriebsgrößen wird von den aufgeführten Quellen nicht dargelegt. 
Einige Quellen geben Hinweise auf die zu Grunde liegenden Abhängigkeiten. Die Vermutung 
liegt nahe, dass die kinematischen Beziehungen zwischen bewegter Zahnkontur und eindringen- 
dem Ölstrahl einen erheblichen Einfluss auf die Kühlwirkung des Öls haben. Akin et al. (1975) 
greifen diesen Gedanken auf, indem sie den Einfluss der Eindringtiefe des Ölstrahls in radialer 
Richtung auf den Wärmeübergang betrachten. Die Autoren stellen ein kinematisches Modell des 
Einspritzvorgangs auf, das in Abbildung 3.4 dargestellt ist. 
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Abbildung 3.4: Kinematisches Olaufprallmodell nach Akin et al. (1975) 


Das Modell von Akin et al. (1975) erlaubt die Berechnung der Eindringtiefe unter Vereinfachung 
des Einspritzvorgangs. Dazu wird die Rotationsbewegung des Zahnrads als lineare Bewegung 
modelliert. Es wird weiterhin vorausgesetzt, dass Ol radial eingespritzt wird. Die Geschwindig- 
keit des Olstrahls wird als homogen angenommen und abhingig vom Versorgungsdruck unter 
Annahme einer reibungsfreien Beschleunigung berechnet. Zudem wird der Einfluss der vom 
Zahnrad mitbewegten Luft, die den Ölstrahl quer anströmt, als gering angesehen und daher nicht 
berücksichtigt. Anhand der Ölstrahlgeschwindigkeit und Zahngeschwindigkeit bestimmen die 
Autoren den Eindringwinkel ß, unter dem der Ölstrahl in die Verzahnung eindringt. Townsend 
und Akin (1983) berechnen anhand des Modells von Akin et al. (1975) die Eindringtiefe des 
Ölstrahls für die experimentell untersuchten Bedingungen. Für die Betriebsbedingungen, bei 
denen der Ölstrahl den Zahnfuß erreicht, werden die geringsten Temperaturen und somit die 
größte Kühlwirkung beobachtet. 


Erweiterung des kinematischen Modells zur Berücksichtigung des Diisenwinkels a 


Das Modells von Akin et al. (1975) wird im Folgenden erweitert, um die betrachtete Ein- 
spritzkühlung zu beschreiben. Dazu wird die Relativgeschwindigkeit des Ölstrahls, w., als 
Vektorsumme aus der Geschwindigkeit des Zahnrads, u,, und der Ausbreitungsgeschwindigkeit 
des Ölstrahls in Richtung der Düsenachse, üs, betrachtet. Um den Einfluss des Düsenwinkels 
a zu berücksichtigen, wird u, unter dem Winkel & angesetzt. Eine schematische Darstellung 
dieses Modells ist in Abbildung 3.5 gegeben. 
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Abbildung 3.5: Kinematik des Ölstrahls unter Berücksichtigung des Düsenwinkels 


Zur Interpretation der Ergebnisse dieser Arbeit wird der Winkel ß aus Abbildung 3.5 verwendet, 
der in Anlehnung an die Literatur im Folgenden als Aufprallwinkel bezeichnet wird. Zur Veran- 
schaulichung des Wertebereichs von ß sind in Abbildung 3.6 vier signifikante Aufprallwinkel 
eingezeichnet. Bei Werten für ß von unter 28,7° erreicht der Ölstrahl auch die strömungs- 
abgewandte Flanke. Bei einem Einfall des Ölstrahls mit ß = 33,8° erreicht der Ölstrahl die 
Zahnfußrundung und damit die gesamte strömungszugewandte Zahnflanke. Bei einem Aufprall- 
winkel von ß = 60,3° erreicht der Ölstrahl etwa die Hälfte der Zahnhöhe. Der Strahl triff dabei 
an seiner Front orthogonal auf die Zahnflanke auf. Ab dem Punkt des orthogonalen Auftreffens 
und bei größeren Werten für B besitzt das Öl keine Geschwindigkeitskomponente in Richtung 
des Zahnfußes. Nach dem Aufprall ist in diesem Fall nicht mit einem weiteren Eindringen des 
Öls in Richtung des Zahnfußes zu rechnen. 


Abbildung 3.6: Diskrete Werte des Aufprallwinkels an vier signifikanten Stellen der Flanke 


Um eine vektorielle Betrachtung des Ölaufpralls anstellen zu können, muss die Verteilung der 
axialen Geschwindigkeit des Ölstrahls als homogen betrachtet werden. Zur Berechnung der 
axialen Geschwindigkeit wird in der Literatur meist die mittlere Geschwindigkeit in der Düse 
V 
Us 


ae (3.4) 
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aus Volumenstrom V und Düsenquerschnitt Ag herangezogen. Im vorliegenden Fall ist weder 
bei laminaren noch bei turbulenten Strahlen eine homogene Geschwindigkeitsverteilung zu 
erwarten. Zusätzlich kann der Strahldurchmesser bei Ablösung der Strömung in der Düse deutlich 
reduziert sein, wodurch aus Kontinuitätsgründen höhere Strahlgeschwindigkeiten auftreten. Das 
kinematische Modell liefert unter Verwendung der Geschwindigkeit u, nur Näherungswerte. 
Eine genauere Beschreibung erfordert die Charakterisierung der Düsenströmung der eingesetzten 
Düse, die im Rahmen dieser Arbeit nicht erfolgt. Arbeiten wie die von Townsend und Akin 
(1983) und die Betrachtung der Ergebnisse in Kapitel 6 zeigen jedoch, dass die vereinfachte 
Betrachtung zur Beurteilung der auftretenden Wärmeübergangsverteilung bereits ausreichend 
ist. 

Unter den getroffenen Annahmen sind der Aufprallwinkel B und der spezifische Auftreffimpuls n 
nur von den geometrischen Verhältnissen und nicht von absoluten Geschwindigkeiten abhängig. 
Daher kann zur Vereinfachung der Betrachtung das Geschwindigkeitsverhältnis 


Y= — (3.5) 


aus Ölstrahlgeschwindigkeit und Tangentialgeschwindigkeit des Zahnrads gebildet werden. Für 
den Winkel ß gilt: 


ß = atan (=) : (3.6) 


Urr 


Im geschilderten Modell ist der Winkel œ positiv, wenn die Düse in die Drehrichtung des 
Zahnrads gerichtet ist. Somit gilt für die Tangentialkomponente der Relativgeschwindigkeit 


Urt = Uz — Us,t (3.7) 


und somit 


p = atan (= = fs) ; (3.8) 
Us,r 


Durch Bestimmen und Einsetzen der Tangential- und Radialkomponenten der Strahlgeschwin- 
digkeit kann der Term in Gl. 3.8 umgeformt werden zu 


P = atan (Sa) f (3.9) 


cos(&) Us 


Durch Einsetzen des Geschwindigkeitsverhältnisses © kann B bestimmt werden durch 
1 . 
= — sin(a 
p = atan (o) ; (3.10) 


Der Zusammenhang zwischen Ł, œ und ß aus Gleichung 3.10 ist in Abbildung 3.7 dargestellt. 
Bei £ = 1 führt eine Änderung des Düsenwinkels zu einer linearen Änderung des Aufprallwinkels. 
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Für die radiale Einspritzung mit & = 0° beträgt der Aufprallwinkel B = 45°. Bei © < 1 existiert 
ein Minimum des Aufprallwinkels, mit: 


Bmin = arcsin (X) (3.11) 


bei dem Diisenwinkel 

Op) = 90° — arcsin(Z). (3.12) 
Das Minimum des Aufprallwinkels ist in Abbildung 3.7 als gestrichelte Linie gekennzeichnet. 
Eine Erhöhung des Düsenwinkels zu höheren Werten führt bei & > ag, zu flacheren Aufprall- 
winkeln und zusätzlich geringeren Relativgeschwindigkeiten. Daher kann Œg als obere Grenze 
für den Düsenwinkel gesehen werden, um die Kühlwirkung zu maximieren. 


Bei Geschwindigkeitsverhältnissen £ > 1 existiert diese Grenze nicht und der Aufprallwinkel 
2 
sinkt monoton mit & und SE <0. 


80 


60 
40 


Abbildung 3.7: Einfluss des Diisenwinkels auf den Aufprallwinkel B fiir verschiedene Geschwin- 
digkeitsverhältnisse & 


Aus der Übersicht zur Prallkühlung in Abschnitt 2.4 wird deutlich, dass der Impuls des Prallstrahls 
einen entscheidenden Einfluss auf den Wärmeübergang hat. Auch im vorliegenden Fall kann 
der relative Impuls herangezogen werden, um das beobachtete Verhalten des Wärmeübergangs 
zu interpretieren. Analog zur Bestimmung des Aufprallwinkels kann u, aus der kinematischen 
Betrachtung abgeleitet werden. Um den Impuls in Abhängigkeit von X und & zu bestimmen, 
wird der spezifische Impuls gebildet: 


“ = /1+2?-2%sin(o). (3.13) 
Uz 

Für die Herleitung des Zusammenhangs sei auf den Anhang verwiesen. Aus der Darstellung des 
spezifischen Impulses aus Gleichung 3.13 in Abbildung 3.8 ist ersichtlich, dass eine Erhöhung 
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Abbildung 3.8: Einfluss des Düsenwinkels auf den Relativimpuls a = Y1+2?- 2¥sin(q) für 
verschiedene Geschwindigkeitsverhältnisse Ł 


des Düsenwinkels stets zu einer Reduzierung der spezifischen Relativgeschwindigkeit führt. Je 
geringer das Geschwindigkeitsverhältnis ÈX ist, desto geringer ist der Einfluss des Diisenwinkels 
auf den Impuls. Dies ist für den Grenzfall £ — 0 einfach nachvollziehbar, bei dem der Ölstrahl 
keine Eigengeschwindigkeit und somit stets die Relativgeschwindigkeit u, = uz besitzt. 


Aus der Betrachtung des Aufprallwinkels wurde deutlich, dass ein positiver Düsenwinkel, 
also Einspritzung in Richtung der Laufrichtung des Zahnrads, zu einem tieferen Eindringen 
des Ölstrahls führen kann. Tiefes Eindringen wird in der Literatur mit vergleichsweise hoher 
Kühlwirkung in Verbindung gebracht und ist somit für die effiziente Nutzung des Kühlmittels 
vorteilhaft. Die Betrachtung des spezifischen relativen Impulses zeigt jedoch, dass der Impuls des 
Ölstrahls mit dem Düsenwinkel absinkt. Unter Verwendung der Erkenntnisse zur Prallkühlung 
kann dies als abträglich für hohe Wärmeübergänge betrachtet werden. Somit liegen zwei 
gegenläufige Einflüsse vor, die den Öl-Zahnrad-Wärmeübergang beeinflussen. Eine Bewertung 
der jeweiligen Einflüsse ist durch das vereinfachte Modell nicht möglich und kann erst durch die 
Interpretation der Ergebnisse in Kapitel 6 erfolgen. 
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4 Experimentelle Methode zur Bestimmung des 
Öl-Zahnrad-Wärmeübergangs 


In Kapitel 2 wird dargelegt, dass die genaue Kenntnis des Wärmeübergangs an einspritzgekühlten 
Zahnrädern zur Entwicklung von effizienten hochbelasteten Luftfahrtgetrieben unerlässlich ist. 
Aus der Analyse der Literatur geht hervor, dass Daten bisheriger Untersuchungen ungenügend 
sind, da der Wärmeübergang durch Kühlung und die Reibung nicht voneinander getrennt werden 
können. Die experimentelle Methode, mit der im Rahmen dieser Arbeit diese Einschränkung 
aufgehoben und erstmals der Öl-Zahnrad-Wärmeübergang ermittelt wird, soll im Folgenden 
vorgestellt werden. Zunächst wird die Erzeugung des Wärmeübergangs für ein Einzelzahnrad 
und die numerische Methode zur Auswertung der gemessenen Temperaturen beschrieben. An- 
schließend wird die experimentelle Umsetzung der Methode geschildert. Abschließend wird 
die Messgenauigkeit der Methode für die untersuchte Geometrie und Messstellenanordnung 
statistisch analysiert. 


4.1 Erzeugung des Wärmeübergangs 


In Untersuchungen wie der von Townsend und Akin (1981), Leoni (1991) oder Handschuh 
(1992) werden experimentelle Daten zum Wärmeübergang an einspritzgekühlten Zahnrädern 
erhoben. In allen bekannten Fällen wird dabei im Zahnradeingriff Leistung übertragen und 
Verlustwärme freigesetzt. Durch verschiedene Methoden, z.B. Pyrometrie, Einsatz von Thermo- 
elementen, werden Temperaturen auf der Zahnradoberfläche oder im Zahnradinneren gemessen. 
Die anschließende Analyse der Temperaturen liefert qualitative Informationen zum Einfluss der 
verschiedenen Kühlparameter wie Ölvolumenstrom oder Ölstrahlkonfiguration. Die Bilanzierung 
der Wärmequellen und -senken ist jedoch nicht möglich, da sowohl die entstehenden Verluste 
als auch die Randbedingungen nicht eindeutig und örtlich aufgelöst bestimmt werden können. 
Somit kann der Wärmeübergang, der durch die Einspritzkühlung hervorgerufen wird, durch diese 
Untersuchungen nur als stark fehlerbehaftete und globale Größe quantifiziert werden. 


In der vorliegenden Arbeit wird der Wärmeübergang durch Einspritzkühlung isoliert von anderen 
Wärmeströmen gemessen. Dazu wird im Folgenden eine experimentelle Methode vorgestellt, 
mit der der zur Messung des Wärmeübergangs erforderliche Wärmestrom nicht durch die Rei- 
bungswärme eines Zahnradeingriffs erzeugt wird. Vielmehr wird ein einzelnes hohles Zahnrad 
eingesetzt, auf dessen Außenmantel heißes Öl gespritzt und dessen Innenkontur durch Prall- 
kühlung mit Luft gekühlt wird. Dadurch stellt sich ein Temperaturfeld über das Zahnrad ein, 
das von den Wärmeübergängen an Außen- und Innenfläche, h(X) und hx, sowie der Geometrie 
und der Wärmeleitfähigkeit A des Zahnradmaterials abhängig ist. Durch die Bestimmung des 
Temperaturfelds kann die Größe h(x) bestimmt werden. Im Folgenden soll dargelegt werden, 
unter welchen Umständen dies möglich ist. 


In Abbildung 4.1 ist ein Sektor eines hohlen Zahnrads dargestellt, an dessen in Drehrichtung 
orientierten Zahnflanken Wärme eingetragen und dessen Innenmantel gekühlt wird. Die Wärme- 


46 Experimentelle Methode zur Bestimmung des Öl-Zahnrad-Wärmeübergangs 


übergangskoeffizienten der Zahnflanken und der Zahnradinnenkontur wurden dabei als homogen 
angenommen. Die der Drehrichtung abgewandten Flanken, die Zahnspitze und die Zahnfuß- 
rundung sowie die Seitenflächen wurden zunächst als adiabate Wände betrachtet. Das vordere 
Viertel des Sektors ist ausgeschnitten, um den qualitativen thermischen Zustand innerhalb des 
Festkörpers darzustellen. Der Verlauf der Isothermen ist durch durchgezogene Linien, die 
Wärmestromlinien sind als Strichlinien dargestellt. 


Abbildung 4.1: Qualitative Temperaturverteilung eines Zahnradsektors mit vereinfachten thermi- 
schen Randbedingungen 


Trotz der angenommenen homogenen Randbedingungen stellt sich im Festkörper ein nicht 
triviales Temperaturfeld ein. Das Temperaturfeld ist abhängig von den Randbedingungen, der 
Geometrie des Zahnrads und dem gewählten Material des Festkörpers. Für die Wärmeleitung im 
Prüfkörper gilt allgemein: 


d 

77 Ot) — aA (Kt) = qy (X,t) (4.1) 
mit der Temperatur v, der Temperaturleitfähigkeit a und der volumetrischen Wärmestromdichte 
qy an der Stelle X zum Zeitpunkt t. Die Wärmestromdichte ist in Gleichung 4.1 auf die spezifische 
Wärmekapazität bezogen. In den vorliegend geschilderten Experimenten werden stets stationäre 


Temperaturzustände des Versuchszahnrads betrachtet, somit gilt: 


Zusätzlich gilt für die im Rahmen dieser Arbeit betrachteten Experimente Quellen- und Senken- 
freiheit, da im Zahnrad keine Wärme freigesetzt wird. Es gilt: 


qv (x) = 0. (4.3) 


Dadurch wird Gleichung 4.1 zu: 
— aAd(xX) = 0. (4.4) 
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Gleichung 4.4 ist eine elliptische partielle Differentialgleichung zweiter Ordnung, die Laplace- 
Gleichung genannt wird. Durch die Quellen- und Senkenfreiheit, die im vorliegenden Fall stets 
gilt, ist die Laplace-Gleichung eindeutig lésbar, sofern ihre Randwerte bekannt sind und die 
Temperaturleitfahigkeit als konstant betrachtet werden kann. Dies bedeutet, dass die Kennt- 
nis der Temperaturverteilung am Rand des Gebiets, also der Festkörperoberfläche, ausreicht, 
um die Verteilung der Wärmestromdichte im Gebiet zu bestimmen. Die Ausnutzung dieses 
Zusammenhangs ist der Kerngedanke der beschriebenen Methode. 


Durch die oben beschriebenen Zusammenhänge kann der Wärmestrom im Zahnrad eindeutig 
durch Messung der Oberflächentemperatur des Zahnrads bestimmt werden. Da an der Oberfläche 
des Zahnrads nur Wärmeleitung auftritt und keine Wärme freigesetzt wird muss gelten: 


Gkonv (X) lo = Gkond (X) llo- (4.5) 


Der Wärmestrom an das Fluid gkonv (X) entspricht demnach dem Wärmestrom an der Oberfläche 
des Festkörpers gGxona(X). Der Wärmeübergangskoeffizient kann bestimmt werden als: 


h(&) = Ion) (4.6) 


(4.7) 


mit der Warmeleitfahigkeit des Öls, A, und der charakteristischen Länge, für die hier der Modul 
des Zahnrads, m, eingesetzt wird. 


Im Vorangegangenen wird stets davon ausgegangen, dass ein Temperaturgradient vorliegt, der 
zu einem Wärmeübergang vom eingespritzten Öl in das Zahnrad führt. In der Anwendung in 
einem Höchstleistungsgetriebe ist eine entgegengesetzte Wärmestromrichtung zu erwarten. Die 
Umkehrung der Wärmestromrichtung ist aufgrund der experimentellen Umsetzung im Rahmen 
dieser Arbeit erforderlich, wie in Kapitel 5 erläutert wird. Für die geschilderte Berechnung des 
Wärmestroms ist die Richtung des Wärmestroms unerheblich. Aufgrund der hervorgerufenen 
Temperaturänderung im Öl beeinflusst die Richtung des Wärmestroms jedoch die Strömungs- 
vorgänge durch Prallkühlung. Im Fall von Fluiden, deren Viskosität temperaturunabhängig 
ist, ist der Wärmeübergangskoeffizient temperaturunabhängig, wodurch die Richtung des Wär- 
meübergangs keine Rolle spielt. Diese Annahme kann für das verwendete Öl nicht getroffen 
werden. Jedoch wird durch die Auslegung des Versuchszahnrads im Folgenden sichergestellt, 
dass nur vergleichsweise geringe Temperaturänderungen auftreten und so die Übertragbarkeit 
der Ergebnisse auf den Anwendungsfall gewahrt bleibt. 


Zur experimentellen Bestimmung des Öl-Zahnrad-Wärmeübergangs ist ein Versuchszahnrad 
erforderlich, dessen Temperaturen im Versuchsbetrieb unter definierten Betriebsbedingungen 
gemessen werden. Die Geometrie und der Wärmeleitwert des eingesetzten Materials müssen 
bekannt sein. Sind diese Bedingungen erfüllt, kann der Wärmeübergang durch die Lösung von 
Gleichung 4.4 bestimmt werden. Aufgrund der Geometrie des Zahnrads, die sich einer analyti- 
schen Beschreibung entzieht, muss die Lösung numerisch bestimmt werden. Die Beschreibung 
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des Lösungsverfahrens erfolgt im folgenden Kapitel. Die Geometrie des Versuchszahnrads ist in 
Kapitel 4.4 beschrieben. Auf den Prüfstand, mit dem der Zahnradwärmeübergang unter Einfluss 
variabler Betriebsbedingungen untersucht wird, wird in Kapitel 5 eingegangen. 


4.2 Modellvereinfachung zur Erhöhung der Auflösung 


In vielen Wärmeübergangsuntersuchungen werden hochaufgelöste flächige Temperaturmessun- 
gen durch Thermographie durchgeführt, um die Temperaturverteilung von Oberflächen und somit 
die thermischen Randbedingungen zu bestimmen (u.a. Kneer (2017), Ochs (2011), Pütz (2018)). 
Die örtlich hoch aufgelöste Erfassung der Oberflächentemperaturen des Zahnrads ist jedoch 
nicht möglich. Im vorliegenden Fall verhindern das beteiligte Öl, die hohen Oberflächenge- 
schwindigkeiten und die eingeschränkte optische Zugänglichkeit den Einsatz von Thermographie. 
Die verfügbaren Öle weisen bereits ab geringen Filmdicken Emissionsgrade von etwa € = 1 
auf, weshalb eine infrarotthermographische Messung der Temperatur an den unter dem Ölfilm 
befindlichen Oberflächen verhindert wird. Zudem ist bei hohen Drehzahlen Öl in der Form von 
Tropfen und Nebel in der Luft vorhanden, wodurch keine sinnvolle optische Betrachtung des 
Vorgangs mehr möglich ist. 


Für die Bestimmung der Oberflächentemperaturen in dieser Arbeit werden daher Temperaturen 
durch den Einsatz von Thermoelementen auf dem rotierenden Zahnrad und drahtloser Daten- 
übertragung gemessen. Optische Zugänglichkeit und zeitlich hochaufgelöste Messungen sind 
hierfür nicht erforderlich. Allerdings ist die Anzahl an Temperaturmesspunkten durch die gegen- 
wärtig verfügbaren Telemetriesysteme eingeschränkt. Um dennoch örtliche Effekte auflösen und 
beobachten zu Können, wird diesem Umstand durch Ausnutzung von Symmetrie- und Periodi- 
zitätseigenschaften sowie zielgerichtete Gestaltung des Versuchszahnrads Rechnung getragen. 
Dazu werden folgende Modellvereinfachungen getroffen: 


e Periodizität: Nahezu alle Thermalmodelle in der Literatur nutzen die geometrische Periodi- 
zität von Zahnrädern, um den Modellierungsaufwand auf einen Einzelzahn zu beschränken. 
Diese Eigenschaft kann auch messtechnisch ausgenutzt werden. Die Schnittebenen eines 
periodischen Elements können in einfachsten Fällen als adiabate Wand modelliert werden. 
Wärmeströme zwischen zwei Zahnelementen sind jedoch grundsätzlich möglich und kön- 
nen durch die Implementierung periodischer Ränder berücksichtigt werden. Dabei gilt an 
der Schnittebene: 


ga(X) = —ġB(¥) (4.8) 
sowie 


Ba (x) = vB (X). (4.9) 


Anschaulich bedeutet die Forderung aus Gleichung 4.8, dass an korrespondierenden Punk- 
ten x auf den periodischen Rändern A und B die gleichen Wärmeströme ein- und ausfließen 
sowie gleiche Temperaturen herrschen (Gleichung 4.9). Im Fall des geradverzahnten Zahn- 
rads bieten sich Schnittebenen konstanten Drehwinkels an (siehe Abbildung 4.2). 
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Symmetrie: Die in dieser Untersuchung eingesetzte Ölbrause weist vier identische Einzel- 
düsen auf, die in regelmäßigen Abständen angeordnet sind. Aufgrund dieser Anordnung 
sind in axialer Richtung periodische Wärmeübergänge zu erwarten. Das Zahnelement kann 
daher zunächst zu einem Viertel mit periodischen Rändern vereinfacht werden. Durch die 
Annahme, dass der Wärmeübergang innerhalb eines Viertelzahns zudem symmetrisch ist 
(siehe Strichlinie in Abbildung 4.2), werden die periodischen Ränder zu Symmetrieebe- 
nen. Symmetrieebenen des Wärmeübergangs lassen sich numerisch als adiabate Wände 
abbilden. 


Biotzahl: Das Verhältnis aus Wärmeübergangskoeffizient und Wärmeleitfähigkeit des 
Zahnrads ist entscheidend für die Temperaturverteilung innerhalb des Zahnrads. Dabei 
sind sowohl zu hohe als auch zu niedrige Verhältnisse ungünstig, da sie die Messgenau- 
igkeit negativ beeinträchtigen. Durch geeignete Wahl des Zahnradmaterials, auf die in 
Abschnitt 4.4 eingegangen wird, kann sichergestellt werden, dass trotz weniger diskreter 
Messstellen Wärmeübergänge zwischen den Messpunkten berücksichtigt werden. 


Symmetrieebenen 


Ölstrahlen 


Stirnwand 


I periodischer Rand B 


Achtelzahn 


periodischer Rand A 


Abbildung 4.2: Vereinfachung des Zahnrads als Viertelzahn - Symmetrieebenen und periodische 
Ränder 


Durch die vorangegangen Vereinfachungen ist es möglich, den Zahnradwärmeübergang durch 


die Messung diskreter Oberflächentemperaturen auf einem Achtelzahn zu bestimmen. Dadurch 
können die Messpunkte, die aus technischen Gründen in ihrer Anzahl beschränkt sind, dichter 
verteilt werden. Der Wärmeübergang an den Stirnwänden des Zahnrads blieb bisher unberück- 
sichtigt. Anhand von Abschätzungen durch Korrelationen lässt sich zeigen, dass das Niveau 
des Wärmeübergangs der Stirnseiten weit unterhalb des erwarteten Wärmeübergangs durch 
Ölkühlung liegt und somit nur einen Einfluss im Randbereich des Zahnrads haben wird. Dies ent- 
spricht der Beobachtung von Patir und Cheng (1979). Um den Einfluss der Seitenwände auf die 
gemessenen Wärmeströme gering zu halten, wird im vorliegenden Fall einer der innenliegenden 
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Achtelzähne mit Thermoelementen instrumentiert. Um die Annahme der Symmetrie innerhalb 
des Zahnviertels experimentell überprüfen zu Können, wird der benachbarte Achtelzahn ebenfalls 
mit Messpunkten versehen. Auf die Verteilung der Messpunkte wird in Abschnitt 4.4 genauer 
eingegangen. 


4.3 Iterative Berechnung des Wärmeübergangs 


Das in Abschnitt 4.1 beschriebene Randwertproblem kann dann gelöst werden, wenn die Rand- 
bedingungen vollständig bestimmt sind. Dies setzt voraus, dass die Temperatur jeden Punktes 
an der Berandung des Festkörpers, an dem keine Neumann-Randbedingung gilt, bekannt ist. 
Die Temperaturen an der Oberfläche des Zahnrads werden jedoch im Rahmen dieser Arbeit an 
diskreten Punkten bestimmt. Die Temperaturverteilung zwischen diesen Punkten ist zunächst 
nicht bekannt. Aus Abbildung 4.1 wird deutlich, dass die Temperaturverteilungen, die sich über 
das Zahnrad einstellen, aufgrund der Geometrie des Zahnrads komplex sind. Wird diese Tempe- 
raturverteilung durch lineare Interpolation angenähert und das Randwertproblem gelöst, ergeben 
sich unphysikalische Lösungen des Wärmeübergangs. Dies liegt daran, dass die Temperatur- 
verteilung im Fall der linearen Interpolation an den Stützstellen nicht stetig differenzierbar ist. 
Aber auch die Splineinterpolation „zwingt” die Finite-Elemente-Methode zu unphysikalischen 
Lösungen, da ein stückweise polynomieller Temperaturverlauf nicht den sich in der Realität 
einstellenden Verläufen entspricht. 


Physikalisch sinnvolle Lösungen können bestimmt werden, indem das Randwertproblem iterativ 
gelöst wird. Um die Interpolation von Temperaturwerten zu vermeiden, werden anstatt von 
Temperaturwerten Wärmeübergangskoeffizienten als Randbedingung verwendet. Das gemischte 
Randwertproblem wird so in ein Neumann-Randwertproblem überführt, dessen Lösung ebenfalls 
eindeutig ist. Da nun h(X), der a priori nicht bekannt ist, sowohl Randbedingung als auch 
Lösung des Wärmeleitungsproblems ist, muss zur iterativen Lösung ein Gütekriterium und eine 
Iterationsvorschrift gefunden werden. Die gesuchte Lösung erfüllt die Bedingung 


of = oM., (4.10) 


Anschaulich bedeutet die Forderung aus Gleichung 4.10, dass die Temperaturen der numerischen 
Lösung V” den gemessenen Temperaturen 0” an der Festkörperoberfläche entsprechen sollen. 
Daraus lässt sich unmittelbar das Gütekriterium der Iteration ableiten: 


AD! = oF — 9M, (4.11) 
In der Arbeit von Stripf (2007) wird zur Losung folgender Iterationsansatz vorgestellt: 


0 


ar age — S i 
i i 9, — 07 A 


(4.12) 
Voruntersuchungen haben gezeigt, dass der Ansatz nach Gleichung 4.12 im vorliegenden Fall 
zu unbefriedigendem Konvergenzverhalten führt. In den Wärmeübergangsmessungen von 
Stripf (2007) sind die Temperaturmesspunkte vergleichsweise dicht verteilt, der Verlauf der 
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Wärmeübergangskoeffizienten unterliegt geringen lokalen Schwankungen. Im vorliegenden Fall 
können jedoch unter bestimmten Bedingungen örtlich sehr begrenzte hohe Wärmeübergänge 
auftreten. Daher wird folgender Iterationsansatz abgeleitet: 

AV? 


BH — k +a, (4.13) 
W — 0” 


mit dem arithmetischen Mittel der Wärmeübergangskoeffizienten 


es a, 
Wen (4.14) 
z 
und der Temperaturen 
> Z mM 
ga ELL, (4.15) 
z 


mit der Anzahl der Messpunkte z. 


Die Implementierung der beschriebenen iterativen Berechnung der Wärmeübergangskoeffizien- 
ten ist in Abbildung 4.3 dargestellt. Zur Initialisierung der Berechnung werden Startwerte h? 
an den Messpunkten vorgegeben. Der Wärmeübergang an der Zahnradaußenkontur wird auf 
Grundlage der Startwerte h? interpoliert. Unter Verwendung dieser Startverteilung des Wärme- 
übergangs, der gemessenen Temperaturen am Innenradius des Zahnrads, der Zahnradgeometrie 
sowie der Wärmeleitfähigkeit des Zahnrads wird der Wärmetransport im Zahnrad simuliert. 
Aus der berechneten Temperaturverteilung werden die Temperaturen an den Messpunkten, 0’, 
extrahiert. Durch Subtraktion mit den gemessenen Messpunkttemperaturen V% kann die Tempe- 
raturabweichung Av” bestimmt werden. Mit der Iterationsvorschrift (Gleichung 4.13) können 
anschließend korrigierte Wärmeströme berechnet werden. Bei Unterschreiten einer Schwelle 
für AV” von 0,25 K für alle i werden die zuletzt verwendeten Wärmeübergangskoeffizienten h? 
und die daraus resultierende, interpolierte Wärmeübergangsverteilung als Berechnungsergebnis 
verwendet. 
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Startwerte h? O Interpolation h(x) FEM-Modell 
Sampling 3? 0" (x) 
ntl — pe pn LH Experiment 
0-0" 
nein 
Ad!" <0,25K? Ag? = on — 9M or 


Interpolation h 


Ergebnis 


Abbildung 4.3: Schematische Darstellung der iterativen Berechnung des Wärmeübergangs aus 
diskreten Temperaturmesswerten 


4.4 Experimentelle Umsetzung 


Zur experimentellen Umsetzung der im Vorangegangenen beschriebenen Methode muss ein 
Versuchszahnrad mit diskreten Messpunkten der Temperatur definiert werden. Die Gestaltung 
des Versuchszahnrads, das im Rahmen dieser Arbeit eingesetzt wird, sowie die Instrumentierung 
dieses Zahnrads sollen im Folgenden geschildert werden. 


4.4.1 Optimierte Gestaltung eines Versuchszahnrads 


Die Wahl des Prüfkörpermaterials hat maßgeblichen Einfluss auf die Messgenauigkeit des 
Wärmeübergangskoeffizienten. Die Außenkontur des Prüfkörpers ist durch die Definition der 
Zahnradgeometrie und der Düsenanordnung in Kapitel 3 bestimmt. Die Dicke des Zahnrad- 
mantels kann jedoch geeignet gewählt werden, um günstige Voraussetzungen für eine möglichst 
geringe Messunsicherheit zu schaffen. Dabei muss ein Kompromiss getroffen werden, um mehre- 
re Anforderungen zu erfüllen. Die Viskosität des verwendeten Öls ist stark temperaturabhängig. 
Eine zu große Abkühlung des Öls durch den Wärmeübergang am Zahnrad ist einerseits zu 
vermeiden, um die Strömungsvorgänge nur wenig zu beeinflussen. Die Temperaturdifferenz 
zwischen Öl und Zahnradoberfläche geht in den Nenner des Wärmeübergangskoeffizienten ein. 
Daher sind andererseits Temperaturdifferenzen von mehreren Kelvin notwendig, um den Einfluss 
der Unsicherheit der Temperaturmessung gering zu halten. Zudem sind bei ungünstiger Wahl des 
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Materials lokale Wärmeübergänge nicht messbar. Dies gilt sowohl für Materialien zu geringer 
als auch zu hoher Wärmeleitfähigkeit. 


Für die anschauliche Betrachtung des Zusammenhangs kann die Biotzahl herangezogen werden, 


definiert als: 
5 hLpi 


À ri 
mit dem Wärmeübergangskoeffizienten h, der charakteristischen Länge Lgi und der Wärmeleitfä- 
higkeit A. Für Werkstoffe geringer Wärmeleitfähigkeit nähern sich die Oberflächentemperaturen 
des Prüfkörpers für Bi — © der Temperatur der Umgebung bzw. des Öls. Damit werden die 


Bi 


(4.16) 


gemessenen Temperaturen für hohe Biotzahlen unabhängig vom Wärmeübergang, der Einfluss 
des stets vorhanden Messfehlers steigt gegen unendlich. Im entgegengesetzten Fall von Werk- 
stoffen sehr hoher Wärmeleitfähigkeit wird das Temperaturfeld des Prüfkörpers für Bi + 0 an 
jeder Stelle den gleichen Wert annehmen, da für A — œ auch kleinste Temperaturgradienten 
zu einem unendlich hohen Ausgleichswärmestrom führen würden. Auch in diesem Fall kann 
der lokale Wärmeübergang nicht sinnvoll gemessen werden. Jacobsen (1987) und Glahn (1995) 
kommen aufgrund ähnlicher Überlegungen für den Wärmeübergang in Labyrinthdichtungen 
und Lagerkammern zu dem Schluss, dass die Wärmeleitfähigkeit des Prüfkörpers und seine 
Dimensionierung so gewählt werden muss, dass Bi ~ 1 gilt. 


Im vorliegenden Fall kann für Lg; der Abstand der Messpunkte eingesetzt werden, deren An- 
ordnung im folgenden Abschnitt geschildert wird und etwa 4 mm beträgt. Die Forderung nach 
Bi = 1 erfordert die Kenntnis des Wärmeübergangs an der Oberfläche des Zahnrads. Dieser ist 
im Vorfeld nicht bekannt, kann jedoch durch die Korrelation von DeWinter und Blok (1974) 
sowie aufgrund der in Abschnitt 2.4 geschilderten Untersuchungen des Wärmeübergangs durch 
Prallkühlung abgeschätzt werden. Daraus resultieren abgeschätzte mittlere Wärmeübergänge 
in der Größenordnung von 1000 W/(m?K). Im vorliegenden Fall wurde die Titanlegierung 
TiAl6V4 als Material gewählt, die mit A = 7W/(mK) und Bi = 0,57 die Forderung erfüllt. 
Durch die sehr geringe Temperaturabhängigkeit der Wärmeleitfähigkeit von TiAl6V4 ist zudem 
die Gültigkeit von Gleichung 4.4 sichergestellt. 


Um die Abhängigkeit der mittleren Temperatur von der Wärmeleitfähigkeit des Versuchszahnrads 
und dem Wärmeübergangskoeffizienten darzustellen, kann der Einzelzahn als eindimensional 
betrachtet werden (siehe Abbildung 4.4). Mit der Fourier’schen Wärmeleitungsgleichung, dem 
Wärmeübergang durch Kühlung am Innenradius des Zahnrads und dem abgeschätzten Öl- 
Zahnrad-Wärmeübergang ergibt sich die Temperatur der Zahnaußenfläche: 

Vs — OK 


hy + he 1 
mit der Oltemperatur ®,, der Temperatur der Kühlluft und dem Wärmeübergangskoeffizienten 
auf der Innenseite Vg und hx, der Dicke des eindimensionalen Vergleichskörpers L sowie dem 
Wärmeübergangskoeffizienten durch Einspritzkühlung, h. 


Zur Kühlung der Innenseite des Zahnrads wird ein Runddüsenfeld eingesetzt, das auf dem 
Außenmantel eines Zylinders angeordnet ist. Unter Verwendung der umfangreichen Literatur zur 
Berechnung des Wärmeübergangskoeffizienten für Runddüsenfelder kann der Wärmeübergang 
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auf der Innenseite des Zahnrads abgeschätzt werden und beträgt etwa 350W/(m?K) (VDI 
(2013)). 


Die Oberflächentemperatur nach Gleichung 4.17 ist in Abbildung 4.5 für verschiedene Prüfkör- 
perdicken dargestellt. Dabei wurde eine Öltemperatur von 0, = 80°C und Kühllufttemperatur 
von Ux = 20°C angenommen. Bei geringen Dicken des Prüfkörpers sind die Temperaturdifferen- 
zen zwischen einspritzgekühlter und luftgekühlter Oberfläche gering, die Oberflächentemperatur 
ist niedrig. Dadurch sind größere Viskositätsänderungen des Öls und somit geänderte Strömungs- 
vorgänge zu erwarten. Bei einer Prüfkörperdicke von 50 mm ist die Oberflächentemperatur im 
erwarteten Bereich wesentlich schwächer abhängig von h, was die Messunsicherheit erhöht (siehe 
rot markierter Bereich in Abbildung 4.5). Aus diesen und aus konstruktiven Gründen wurde der 
Innendurchmesser des hohlen Zahnrads zu d; = 230mm gewählt, was einer Prüfkörperdicke von 
15mm entspricht. 


4.4.2 Instrumentierung des Versuchszahnrads 


Um die Anordnung der Messpunkte und die Ergebnisse dieser Arbeit anschaulich darzustellen, 
wird ein Koordinatensystem eingeführt, das die Zahnflanken beschreibt. Die Anordnung der 
Koordinaten s und y, mit denen die dreidimensionale prismatische Oberfläche des Zahnrads in 
eine zweidimensionale Darstellung überführt werden kann, sind in Abbildung 4.6 dargestellt. Die 
Koordinate s verläuft entlang der Zahnkontur und hat ihren Ursprung an der Mitte des Zahnfußes. 
Die geometrische Symmetrieebene zwischen zwei Zähnen, die auch den periodischen Rand 
darstellt, verläuft durch s = 0. Das Ende des Koordinatenbereiches liegt in der Mitte des 
gegenüberliegenden Zahnfußes bei s = Smax. Damit kann die dimensionslose Koordinate 
ee (4.18) 


Smax 


eingeführt werden. Die axiale Koordinate y verläuft von der vorderseitigen Stirnwand bei y = 0 
bis zur rückseitigen Stirnwand bei y = W. Analog zu s* gilt für die dimensionslose axiale 


Koordinate 
(4.19) 


Abbildung 4.4: Eindimensionale Betrachtung des Einzelzahns 
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Abbildung 4.5: Verlauf der Temperaturen an den beiden Oberflächen des Ersatzkörpers, 0 und 
Ox, über einen Bereich von h für verschiedene Prüfkörperdicken L in m 


Abbildung 4.6: Flächentreues, abgewickeltes (s, y)- Koordinatensystem eines Einzelzahns 


In Abbildung 4.6 ist die Transformation in eine zweidimensionale Darstellung anschaulich 
dargestellt. Sie kommt einer Abwicklung der Zahnflanken auf eine Ebene gleich. Der Vorteil 
dieser Transformation gegenüber der Projektion liegt darin, dass sie flächentreu ist. Die ex- 
perimentell bestimmten Wärmeübergänge sind ortsaufgelöst. Die Bestimmung des mittleren 
Wärmeübergangs kann im (s, y)-Koordinatensystem durch einfache arithmetische Mittelung 
erfolgen, während dies für projizierte Koordinatensysteme eine Gewichtungsfunktion erfordern 
würde. 


Aus Gründen, die in Abschnitt 4.2 geschildert werden, sind die Temperaturmesspunkte des 
Versuchszahnrads über einen Viertelzahn verteilt. Dabei ist eine Hälfte dieses Viertelzahns mit 
18 Messpunkten versehen, die in drei Reihen entlang der Koordinate s* verteilt sind. Entlang 
der Mitte der anderen Hälfte ist eine Referenzreihe mit Messpunkten versehen. Die Verteilung 
der Messpunkte ist in Abbildung 4.7 dargestellt. 
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Abbildung 4.7: Verteilung der Messpunkte über einen Viertelzahn im (s*,y*)- 


Koordinatensystem 


Aus Abbildung 4.7 ist ersichtlich, dass der Bereich zwischen y* = 0,375 und yt = 0,5 mit drei 
Reihen mit jeweils sechs Messpunkten instrumentiert ist. Dieses Zahnachtel wird durch die 
iterative Berechnung ausgewertet. Bei y* = 0,3125 liegt eine weitere Reihe Messpunkte, die 
die Referenzmessreihe darstellt. Die dicht instrumentierte Hälfte des Viertels ist näher an der 
Zahnradmitte, um den Einfluss der Stirnseiten möglichst gering zu halten. Die Referenzmessreihe 
befindet sich mittig auf dem weiter außerhalb der Mitte befindlichen Achtel. Aufgrund der 
Beschränkung auf ein Zahnviertel wird in allen folgenden Abschnitten ausschließlich das Intervall 
yt = [0,25;0,5] betrachtet. 


Die Messpunkte des Versuchszahnrads werden durch Mantelthermoelemente vom Typ K mit 
einem nominellen Manteldurchmesser von 0,25 mm realisiert. Zur Applizierung der Thermo- 
elemente ist das Versuchszahnrad mit Nuten versehen, die axial entlang der Zahnflanken vom 
Ort des jeweiligen Messpunkts bis zur Stirnseite des Zahnrads verlaufen. Die Thermoelemen- 
te werden mit einem Klebstoff angepasster Wärmeleitfähigkeit in die Nuten eingeklebt. Die 
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Wärmeleifähigkeit des Klebstoffs ist mit 5,7W/(mK) nur geringfügig niedriger als die des 
Zahnradmaterials. Dadurch und durch die sehr geringe Größe der Nuten und der Thermoelemen- 
te kann die Beeinflussung der Temperaturverteilung gering gehalten werden. Zum Einkleben 
eines Thermoelements wird die Nut zunächst vollständig mit Klebstoff gefüllt. Anschließend 
wird das Thermoelement eingeklebt. Durch Auffüllen der Nut, Aushärten des Klebstoffes und 
anschließendes, manuelles Abschleifen stellen sich nur sehr geringe Geometriebeeinflussungen 
ein. 


Abbildung 4.8: In Nut eingeklebtes Mantelthermoelement 


Aufgrund der Anzahl an Messpunkten und konstruktiven Einschränkungen durch die Lage 
der Anschlusspunkte der Telemetrie können die Messpunkte nicht auf einem einzelnen Zahn 
angeordnet werden. Daher sind die Messpunkte über 24 Zähne verteilt. Das Temperaturfeld 
eines Zahns wird so nur durch eine Nut mit Thermoelement beeinflusst. Die Verteilung der 
Messpunkte ist in Tabelle 4.1 dargestellt. 
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Tabelle 4.1: Positionen der Messpunkte und Verteilung über das Zahnrad 


Messpunkt st yt Zahn 


1 0,0000 0,4792 60 
2 0,0000 0,4375 6l 
3 0,0000 0,3958 62 
5 0,1811 0,4792 65 
6 0,1811 0,4375 1 
7 0,1811 0,3958 2 
9 0,3621 0,4792 4 
5 
6 


10 0,3621 0,4375 
11 0,3621 0,3958 
13 0,5000 0,4792 29 
14 0,5000 0,4375 30 
15 0,5000 0,3958 31 
17 0,6605 0,4792 33 
18 0,6605 0,4375 34 
19 0,6605 0,3958 35 
21 0,8449 0,4792 37 
22 0,8449 0,4375 38 
23 0,8449 0,3958 39 
0,0000 0,3125 63 
8 0,1811 0,3125 3 
12 0,3621 0,3125 28 
16 0,5000 0,3125 32 
20 0,6605 0,3125 36 


24 0,8449 0,3125 40 


4.5 Analyse der Messgenauigkeit 


Die experimentelle Umsetzung der im Vorangegangenen beschriebenen Methode unterliegt 
mehreren Störeinflüssen, die die Messgenauigkeit beeinflussen. Dazu gehören unter anderem 
Unsicherheiten in der Bestimmung der Wärmeleitfähigkeit des Materials, Messfehler in der 
Temperaturerfassung durch Thermoelemente und geometrische Abweichungen der Zahnradform 
und der Position der Messpunkte. In Arbeiten wie der von Jacobsen (1987) oder Stripf (2007) 
wird die aus den verschiedenen Störquellen resultierende Unsicherheit der Wärmeübergangsmes- 
sung in zwei Schritten quantifiziert. Dazu werden zunächst die Unsicherheiten der beteiligten 
Größen, wie etwa Temperatur oder Wärmeleitfähigkeit, durch Herstellerangaben oder Kalibrier- 
messungen festgestellt. Anschließend wird der Einfluss dieser Einzelfehler auf den gemessenen 
Wärmeübergang bestimmt, indem die Wärmeleitungsgleichung entsprechend ausgewertet wird. 


Eine derartige analytische Betrachtung der Fehlerfortpflanzung wird im vorliegenden Fall da- 
durch verhindert, dass die Geometrie des Versuchskörpers komplex ist und die Messpunkte dicht 
beieinander liegen. Dadurch beeinflussen Messfehler an einem bestimmten Messpunkt auch die 
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ermittelten Wärmeübergangskoeffizienten an benachbarten Stellen. Dies ist ein Nebeneffekt der 
Auslegung des Versuchskörpers unter Annäherung von Bi = 1. Eine grobe Näherung der Unsi- 
cherheit in der Wärmeübergangsbestimmung könnte erreicht werden, indem das Versuchszahnrad 
modellhaft zu einem eindimensionalen Vergleichskörper vereinfacht wird. Dadurch würde je- 
doch die zu erwartende Messunsicherheit stark überschätzt. Die Wechselwirkung zwischen 
Messpunkten im vorliegenden Fall führt zu einem stochastischen Ausgleich der aufgeprägten 
Messunsicherheiten, wodurch die mittlere Messunsicherheit geringer ausfällt als für eine einzelne 
Temperaturmessung. 


Die Messunsicherheit der beschriebenen Methode kann dennoch bestimmt werden, wenn sie 
als statistisches Problem betrachtet und mit der Monte-Carlo-Methode quantifiziert wird. Die 
Monte-Carlo-Methode soll im Folgenden beschrieben und die so gewonnenen Erkenntnisse 
über die Messgenauigkeit der vorliegenden Wärmeübergangsmessungen anschließend analysiert 
werden. 


4.5.1 Monte-Carlo-Methode 


Die analytische Betrachtung der Fehlerfortpflanzung setzt zwingend die geschlossene Beschreib- 
barkeit des jeweiligen Vorgangs voraus. In Fällen, in denen die Natur eines Vorgangs entweder 
unbekannt aber beobachtbar ist (z.B. Finanzmärkte) oder numerisch berechenbar, aber nicht 
geschlossen lösbar ist (wie im vorliegenden Fall), ist diese Voraussetzung nicht gegeben. Die 
Abhängigkeit solcher Vorgänge von ihren Eingangsparametern kann statistisch bestimmt werden. 
Bei bekannten Unsicherheiten der Eingangsgrößen ist dadurch die Unsicherheit der Ausgangs- 
größe bekannt - im Fall der Wärmeübergangsmessung also die Unsicherheit des Wärmeüber- 
gangskoeffizienten. 


In Abbildung 4.9 ist die Vorgehensweise der Monte-Carlo-Methode zur Analyse der Fehler- 
fortpflanzung schematisch dargestellt. Der zu Grunde liegende Gedanke ist, dass stochastische 
Störungen der Eingangsgrößen eines Vorgangs sich auf eine Art in den Ausgangsgrößen nie- 
derschlagen, die nur von diesem Vorgang selbst abhängen. Um die statistische Verteilung der 
Ausgangsgröße zu bestimmen, wird eine große Anzahl sogenannter Ziehungen getätigt. Eine 
Ziehung ist eine Diskretisierung aller als Wahrscheinlichkeitsgrößen beschriebenen Eingangsgrö- 
ßen. Jede Ziehung wird durch das Modell auf die Ausgangsgröße abgebildet. Die so entstandene 
Ausgangsgrößenverteilung kann statistisch analysiert werden, um die Messunsicherheit der 
Methode zu bestimmen. 


Dazu werden die Eingangsgrößen des Vorgangs als stochastische Größen erfasst. Im vorliegenden 
Fall bedeutet dies die Quantifizierung der Messunsicherheit aller Temperaturmessungen sowie 
der Unsicherheit der Warmeleitfahigkeit des Zahnradmaterials. Die numerische Berechnung 
der Ausgangsgrößenverteilung des Messvorgangs wurde in Abschnitt 4.3 beschrieben. Sie 
kann im Rahmen der Monte-Carlo-Methode als Zuordnungsvorschrift einer Einzelziehung von 
Eingangsgrößen zu einer Auswertung verwendet werden. 


Aufgrund der so gewonnenen Daten kann die Unsicherheit des mittleren Wärmeübergangs an der 
Zahnoberfläche bestimmt werden. Dadurch, dass der Wärmeübergang nicht nur als Mittelwert, 
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Abbildung 4.9: Fehlerfortpflanzungsanalyse durch die Monte-Carlo-Methode nach Wiibbeler 
et al. (2008) 


sondern örtlich aufgelöst bestimmt wird, kann die Messunsicherheit an jeder einzelnen Messstelle 
bestimmt werden. Durch die Verfiigbarkeit einer hohen Anzahl an Auswertungen zufalliger 
Messwertvariationen können zusätzlich die Abhängigkeiten der Ausgangsgrößen von den Ein- 
gangsgrößen durch ihre jeweiligen Korrelationskoeffizienten ausgedrückt werden. Dadurch kann 
die Robustheit der Methode geprüft werden, indem die beobachteten Abhängigkeiten analysiert 
werden. Im vorliegenden Fall sind dies die Abhängigkeiten des Wärmeübergangskoeffizienten 
an einem Messpunkt von den gemessenen Temperaturen an diesem und den ihn umgebenden 
Messpunkten. 


Zur Durchführung der Monte-Carlo-Simulation muss eine Wärmeübergangsverteilung angenom- 
men werden, um eine Temperaturverteilung zu berechnen, aus der lokale Messtemperaturen ex- 
trahiert werden können. Diese berechneten Temperaturen werden anschließend für jede Ziehung 
mit zufälligen Messfehlern überlagert und ausgewertet. Die angenommene Wärmeübergangsver- 
teilung ist in Abbildung 4.10a dargestellt. Die Verteilung des Wärmeübergangskoeffizienten auf 
der strömungszugewandten Zahnflanke entspricht entlang der y-Koordinate einer Sinuskurve, 
deren Maximum von 1000 W/(m?K) in der Mitte des Viertelzahns liegt und zu den Rändern 
hin auf einen Wert von Null abfällt. Auf dem Zahnkopf bei s* = 0,5 beträgt der Wärmeüber- 
gangskoeffizient 200 W/(m? K), auf der strömungsabgewandten Zahnflanke 100 W/(m?K). Die 
daraus resultierende, ungestörte Temperaturverteilung ist in Abbildung 4.10b dargestellt. 


Die ungestörten Temperaturen, die sich aus der angenommenen Wärmeübergangsverteilung 
ergeben, werden für jede Ziehung mit zufälligen Messfehlern überlagert. Dabei sind für die 
Temperaturmessung des Öls und die Temperaturmessungen auf dem Zahnrad unterschiedliche 
Messunsicherheiten zu erwarten. 


Die Messunsicherheit der Öltemperatur, ôW, setzt sich zusammen aus der Messunsicherheit 
der Kaltstellentemperatur, 60x, und der Messunsicherheit V, r, die durch herstellungsbedingte 
Fehler des Thermoelements und die Signalwandlung entstehen. Die Häufigkeitsverteilung der zu 
erwartenden Messunsicherheit wird aufgrund von Herstellerangaben als Normalverteilung mit 
einer Standardabweichung von o = 0,64K abgeschätzt. 


Die Signalwandlungskette der Telemetrie, mit Hilfe derer die Temperaturen am Zahnrad ge- 
messen werden, ist wesentlich komplexer als die des Signalwandlers der Öltemperatur. Durch 
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(a) Angenommener Wärmeübergang in W/(m? K) (b) Berechnete Temperaturverteilung in K 


Abbildung 4.10: Angenommener Wärmeübergang und berechnete Temperaturverteilung für die 
Monte-Carlo-Simulation 


Kalibrierdaten des Herstellers kann jedoch die Unsicherheit Urg bestimmt werden. Die Unsi- 
cherheit der Messung der Temperaturen an der Zahnoberfläche wird durch eine Wahrscheinlich- 
keitsdichtefunktion, die aus Kalibrierdaten gewonnen wird, abgebildet. Die Standardabweichung 
der Messung beträgt o = 0,21 K. 


Aufgrund des Rechenaufwands, der pro Ziehung bei etwa 15 Minuten liegt, und der Vielzahl 
an notwendigen Ziehungen wird im vorliegenden Fall das Stein’sche Schema zur Reduktion 
des Rechenaufwands verwendet (Wiibbeler et al. (2010)). Bei Anwendung dieses Verfahrens 
werden zunächst einige Blöcke, bestehend aus etwa 10! bis 10° Ziehungen, berechnet. Aufgrund 
der statistischen Verteilung der Blockmittelwerte kann anschließend eine Aussage über die 
benötigte Gesamtzahl an Ziehungen getroffen werden, die zum Erreichen eines definierten 
Konfidenzintervals der Schätzung der Messunsicherheit ausgewertet werden müssen (Wübbeler 
et al. (2010)). Im vorliegenden Fall ergibt die Auswertung von 8 Blöcken mit jeweils 50 
Ziehungen eine benötigte Gesamtzahl von 5650 Ziehungen, um die Unsicherheit der Methode zu 
charakterisieren. Dabei wird eine Charakterisierung der Unsicherheit mit einer Abweichung von 
1% und einem Konfidenzintervall von 95% angestrebt. 


4.5.2 Unsicherheit der Wärmeübergangsmessungen 


Auf Grundlage der zuvor genannten Unsicherheiten der Eingangsgrößen kann durch die Monte- 
Carlo-Methode die Unsicherheit der Wärmeübergangsmessungen bestimmt werden. Dabei 
können sowohl die Unsicherheit des mittleren Wärmeübergangs als auch die lokalen Messunsi- 
cherheiten an den Messpunkten betrachtet werden. 


Zunächst soll die Unsicherheit des gemittelten Wärmeübergangs betrachtet werden. In der Arbeit 
von Niehues (2018) wird der verwendete zweistufige Stein’sche Ansatz zur Bestimmung der 
Messunsicherheitsverteilung beschrieben, dessen Ergebnisse im Folgenden ausgewertet werden 
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sollen. 


Die in das Zahnrad übertragenen lokalen Wärmeströme, die mit der iterativen Methode berechnet 
werden, bilden zusammen mit der Öltemperatur und der lokalen Oberflächentemperatur den 
Wärmeübergangskoeffizienten h. Die Unsicherheit der Messung der Öltemperatur hat dadurch 
einen Einfluss auf die Unsicherheit des mittleren Wärmeübergangs. Die Häufigkeitsverteilung 
des mittleren Wärmeübergangskoeffizienten h aus allen Ziehungen der Monte-Carlo-Simulation 
istin Abbildung 4.11a dargestellt. Die Verteilung entspricht einer Normalverteilung mit dem 
Mittelwert Am = 298,80 W/(m? K) und einer Standardabweichung von o = 11,24W/(m?K) 
oder o = 3,76% des Mittelwerts. Aus Abbildung 4.11b ist ersichtlich, dass die Öltemperatur ¥, 
stark mit dem mittleren Wärmeübergangskoeffizienten h korreliert. Der Korrelationskoeffizient 
beträgt r= —0,9815. 
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Abbildung 4.11: Messunsicherheit des mittleren Wärmeübergangs 


Aufgrund der Geometrie des Zahnrads und der Auslegung zu Bi ~ 1 sind die Unsicherheiten 
der lokalen Wärmeübergangsmessung von den Messunsicherheiten aller Eingangsparameter 
abhängig. Jedoch kann gezeigt werden, dass die Öltemperatur und die lokale Temperatur einen 
dominierenden Einfluss auf die Messunsicherheit haben und die benachbarten Messpunkte sie 
nur nachrangig beeinflussen. 


Die Korrelationskoeffizienten der berechneten lokalen Wärmeübergangskoeffizienten h; und 
der Öltemperatur sind in Abbildung 4.12 dargestellt. Die Öltemperatur korreliert dabei stets 
negativ mit den Wärmeübergangskoeffizienten. Messpositionen hohen Wärmeübergangs weisen 
eine stärkere Abhängigkeit von der Öltemperatur auf, da die Temperaturdifferenz zwischen 
Oberfläche und Öltemperatur an diesen Positionen niedriger ist. Die Temperaturdifferenz, die in 
den Nenner des Wärmeübergangskoeffizienten eingeht, ist somit kleiner (vgl. Gleichung 4.6). 
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Abbildung 4.12: Abhängigkeit des lokalen Wärmeübergangskoeffizienten von der Messunsicher- 
heit der Öltemperatur 


Die Streuung der Wärmeübergangskoeffizienten an Messpunkt 9 ist in Abbildung 4.13a darge- 
stellt. Der betrachtete Messpunkt liegt in der Nähe der Zahnmitte und nah am Zahnkopf auf der 
strömungszugewandten Zahnflanke. Die Verteilung der gemessenen Wärmeübergangskoeffizien- 
ten entspricht einer Normalverteilung mit dem Mittelwert ho m = 1057,3 W/(m? K) und einer 
Standardabweichung von o = 112,9W/(m?K). Aus Abbildung 4.13b wird deutlich, dass der 
ermittelte Wärmeübergangskoeffizient stark von der lokalen Temperatur abhängt. Mit einem 
Korrelationskoeffizienten von r = 0,8189 hat die lokale, gemessene Temperatur den größten 
Einfluss auf die Messunsicherheit des lokalen Wärmeübergangskoeffizienten. 
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Abbildung 4.13: Unsicherheit des lokalen Wärmeübergangs an der Messposition 9 


Durch die Darstellung der Korrelationskoeffizienten des gemessenen Wärmeübergangskoef- 
fizienten an Position 9 und aller gemessenen Temperaturen in Abbildung 4.14 wird deutlich, 
dass auch die Temperaturen der nächstgelegenen Temperaturmesspunkte 6, 8 und 12 den an 
Position 9 gemessenen Wärmeübergangskoeffizienten beeinflussen. Dies ist eine Auswirkung 
des beabsichtigten örtlichen „Mittelungseffektes“ bei Bi ~ 1. Die weiter entfernt gelegenen 
Messpunkte beeinflussen den lokalen Wärmeübergangskoeffizienten kaum. 
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Abbildung 4.14: Korrelationskoeffizienten des Wärmeübergangs an Position 9 und der Tempera- 
turen an den Messpunkten 


Die beobachteten Einflüsse auf die Unsicherheit der lokalen Wärmeübergangsmessung an Positi- 
on 9 werden auch an allen anderen Messpunkten beobachtet. Die näher am Zahnfuß gelegenen 
Messpunkte zeigen dabei tendenziell weiträumigere Abhängigkeit von anderen Messstellen. Die 
ausführliche Diskussion der jeweiligen Unsicherheiten und Korrelationskoeffizienten soll an 
dieser Stelle jedoch nicht erfolgen. In Abbildung 4.15 ist |r| für alle Kombinationen von ®,; 
und h; dargestellt. Werte der Hauptdiagonale stellen darin die Abhängigkeit eines berechneten, 
lokalen Wärmeübergangskoeffizienten h; von der gemessenen Temperatur v; an der Position 
i dar. Durch diese Form der Darstellung wird deutlich, dass die lokale Temperatur stets den 
größten Einfluss auf den berechneten Wärmeübergangskoeffizienten hat. 
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Abbildung 4.15: Korrelationskoeffizienten aller lokalen Wärmeübergangskoeffizienten von allen 
Oberflächentemperauren 
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4.5.3 Zusammenfassung 


Die vorangegangene Betrachtung liefert die Erkenntnis, dass die Messunsicherheit des mittleren 
Wärmeübergangs vergleichsweise gering ist. Für den verwendeten Testfall wird eine Standard- 
abweichung von © = 3,76% ermittelt. Die Verteilung und der Wertebereich des Testfalls sind 
den gemessenen Wärmeübergängen dieser Arbeit ähnlich. Die Unsicherheit der Öltempera- 
turmessung hat den größten Einfluss auf die Messunsicherheit des Wärmeübergangs. Durch 
Kalibrierung der verwendeten Temperaturmesskette für die Öltemperatur wird eine wesentlich 
geringere Unsicherheit der Öltemperatur festgestellt. Es ist daher davon auszugehen, dass die 
tatsächliche Unsicherheit des mittleren Wärmeübergangs ebenfalls wesentlich geringer als die 
durch die Monte-Carlo-Simulation bestimmte Unsicherheit ist. Die Unsicherheiten der lokalen 
Wärmeübergangsmessungen übersteigen die des mittleren Wärmeübergangs erheblich. Sie 
werden am stärksten von den lokalen Temperaturmessungen und der Messung der Öltemperatur 
beeinflusst. Die Unsicherheiten der jeweils benachbarten Messpunkte haben ebenfalls einen 
eindeutigen, jedoch geringeren Einfluss. 
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Experimentelle Methode zur Bestimmung des Öl-Zahnrad-Wärmeübergangs 


5 Versuchsaufbau und Messtechnik 


Aus Kapitel 2 geht hervor, dass die Charakterisierung des Öl-Zahnrad-Wärmeübergangs für 
das Erreichen höchster Leistungsdichten in Getrieben unabdingbar ist. Die im Rahmen der 
vorliegenden Arbeit untersuchte Ausgestaltung einer Einspritzkühlung wird in Kapitel 3 be- 
schrieben. In Kapitel 4 wird die experimentelle Methode beschrieben, die die Messung des 
Wärmeübergangs zwischen Ol und Zahnrad ermöglicht. Die zur Erfüllung der Zielsetzung dieser 
Arbeit notwendige Implementierung dieser Methode in Form eines Versuchsstands und einer 
Versorgungsanlage soll im Folgenden beschrieben werden. 


5.1 Implementierung in einem Versuchsstand 


Zur Umsetzung der in Kapitel 4 beschriebenen Methode sind folgende Anforderungen an Ver- 
suchsstand und Versorgungsanlagen zu erfüllen: 


e Einspritzen des beheizten Ols 

e Variation des Düsenwinkels 

e Auffangen und Abführen des Ols 

e Kühlen des Versuchszahnrads 

e Antrieb des Versuchszahnrads 

e Ermöglichen von Gegenlichtaufnahmen 


e Messen von Temperaturen im rotierenden System 


Aus der Beschreibung der Methode in Kapitel 4 wird deutlich, weshalb ein Temperaturgradient 
zur Messung des Wärmeübergangs an einem einzelnen Zahnrad erzeugt werden muss. Dies 
geschieht im vorgestellten Versuchsstand durch Heizen des Öls und Kühlen des Zahnradinnenra- 
dius, was eine Umkehr der Wärmestromrichtung im Vergleich zum Anwendungsfall darstellt. 
So können die Luftversorgung, in der keine Heizung benötigt wird, und die Prüfstandskonstruk- 
tion, die geringeren Temperaturen standhalten muss, einfacher gestaltet werden. Es werden 
Betriebsbedingungen angestrebt, die denen existierender Höchstleistungsgetriebe in Triebwerken 
entsprechen. Die Betriebsgrenzen des Versuchsstands sind in Tabelle 5.1 zusammengefasst. 


Um die genannten Anforderungen zu erfüllen, wurde der in Abbildung 5.1 dargestellte Versuchs- 
stand entworfen und am Institut für Thermische Strömungsmaschinen gefertigt. 


Das Versuchszahnrad, das in Kapitel 4 beschrieben wird, ist durch eine Wälzlagerung gelagert 
und wird vom Antriebsmotor über eine Drehmomentmesswelle direkt angetrieben. Durch die 
„fliegende“ Lagerung der Drehmomentmesswelle führen Ausrichtungsfehler, die nicht vollständig 
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Tabelle 5.1: Betriebsgrenzen des Versuchsstands 


Parameter Symbol Wert 


Drehzahl n 6500/min 
Olvolumenstrom V 151/min 
Düsenwinkel a —45° bis 45° 
Oltemperatur Us 100°C 
Olbrause 
| Verstellantrieb 
| Motor 
Drehmomentmesswelle 


Lagerung 


Düsenkasten 


Telemetrieantenne Zahnrad 


Telemetriesender 


Abbildung 5.1: Schematische Darstellung des eingesetzten Versuchsstands zur Messung des 
Öl-Zahnrad-Wärmeübergangs 


zu vermeiden sind, nicht zu umlaufenden Biegemomenten im Antriebsstrang, was zu einem 
äußerst schwingungsarmen Betrieb führt. 


Eine stirnseitige Labyrinthdichtung dichtet die ölhaltige Umgebung, in der sich die verstellbare 
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Olbrause befindet, gegen die Umgebung des Prüfstands ab. Innerhalb des hohlen Zahnrads befin- 
det sich ein zylinderförmiger Düsenkasten, der zur Prallkühlung der Innenfläche des Zahnrads 
durch einen Seitenkanalverdichter mit Luft versorgt wird. Die Vorderseite des Versuchsstands ist 
mit einem Fensterrahmen versehen, der als integrales Bauteil sowohl den Düsenkasten als auch 
die Telemetrieantenne aufnimmt. Der Telemetriesender ist in das Zahnrad integriert. 


Die Ölbrause ist oberhalb des Zahnrads angebracht. Sie ist schwenkbar gelagert und kann, 
ausgehend von der Neutrallage, durch einen Verstellantrieb während des Versuchs um 45° in 
beide Richtungen geschwenkt werden. Die Neutrallage ist dabei die Lage, in der der Ölstrahl in 
Richtung der Zahnradachse gerichtet ist. Die Ölbrause besitzt vier Auslässe, die in Abständen 
von jeweils 0,25 W gleichmäßig über die Tiefe des Zahnrads W = 90mm angeordnet sind 
(siehe Abbildung 5.2) und einen Auslassdurchmesser von d = 1,15 mm bei einer Bohrungslänge 
von / = 7mm besitzen. Die Düsenbohrungen werden durch Hartmetallbohrer gefertigt, erst 
anschließend wird die Auslassfläche gefräst. Dadurch kann sichergestellt werden, dass der 
Auslass gratfrei ist. Die Ölzuführung zu den Düsen erfolgt über eine Kernbohrung mit einem 
Durchmesser von 10mm, die stirnseitig in die Olbrause eingebracht und mit einem hartgelöteten 
Deckel verschlossen wird. Durch boroskopische Inspektion wird die Freiheit von Riefen und 
Graten innerhalb der Düse sichergestellt. Die Ölbrause wird durch eine flexible Leitung, die 
die Schwenkbewegung ermöglicht, und eine rückseitige Kernbohrung mit Öl versorgt. Die 
vorderseitige und die rückseitige Kernbohrung weisen einen Achsabstand von 2,1 mm auf. 
Durch numerische Untersuchungen kann gezeigt werden, dass durch den Achsabstand und die 
dadurch entstehende Stufe nur sehr geringe Störungen in der Anströmung der Düsen entstehen 
(Kanzow (2017)). Durch ein mittig in den Einlass der Ölbrause eingelassenes Thermoelement 
kann die Öltemperatur unmittelbar vor den Düsen gemessen werden. 


0,125W 0,25W 025W 025W 2,1mm 


- 


Messposition % 


1W 


+ > 


Abbildung 5.2: Schnittdarstellung der Ölbrause mit Zuführbohrung durch die Düsenebene 


Durch ein vorderseitiges und ein rückseitiges Fenster und unter Verwendung einer Gegenlichtquel- 
le können Hochgeschwindigkeits-Schattenbildaufnahmen des Einspritzvorgangs aufgenommen 
werden. Die Ölbrause ist so gestaltet, dass sich der Austritt aus den Öldüsen auf dem größten 
Durchmesser der Ölbrause befindet. Der Ölstrahl, der die Ölbrause verlässt, ist somit in den 
Aufnahmen auf seiner gesamten Länge sichtbar. Durch die Verwendung eines geradverzahnten 
Zahnrads ist das Auftreffen des Ölstrahls auf die Zahnradoberfläche beobachtbar. Die verwendete 
Hochgeschwindigkeitskamera vom Typ LaVision HighSpeedStar 8 ermöglicht Aufnahmese- 
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quenzen mit einer Bildrate von 27 kHz bei Bildausschnitten mit 512x512 Pixeln. Dadurch sind 
Aufnahmen möglich, mit denen die Ausbreitung und das Auftreffen des Olstrahls mit mehreren 
Einzelbildern aufgenommen werden können. Bei n = 3000/min sind so beispielsweise acht 
Einzelbilder eines Auftreffvorgangs erfassbar. 


5.2 Anlagen zur Ölversorgung und Luftkühlung 


Die Bereitstellung des benötigten heißen Öls erfolgt über eine Versorgungsanlage, die in Ab- 
bildung 5.3 schematisch dargestellt ist. Durch den Heizkreislauf, bestehend aus einer Zahn- 
radpumpe und einer elektrischen Heizung, wird das Öl aus dem Tank gefördert, erwärmt und 
in den Tank zurückgeführt. Über einen PID-Prozessregler wird die Temperatur am Heizungs- 
auslass auf einen Sollwert geregelt. Die Versorgung des Versuchsstands erfolgt über den Ver- 
sorgungskreislauf, bestehend aus einer Zahnradpumpe, einem Filter und einem hochgenauen 
Coriolis-Massendurchflusssensor. Die Förderleistung der Zahnradpumpe wird durch einen 
Frequenzumrichter und einen PID-Prozessregler geregelt, um einen voreingestellten Massen- 
durchfluss zu erreichen. Zur Sicherstellung von Blasenfreiheit wird die Dichte des Öls, die 
durch den Massendurchflusssensor gemessen wird, überwacht. Die ölführenden Leitungen der 
Versorgungsanlage sind nahezu vollständig mit Isoliermaterial ummantelt, um Temperatur- und 
Dichteunterschiede sowie thermische Verluste gering zu halten. Eine weitere Zahnradpumpe 
saugt das Öl, das den Boden des Versuchsstands erreicht, durch drei Ablaufbohrungen ab. 


Heizkreislauf Versorgungskreislauf | Versuchsstand 


EEE 


Regelkreis 
Volumenstrom 


„Regelkreis 

Oltemperatur Olbrause 
Zahnrad 
Diisenkasten 


is nah tein is aia es | 


Plenum 


Tank | Rückführpumpe Seitenkanalverdichter 
Abbildung 5.3: Anlagenplan des eingesetzten Versuchsstands zur Messung des Öl-Zahnrad- 
Wärmeübergangs 


Um die notwendige Kühlung des Zahnrads auf seiner Innenfläche, die in Kapitel 4 erläutert 
wird, zu gewährleisten, wird der Zahnradinnenmantel durch einen Düsenkasten mit gleichmäßig 
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angeordneten Auslassbohrungen gekühlt. Der Diisenkasten und das Zahnrad sind konzentrisch 
angeordnet. Jeweils fiinf Auslassbohrungen sind in einer axialen Reihe angeordnet und werden 
durch eine gemeinsame Bohrung mit Luft versorgt. Der Diisenkasten besitzt 40 dieser Reihen, 
wodurch insgesamt 200 Kiihlluftbohrungen mit einem Durchmesser von jeweils 2,5 mm zur 
Kühlung eingesetzt werden. Der radiale Abstand der Auslassbohrungen zum Zahnrad beträgt 
15 mm. Über ein Plenum im Düsenkasten und 20 Schlauchanschlussstücke ist der Düsenkasten 
über einen Luftsammler mit dem Seitenkanalverdichter verbunden. 


5.3 Messtechnik 


Die Messung der Temperaturen am Zahnrad und die Berechnung des Wärmeübergangs sowie die 
dabei herrschende Messunsicherheit werden in Kapitel 4 beschrieben. Die restlichen Betriebs- 
parameter unterliegen ebenfalls messtechnischen Unsicherheiten, die im Folgenden beleuchtet 
werden sollen. Die Genauigkeiten der jeweiligen Messinstrumente sind in Tabelle 5.2 aufgeführt. 


Ölvolumenstrom 


Die Messung des Ölvolumenstroms erfolgt durch einen Durchflussmesser vom Typ Krohne 
Optimass 6400C. Der Sensor arbeitet nach dem Coriolis-Prinzip, durch das bei einphasigen 
Fluiden der Massenstrom und die Dichte direkt gemessen werden. Der Volumenstrom wird aus 
den beiden Größen berechnet. Dadurch, dass die Dichte des Fluids stets gemessen wird, kann 
durch einen Abgleich mit den Stoffdaten des verwendeten Öls sichergestellt werden, dass keine 
Luftblasen im Öl vorhanden sind. Diese können aufgrund des Aufbaus der Versorgungsanlage 
als Kreislauf und durch die Mischung von Öl und Luft prinzipiell auftreten. In den Experimenten, 
die in der vorliegenden Arbeit beschrieben werden, kann jedoch zu keinem Zeitpunkt eine 
Dichteänderung des Öls aufgrund von Luftblasen festgestellt werden. 


Öltemperatur 


Die Temperatur des Öls wird sowohl im Heizkreislauf als auch in der Zuführbohrung der Öldüse 
gemessen. Dabei kommen Mantelthermoelemente vom Typ K mit einem Manteldurchmesser von 
0,5 mm zum Einsatz. Diese sind über Ausgleichsleitungen mit dem Prozessregler bzw. einem 
Signalerfassungsmodul vom Typ National Instruments NI-9211 verbunden. Die Kaltstellenkom- 
pensation erfolgt über integrierte temperaturabhängige Widerstände. 


Drehmoment und Drehzahl 


Die Drehzahl des Zahnrads und das Antriebsdrehmoment werden durch eine Drehmomentmess- 
welle vom Typ Lorenz DR-3000 gemessen. Zur Schwingungsreduktion ist der Sensor zwischen 
zwei Lamellenkupplungen angebracht und nicht fest mit der Bodenplatte des Versuchsstands 
verbunden. 


Düsenwinkel 


Die Ölbrause ist schwenkbar gelagert und wird durch einen Servoantrieb verstellt. Über ein 
Potentiometer, das koaxial zur Antriebswelle des Servos und in den Servo integriert angebracht 
ist, kann der Drehwinkel der Düse ausgelesen werden. Das Spannungssignal des Drehwinkels 
wird über ein Signalerfassungsmodul vom Typ NI-9201 erfasst. 
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Tabelle 5.2: Genauigkeit 


Parameter Gerät Bereich Genauigkeit 
Ölvolumenstrom Krohne Optimass 6400C 0 —201/min +0,05 % 
Drehzahl Lorenz DR-3000 0 —20000/min +1% 
Öltemperatur (Düse) Thermoelement Typ K —270—1300°C ø =0,5°C 
Düsenwinkel Servoantrieb / Potentiometer —45 —45° +2% 


Hochgeschwindigkeitsaufnahmen 


Um die Interaktion des Ölstrahls mit dem Zahnrad zu charakterisieren und das Verhalten des 
Wärmeübergangs einordnen zu können, ist der Versuchsstand mit vorder- und rückseitigen Fens- 
tern ausgestattet. Das Zahnrad und die Ölbrause sowie der Antrieb der Ölbrause sind so gestaltet, 
dass keine optischen Versperrungen zwischen der Austrittsfläche des Öls und der Zahnkontur vor- 
handen sind. Dadurch sind Hochgeschwindigkeitsaufnahmen im Gegenlichtverfahren möglich. 
Zur Aufzeichnung dieser Gegenlichtaufnahmen werden eine LED-Hintergrundbeleuchtung vom 
Typ Phlox HSC und eine Hochgeschwindigkeitskamera vom Typ LaVision HighSpeedStar 8 


eingesetzt. 
Temperaturen im rotierenden System 


Zur Erfassung der Zahnradtemperaturen auf den Zahnflanken und dem Innenradius des Zahn- 
rads ist die kabellose Übertragung der Signale notwendig. Das verwendete Telemetriesystem 
ermöglicht die Erfassung von 28 Temperaturkanälen. Dabei wird die Kaltstellentemperatur über 
einen der Kanäle erfasst. Die Erfassung der Temperaturen erfolgt dabei sequentiell mit einer 
Rate von 1,6Hz. Der Aufbau der Senderelektronik bewirkt eine Mittelung des Thermospan- 
nungssignals. Zur Reduktion des Rauschens wird das digitalisierte Signal empfängerseitig über 
einen Ringspeicher gemittelt. 


6 Ergebnisse 


Um die Zielsetzung dieser Arbeit zu erfiillen, wird im Folgenden der Ol-Zahnrad-Warmeiibergang 
experimentell untersucht. Durch Variation der Betriebsparameter kann in dieser Arbeit erst- 
malig das charakteristische Verhalten des Öl-Zahnrad-Wärmeübergangs in Abhängigkeit der 
Betriebsparameter beschrieben werden. 


Der Wärmeübergang an der Zahnradoberfläche wird nach der in Kapitel 4 beschriebenen Me- 
thode ermittelt. Zusätzlich werden Hochgeschwindigkeitsaufnahmen unter Verwendung einer 
Gegenlichtquelle aufgezeichnet, um Rückschlüsse auf die auftretende Strömung und die Öl- 
Zahnrad-Interaktion zu ziehen. 


Für die Bestimmung des Wärmeübergangs werden die Betriebsbedingungen durch die in Kapi- 
tel 5 beschriebenen Anlagen erzeugt. Variable Betriebsbedingungen im Rahmen dieser Untersu- 
chung sind: 


e Drehzahl n 
e Ölvolumenstrom V 
¢ Diisenwinkel & 


e Oltemperatur w. 


Die untersuchten Betriebspunkte sind in Abbildung 6.1 dargestellt. 
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Abbildung 6.1: Verteilung der untersuchten Betriebspunkte im Parameterraum 


Von insgesamt 129 Betriebspunkten entfallen 110 auf Diisenwinkel œ von 0°, 20° und 45°. Für 
diese Düsenwinkel werden die Drehzahl n in einem Bereich zwischen 1000/min und 4000/min 
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und der Olvolumenstrom V zwischen 3,21/min und 7,2 1/min systematisch variiert. Für œ = 20° 
ist dieser Bereich am feinsten aufgelöst. Die Diskussion des Einflusses von Drehzahl und 
Ölvolumenstrom erfolgt für & = 20° in Abschnitt 6.2.1. Um den Einfluss des Düsenwinkels 
zu charakterisieren, werden in Abschnitt 6.2.2 die zusätzlichen Messungen bei & = 45° sowie 
Messungen bei weiteren Düsenwinkeln zwischen —30° und 45° diskutiert. Eine vollständige 
Liste der Betriebspunkte und die ortsaufgelösten Darstellungen des Wärmeübergangs finden sich 
in Anhang A.1 dieser Arbeit. 


6.1 Messung und Auswertung des Wärmeübergangs 


Im Folgenden soll der Einfluss der variierten Parameter auf den Wärmeübergang beschrieben 
werden. Zunächst soll jedoch ein Betriebspunkt exemplarisch analysiert werden, um die Vorge- 
hensweise zu verdeutlichen. 


Alle gemessenen Größen werden während der Versuche aufgezeichnet, einige werden zur Rege- 
lung verwendet. Die Berechnung der Wärmeströme erfolgt auf Grundlage der so gewonnenen 
Daten. 


Die benötigte Zeit, um das Öl, die Versorgungsanlage und den Prüfstand aufzuheizen und 
in einen stationären Zustand zu bringen, ist lang im Vergleich zu der Zeit, die zum Wechsel 
zwischen zwei Betriebspunkten benötigt wird. Da aus der Literatur nur wenige Informationen zur 
Beschreibung des Öl-Zahnrad-Wärmeübergangs zur Verfügung stehen, wird ein großer Bereich 
des Parameterraums durch Messungen abgedeckt, um alle wesentlichen Einflüsse aufzudecken 
(siehe Abbildung 6.1). Um die Erfassung der dafür benötigten großen Anzahl von Messpunkten 
effizient zu gestalten, werden mehrere Messungen in ähnlichen Parameterbereichen unmittelbar 
hintereinander durchgeführt. In Abbildung 6.2 sind die Verläufe der gemessenen Temperaturen, 
der Drehzahl und des Ölvolumenstroms über die Zeit für einen typischen Versuchsverlauf 
dargestellt. Die dargestellten Temperaturen im oberen Diagramm sind die des Öls in der Düse 
sowie ausgewählter Messpositionen an der Zahnradkontur und am Innenradius des Zahnrads. 


Nach dem Beginn des Aufheizen des Öls bei t ~ 0,3 - 10°s ist ein Ansteigen aller Temperaturen 
zu erkennen. Bei t ~ 3-10°s erreicht die Öltemperatur den Sollwert von 80°C. Zu diesem 
Zeitpunkt sind die Temperaturen des Zahnrads noch nicht stationär. Aufgrund der graduellen 
Änderung der sich einstellenden Temperaturen wird der Betriebspunkt gehalten, bis keine Än- 
derungen der Temperaturen am Zahnrad und in der Öldüse mehr feststellbar sind, die über das 
Rauschen der Messwerte hinausgehen. Bei Erreichen des stationären Zustands wird dieser min- 
destens eine Minute gehalten, um eine ausreichende Datenbasis für die Mittelung zu erzeugen. 
Die letzten 30 Messwerte werden zur Mittelwertbildung herangezogen. Aufgrund der eingestell- 
ten Datenrate von 1 Hz entspricht dies einer Messung über einen Zeitraum von 30 Sekunden. 
Anschließend wird der nächste Parametersatz eingestellt und gehalten, bis erneut stationäre 
Zustände beobachtet werden. In Abbildung 6.2 ist dies durch den sequentiellen Wechsel der 
Kombinationen aus Drehzahl und Ölvolumenstrom und den sich jeweils einstellenden stationären 
Temperaturzuständen zu erkennen. Der Düsenwinkel bleibt in diesem Fall unverändert. 


Am Verlauf der Öltemperatur in der Öldüse in Abbildung 6.2 ist erkennbar, dass die Öltemperatur 
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Abbildung 6.2: Exemplarischer Verlauf der Betriebsparameter und Temperaturen während einer 
Messung 


in der Öldüse geringfügigen Schwankungen unterliegt. Diese Schwankungen treten vor allem 
bei einer Änderung des Ölvolumenstroms auf. Die Erhöhung des Ölvolumenstroms hat eine 
Erhöhung der Öltemperatur zur Folge, wie durch den vergrößerten Ausschnitt in Abbildung 6.2 
verdeutlicht wird. Dies liegt daran, dass die Verweilzeit in der Zuleitung, die stellenweise nicht 
isoliert ist, bei höheren Ölvolumenströmen kürzer ist und somit das Ol zwischen Tank und Oldiise 
weniger abkühlt. Da die gemessene Öltemperatur in der Berechnung des Wärmeübergangs 
berücksichtigt wird und die Viskositätsänderung für die beobachteten Temperaturschwankungen 
sehr gering ist, kann dieser Effekt vernachlässigt werden, sofern die Öltemperatur stationäre 
Zustände erreicht. 


Im Rahmen der exemplarischen Auswertung soll der Betriebspunkt 80 (BP80) betrachtet werden, 
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dessen Messzeitpunkt in Abbildung 6.2 eingezeichnet und durch die in Tabelle 6.1 aufgelisteten 
Parameter bestimmt ist. Aufgrund des vergleichsweise langen Zeitraums der gesamten Messung 
im Vergleich zur Messung eines Betriebspunkts ist die Messung von BP80 zu einem Zeitpunkt 
dargestellt. 


Tabelle 6.1: Parameter BP80 


Parameter Sollwert Istwert 


Volumenstrom 3,21/min 3,257 l/min 


Öltemperatur 80°C 79,87 °C 
Drehzahl 1500/min  1492/min 
Düsenwinkel 0° —0,3° 
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Abbildung 6.3: Gemessene Temperaturen auf der Zahnkontur für BP80 


Für Betriebspunkt 80 wird die in Abbildung 6.3 dargestellte Temperaturverteilung gemessen, 
wobei die Temperaturen an Messpunkten gleicher axialer Position als Linien dargestellt sind. Der 
Temperaturverlauf entlang der Koordinate s* zeigt für alle vier Messpunktzeilen einen qualitativ 
ähnlichen Verlauf. Die maximale Temperatur tritt für alle Zeilen an der Reihe 3, an der Position 
des größten Radius auf der strömungszugewandten Zahnflanke, auf. Eine Ausnahme bildet 
der Messpunkt 5 der Referenzmesszeile (Zeile 4), an dem geringfügig höhere Temperaturen 
gemessen werden als an Messpunkt 3. Das mittlere Temperaturniveau ist am höchsten für die 
Messreihe mit dem geringsten axialen Abstand zur Öldüse. Dies ist durch das dort aufprallende 
Öl und die damit verbundenen, hohen Wärmeübergänge zu begründen. 


Zur Deutung der Abweichung des Messpunkt 5 der Referenzmesszeile kommen mehrere Ur- 
sachen in Frage. Eine Asymmetrie der Wärmeübergangsverteilung ist grundsätzlich denkbar. 
Allerdings wird eine deutliche Temperaturabweichung nur an diesem Messpunkt beobachtet. In 
Messungen bei anderen Betriebspunkten werden ebenfalls Temperaturabweichungen beobachtet, 
diese fallen jedoch meist geringer aus. Des weiteren besteht die Möglichkeit, dass das Thermo- 
element dieser Messstelle eine Abweichung in seiner Einbettungstiefe ausweist. Eine solche 
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Abweichung ist jedoch allenfalls in einem geringen Maße denkbar, da die Thermoelementnut 
nur geringfügig größer ist als der Durchmesser des Thermoelements und ein Herausragen des 
Thermoelements deutlich sichtbar wäre. Dies lässt entweder eine herstellbedingte Abweichung 
im verwendeten Thermoelement oder in der Signalkette der Telemetrie als wahrscheinlichere 
Ursache erscheinen. 


Um aus den gemittelten Temperaturdaten zunächst Wärmeströme im Testzahnrad und anschlie- 
ßend, mit bekannter Öltemperatur, Wärmeübergänge zu berechnen, wird das in Kapitel 4 be- 
schriebene iterative Verfahren angewendet. Dabei werden die Messwerte auf einem Achtelzahn, 
der Hälfte eines periodischen Segments, zur Berechnung verwendet. Da die Interpretation eines 
ganzen Segments jedoch wesentlich anschaulicher ist, wird in allen folgenden Darstellungen die 
Wärmeübergangsverteilung entlang der Mitte eines Segments gespiegelt und zwischen yt = 0,25 
und y* = 0,5 betrachtet. Die zweidimensionale Verteilung des Wärmeübergangs für BP80 ist in 
Abbildung 6.4 dargestellt. 


Wie aus der Interpretation der Temperaturen bereits vermutet, Konzentriert sich das Wärmeüber- 
gangsmaximum an der axialen Position der Düse und an der strömungszugewandten Flanke. Dort 
werden Werte im Bereich von 1000 W/(m? K) erreicht. Entlang der Symmetrielinie yt = 0,375 
nimmt der Wärmeübergang zum Zahnfuß hin ab und erreicht dort Werte von ca. 600 W/(m? K). 
Auf der stromungsabgewandten Zahnflanke werden insgesamt geringere Werte des Wärmeüber- 
gangs gemessen, dessen Anteil am gesamten Wärmeübergang jedoch nicht zu vernachlässigen 
ist. Der Wärmeübergang in axialer Richtung ist gleichmäßiger verteilt als auf der strömungs- 
zugewandten Zahnflanke. Dies kann ein Hinweis auf das Abprallen des Strahls sein, das von 
Keller et al. (2017) beobachtet wurde. 


Negative Werte für h werden in flächenmäßig kleinen Bereichen in der Nähe des Zahnkopfes, 
außerhalb der Mitte beobachtet. Das Auftreten negativer Wärmeübergänge ist unphysikalisch, 
da diese einen Wärmestrom in Richtung eines positiven Temperaturgradienten erfordern und 
somit den zweiten Hauptsatz der Thermodynamik verletzen würden. Im Fall der vorliegenden 
Untersuchungen ist die Ursache dieser negativen Wärmeströme jedoch in der vergleichsweise 
geringen Ortsauflösung der Messpunkte zu suchen. Bereiche geringen Wärmestroms werden 
durch die Interpolation unter Umständen überschätzt. Aus Wärmebilanzgründen ergibt sich 
dadurch an anderen Stellen ein negativer Wärmeübergang. Darüber hinaus wird durch die 
verwendete Definition des Wärmeübergangs der Einfluss der im Prüfstand vorhandenen Luft nicht 
berücksichtigt. Durch Abschätzungen konnte gezeigt werden, dass dieser Einfluss im Vergleich 
zur Prallkühlung gering ist. Jedoch besteht die Möglichkeit, dass die Luft durch Abkühlung an 
den metallischen Wänden des Prüfstands eine geringere Temperatur aufweist als das Öl. Wenn 
das Öl die strömungsabgewandte Flanke nicht erreicht und somit keinen Wärmestrom in den 
Zahn hervorruft, kann der Wärmeübergang durch Luft, der zu einer Wärmeabfuhr führt, zur 
Berechnung von negativen Wärmeübergangskoeffizienten führen. 
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Abbildung 6.4: Wärmeübergangsverteilung für BP80 in W/(m?K) 


6.2 Charakterisierung des Ol-Zahnrad-Warmeiibergangs 


Ziel der vorliegenden Arbeit ist es, den Wärmeübergang an einem einspritzgekühlten Zahnrad 
experimentell zu ermitteln und den Einfluss der Betriebsparameter auf den Wärmeübergang 
zu quantifizieren. Die im Vorfeld bekannten Zusammenhänge und das verfügbare Wissen zur 
Öl-Zahnrad-Interaktion sind äußerst begrenzt. Die häufig herangezogene Wärmeübergangs- 
korrelation von DeWinter und Blok (1974) bildet einzig die Drehzahl als Einflussgröße auf 
den Wärmeübergang ab. Aus Arbeiten von Townsend und Akin (1981), Long et al. (2003) 
und weiteren ist jedoch bekannt, dass der Ölvolumenstrom einen erheblichen Einfluss auf den 
Wärmeübergang hat. 


Um sämtliche Einflüsse der mit dem verwendeten Prüfstand variierbaren Parameter zu erfassen, 
wird eine Vielzahl von Messungen unterschiedlicher Betriebspunkte vorgenommen. Aus Gründen 
der Nachvollziehbarkeit sollen die jeweiligen Parameter nacheinander diskutiert werden. Dabei 
können nur einige Messungen detailliert betrachtet werden. Für eine eingehende Darstellung 
aller Einzelergebnisse sei daher auf den Anhang verwiesen. 


Zunächst soll in Abschnitt 6.2.1 der Einfluss der Drehzahl und des Ölvolumenstroms auf den 
Wärmeübergang diskutiert werden. Anschließend wird in Abschnitt 6.2.2 die Abhängigkeit 
des Wärmeübergangs vom Düsenwinkel betrachtet. Die diskutierten Messungen in diesen 
beiden Abschnitten beziehen sich dabei stets auf Oltemperaturen von 80 °C. Da bei geänderten 
Temperaturen und somit geänderten Prandtlzahlen auch mit geänderten Strömungsvorgängen und 
Wärmeübergängen zu rechnen ist, wird der Einfluss der Prandtlzahl in Abschnitt 6.2.3 ebenfalls 
untersucht. 


Um die charakteristischen Abhängigkeit des Wärmeübergangs von den Betriebsparametern in 
einer möglichst allgemeinen Form zu analysieren, werden die dimensionslosen Größen Ge- 
schwindigkeitsverhältnis ©, Eindringwinkel B und die Relativgeschwindigkeit u, herangezogen. 
Aus der Analyse im Folgenden wird deutlich, dass die wesentlichen Wirkzusammenhänge durch 
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diese Größen erfasst werden. 


Die Beschreibung des charakteristischen Verhaltens des Wärmeübergangs alleine durch diese 
drei Größen greift jedoch zu kurz. Die Drehzahl und der Ölvolumenstrom wirken sich auf 
eine Art auf die gemessenen Wärmeströme aus, die über ihre Bedeutung beispielsweise für © 
hinausgeht. Wie in Abschnitt 3.2 beschrieben, ist zum einen zwischen anliegender und abgelöster 
Strömung zu unterscheiden. Außerdem beeinflusst die Drehzahl das Strömungsfeld um das 
Zahnrad. Daher werden diese beiden Größen ebenfalls in der folgenden Analyse herangezogen. 
Der Ölvolumenstrom wird durch die Strahlreynoldszahl 
pusd usd 


Re; = => wl 
eg 2 = (6.1) 


beschrieben. Die Drehzahl wird durch die Rotationsreynoldszahl 


Reg = —— (6.2) 


ausgedriickt. Als charakteristischer Durchmesser D wird dabei der Kopfkreisdurchmesser ver- 
wendet. 


Die verwendete Methode zur Messung des Öl-Zahnrad-Wärmeübergangs liefert ortsaufgelöste 
Daten. Um den Einfluss einzelner Parameter auf den Wärmeübergang zu charakterisieren, ist es 
sinnvoll, auch den gemittelten Wärmeübergang zu betrachten, als 


Nu= st,yH)aytast 63 
3 rl aN Jays 9a) 


mit dem Normierungsfaktor Nu*, der der höchsten im Rahmen dieser Arbeit gemessenen mittle- 
ren Nußeltzahl entspricht. 


6.2.1 Einfluss der Drehzahl und des Ölvolumenstroms 


Um den Wärmeübergang am Zahnrad zu charakterisieren, soll zunächst der in Messungen 
beobachtete Einfluss der Drehzahl (Reo) und des Ölvolumenstroms (Re,) auf die Nußeltzahl- 
verteilung betrachtet werden. In Abbildung 6.5 ist der Verlauf von Nu* bei Reg = 2639 und 
œ = 20° über einen Bereich von Re, zwischen 48 - 10* und 190- 10* dargestellt. Zusätzlich 
sind das Geschwindigkeitsverhältnis & und der Aufprallwinkel p über den betrachteten Bereich 
eingezeichnet. Durch die Beschreibung und Diskussion der vier markierten Betriebspunkte von 
niedrigen zu hohen Drehzahlen soll der Verlauf von Nu* im Folgenden erklärt werden. Dazu 
werden ergänzend Einzelbilder der Gegenlichtaufnahmen in Abbildung 6.6 herangezogen. 


In Abbildung 6.5 sind für zwei Betriebspunkte Wiederholmessungen verfügbar. Diese weisen 
untereinander geringfügige Abweichungen von Nu* auf. Diese Abweichungen haben ihren 
Ursprung in versuchsbedingten Einflüssen, wie beispielsweise leicht unterschiedliche Ölvor- 
lauftemperaturen durch Messungen zu unterschiedlichen Zeitpunkten. Für eine besser lesbare 
Darstellung des Kurvenverlaufs wird der interpolierte Verlauf von Nu* als Linie dargestellt. An 
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Abbildung 6.5: Einfluss der Drehzahl Re, auf den flächengemittelten Wärmeübergang Nu* bei 
Reg = 2639 und den ortsaufgelösten Wärmeübergang Nu 
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Abbildung 6.6: Gegenlichtaufnahmen bei Reg = 2639 
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Betriebspunkten, an denen mehrere Messungen vorliegen, wird der Mittelwert aller verfiigbarer 
Messungen als Grundlage fiir die Interpolation herangezogen. 


Zwei Einzelbilder der Gegenlichtaufnahme bei Rew = 48 - 104 (BP33) sind in Abbildungen 6.6a 
und 6.6b gezeigt. Bedingt durch das vergleichsweise hohe Geschwindigkeitsverhältnis mit 
E = 1,4 und den positiven Düsenwinkel von & = 20° ergibt sich ein steiler Einfall des Ölstrahls 
mit einem Aufprallwinkel $ von 22,49°. Dies ist daran zu erkennen, dass der Ölstrahl zu ei- 
nem Zeitpunkt kurz vor dem Erreichen der strömungszugewandten Flanke bereits den Zahnfuß 
erreicht hat. Das untere Drittel der Zahnlücke ist zu diesem Zeitpunkt mit Öl gefüllt. Der 
ortsaufgelöste Wärmeübergang bei Rew = 48 - 10* (BP33) ist in Abbildung 6.5b dargestellt. Es 
wird eine Wärmeübergangsverteilung gemessen, die am Zahnfuß konzentriert ist und somit 
die Beobachtungen aus den Gegenlichtaufnahmen bestätigt. Dort werden lokale Nußeltzahlen 
von etwa 30 beobachtet. Dies kann auf das steile Eindringen des Ölstrahls und das dadurch 
hervorgerufene Auftreffen des Öls am Zahnfuß begründet werden. In lateraler Richtung werden 
in der Nähe der Position der Düse, y* = 0,375, höhere Wärmeübergangskoeffizienten als an 
den Rändern des periodischen Elements gemessen. Auf der strömungszugewandten Zahnflanke 
sind höhere Wärmeübergangskoeffizienten zu beobachten als auf der gegenüberliegenden, strö- 
mungsabgewandten Flanke. Die Ursache für diesen Unterschied kann in der höheren Menge Öl 
gesehen werden, die aufgrund der kinematischen Verhältnisse auf der strömungszugewandten 
Zahnflanke auftrifft. Im Bereich des Zahnkopfes wird das geringste Niveau des Wärmeübergangs 
beobachtet. Stellenweise treten hier leicht negative Werte auf. 


Bei Erhöhung der Drehzahl (siehe Verlauf in Abbildung 6.5a) nimmt einerseits das Verhältnis 
aus Strahl- und Tangentialgeschwindigkeit È stetig ab, während der Aufprallwinkel ß zunimmt 
(siehe Abbildung 3.5). Dies ist gleichbedeutend mit einem flacheren Auftreffen des Ölstrahls 
auf die Zahnoberfläche und einer geringeren radialen Eindringtiefe. Andererseits ist durch den 
erhöhten spezifischen Impuls des auftreffenden Öls mit erhöhtem Wärmeübergang zu rechnen. 
Ausgehend von Rew = 48 - 104 (BP33) ist bei Erhöhung der Drehzahl zunächst ein Anstieg von 
Nu* zu beobachten. In diesem Bereich überwiegt demnach der Effekt des erhöhten Impulses. 
Bei Rew = 118- 10* (BP68, Abbildung 6.5c) ist das örtliche Maximum des Wärmeübergangs 
durch die Änderung der kinematischen Verhältnisse in Richtung des Zahnkopfes und in positive 
s*-Richtung verschoben und mit etwa Nu = 40 deutlich höher. Aus den Gegenlichtaufnah- 
men in Abbildungen 6.6c und 6.6d ist ersichtlich, dass der Ölstrahl an diesem Betriebspunkt 
flacher in die Zahnlücke eindringt und den Zahnfuß nicht erreicht. In lateraler Richtung ist, 
ausgehend vom Bereich des Maximums auf der strömungszugewandten Seite, ein starker Ab- 
fall des Wärmeübergangs zu den Rändern des Bereichs zu beobachten. Auf der unbelasteten 
Zahnflanke ist jedoch ein Bereich bei etwa s* = 0,8 zu erkennen, an dem entlang des gesamten 
y*-Bereichs nahezu konstante Nußeltzahlen gemessen werden. Dies ist möglicherweise auf das 
Abprallen des Ölstrahls an der Aufprallstelle, anschließendes Auffächern und erneutes Auftreffen 
auf der strömungsabgewandten Seite zurückzuführen. In den Abbildungen 6.6c und 6.6d ist 
zur Veranschaulichung dieses Zusammenhangs die Zahnkontur der linken benachbarten Zahn- 
lücke hervorgehoben. Hier ist zu erkennen, dass das Öl etwa die untere Hälfte der Zahnlücke 
füllt (Abb. 6.6c) und sich auch in der oberen Hälfte auf die strömungsabgewandte Zahnflanke 
zubewegt (Abb. 6.6c), wie durch die beiden Pfeile angedeutet. 
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Bei weiterer Erhöhung der Drehzahl nimmt Nu* rapide ab. Der Maximalwert Nu* von ca. 0,78 
des mittleren Wärmeübergangs für Re, = 2639 und & = 20° bei Reg © 118 - 10* (BP68) sinkt auf 
etwa 0,57 bei Reo = 143 10* (BP43). In Abbildung 6.5d ist zu erkennen, dass das Maximum des 
Wärmeübergangs weiter in Richtung des Zahnkopfes verschoben ist. Sowohl das Niveau als auch 
die örtliche Ausdehnung des Bereichs hohen Wärmeübergangs auf der strömungszugewandten 
Seite sind vergleichbar mit der Verteilung des Wärmeübergangs für BP68 in Abbildung 6.5c. Auf 
der strömungsabgewandten Zahflanke wird jedoch ein deutlicher Unterschied der Verteilungen 
gemessen. In der Nähe des Zahnkopfes ist der Wärmeübergang auf etwa ein Drittel des Wertes 
von BP68 reduziert, in weiten Bereichen des Zahnfußes ist eine Reduktion auf etwa die Hälfte zu 
beobachten. Als mögliche Erklärung können Abbildungen 6.6e und 6.6f herangezogen werden. 
Dort ist ein ölfreier Spalt an der strömungsabgewandten Zahnflanke zu erkennen, während bei 
BP68 beobachtet wurde, dass Öl die strömungsabgewandte Zahnflanke erreicht (vgl. Abb. 6.6c 
und 6.6d). Durch das Abreißen der Olzufuhr an die strömungsabgewandte Zahnflanke kann hier 
auch keine Kühlwirkung hervorgerufen werden. 


Bei Erhöhung der Drehzahl auf Reo = 191- 10* (BP49) setzt sich die Verschiebung des Bereichs 
hohen Wärmeübergangs in Richtung des Zahnkopfes fort. Dieser ist nun stark in der Nähe des 
Zahnkopfes auf der strömungszugewandten Seite konzentriert (siehe Abbildung 6.5e). Durch die 
Erhöhung des Impulses ist auch das lokale Maximum des Wärmeübergangs erhöht. Während 
dort Maximalwerte von etwa Nu = 50 gemessen werden, ist auf der strömungsabgewandten Seite 
der Wärmeübergang zum Erliegen gekommen. Diese Beobachtungen stellen eine Fortsetzung der 
zuvor beobachteten Effekte dar. Aus Abbildungen 6.6g und 6.6h ist ersichtlich, dass der Ölstrahl 
bei BP49 noch weniger tief in die Zahnlücke eindringt und dass die strömungsabgewandte 
Zahnflanke nahezu frei von Öl bleibt. 


Analog zur vorangegangenen Betrachtung und Abbildung 6.5 soll im Folgenden der drehzahl- 
abhängige Verlauf von Nu* bei reduziertem Ölvolumenstrom und einer Strahlreynoldszahl von 
Re, = 1757 diskutiert werden. Dazu werden der Verlauf des mittleren Wärmeübergangs und 
die Wärmeübergangsverteilungen an vier Betriebspunkten für Reg = 1757 in Abbildung 6.7 
betrachtet. Insgesamt liegt der Wärmeübergang auf einem geringeren Niveau. Während bei 
Reg = 2639 und Rew = 118 - 10* maximale Werte von Nu* œ 0,78 erreicht werden, werden bei 
Reg = 1757 lediglich mittlere Nußeltzahlen von ca. 0,6 bei Rew ~ 40- 10* gemessen. 


Durch die geringere Strahlgeschwindigkeit bei Rey = 1757 ist der Aufprallwinkel ß bei gleichen 
Werten für Reo stets höher als in der vorangegangen Betrachtung bei Reg = 2639. Dies äußert 
sich in einem Abfall von Nu* bei geringeren Drehzahlen, welcher zwischen Rew = 48 - 10* und 
Reo = 95 - 10% beobachtet wird. 


In Abbildung 6.7b ist der ortsaufgelöste Wärmeübergang Reg = 48 - 10* (BP53) dargestellt. Der 
Aufprallwinkel an diesem Betriebspunkt ist mit ß = 38,44° klein, es ist mit einem steilen Einfall 
des Ölstrahls zu rechnen. Dies führt zu einer Konzentration des Wärmeübergangs am Zahnfuß, 
wie sie auch bei Re, = 2639 und kleinen Aufprallwinkeln beobachtet wird. 


Bei Erhöhung der Drehzahl auf Reg = 71- 10* (BP29, siehe Abb. 6.7a) und einem daraus 
resultierenden Aufprallwinkel von ß = 54,32° wird die in Abbildung 6.7c dargestellte Wär- 
meübergangsverteilung gemessen. Entlang der Linie der axialen Düsenposition, y” = 0,375, 
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Abbildung 6.7: Einfluss der Drehzahl Reo auf den flächengemittelten Wärmeübergang Nu* bei 
Reg = 1757 und den ortsaufgelösten Wärmeübergang Nu 
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werden auf der strömungszugewandten Flanke Wärmeübergangskoeffizienten von Nu ~ 20 
gemessen. In Richtung des Zahnkopfes ist eine leichte Erhöhung des Wärmeübergangs im 
Vergleich zum Zahnfuß festzustellen. Diese Tendenz unterscheidet sich qualitativ von der Wär- 
meübergangsverteilung für BP68 (siehe Abbildung 6.5c), an dem mit ß = 57,95° ein nahezu 
gleicher Aufprallwinkel vorliegt. Dort wird eine deutliche Konzentration des Wärmeübergangs 
in der Zahnmitte beobachtet. 


Eine mögliche Ursache dieses Unterschieds liegt in der Strömungsform des Ölstrahls, der 
bei Re, = 1757 eine glaskörperartige, gleichmäßige Struktur aufweist, die charakteristisch 
für laminare, anliegende Düsenströmungen ist. Bei Reg = 2639 wird eine deutlich welligere, 
diffusere Struktur beobachtet, die auf einen abgelösten Strahl hinweist. In diesem Fall liegt durch 
die Einschnürung der Strömung eine höhere Strahlgeschwindigkeit vor, als sie die Berechnung 
anhand 

V 


Us 


erwarten lässt. Folglich liegt bei abgelöster Strömung tatsächlich ein höherer Aufprallwinkel als 
der berechnete vor. 


Der Wärmeübergang auf der strémungsabgewandten Zahnflanke bei Rew = 71 - 104 (BP29, s. 
Abb. 6.7c) ist im Vergleich zu Rew = 48 - 104 bereits deutlich reduziert. Bei weiterer Erhöhung 
auf Reo = 95 - 10* (BP54, s. Abb. 6.7d) kommt es zum beinahe vollständigen Erliegen des 
Wärmeübergangs auf der strömungsabgewandten Flanke. Im Bereich des Zahnkopfes auf der 
strömungsabgewandten Seite werden geringfügige, negative Wärmeübergangskoeffizienten 
gemessen. 


Bei Reg = 95- 10* und B = 62,89° findet der Wärmeübergang nahezu ausschließlich am Zahn- 
kopf der strömungszugewandten Flanke statt. Im Bereich des Zahnfußes sind nur noch geringe 
Nußeltzahlen von ca. 10 zu beobachten. Bei weiterer Erhöhung der Drehzahl auf Rew = 143 - 10* 
(BP56, s. Abb. 6.7e) ist ausschließlich am Zahnkopf signifikanter Wärmeübergang zu beobachten. 
Dies ist auf das sehr geringe Eindringen des Ölstrahls mit einem Aufprallwinkel von ß = 72,21° 
zurückzuführen, bei dem ausschließlich die strömungszugewandte Zahnflanke in der Nähe des 
Zahnkopfes vom Öl erreicht wird. 


Die vorangegangene Betrachtung des Wärmeübergangs bei Reg = 2639 und Reg = 1757 soll im 
Folgenden in knapper Form bei erhöhtem Ölvolumenstrom mit Reg = 3954 angestellt werden. 
In Abbildung 6.8a ist der Verlauf des mittleren Wärmeübergangs bei Reg = 3954 in einem 
Drehzahlbereich von Rew = 48 - 104 bis 240 - 10* dargestellt. Erneut ist ein Bereich geringfügig 
ansteigender Werte von Nu* zu erkennen, der sich von Reo = 48 - 10* bis 145-104 bzw. B = 
50,14° erstreckt. Bei höheren Drehzahlen fällt der mittlere Wärmeübergang ab. Insgesamt ist 
der mittlere Wärmeübergang auf einem höheren Niveau als für Reg = 2639. 


Bei Reo = 71: 10* (BP38, s. Abb 6.7b) und ß = 23,05° ist der Wärmeübergang im Bereich 
des Zahnfußes konzentriert. Entlang der gesamten Linie st = 0,2 nimmt der Wärmeübergang 
hohe Werte von bis zu Nu = 40 an. Auch auf der stromungsabgewandten Seite ist ein breiter 
Bereich hohen Wärmeübergangs zu beobachten. Dies entspricht den beobachteten Verteilungen 
für kleine Aufprallwinkel für Reg = 1757 und Reg = 2639. 
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Mit Erhöhung der Drehzahl auf Reg = 143- 10* (BP36, s. Abb. 6.8c) beträgt der kinematische 
Aufprallwinkel B = 50,14°. Wie auch bei den anderen betrachteten Ölvolumenströmen bzw. Wer- 
ten für Re, führen diese kinematischen Bedingungen zu einem Maximum des Wärmeübergangs 
im unteren Bereich der strömungszugewandten Flanke. Das Niveau des Wärmeübergangs ist 
im Vergleich zu Betriebspunkten bei niedrigeren Werten von Re, und Reg, die einen ähnlichen 
Aufprallwinkel aufweisen, erhöht (siehe Abbildung 6.8c). 


Bei Reo = 180-104 (BP18, s. Abb. 6.8d) ist bereits ein geringfügiger Rückgang von Nu* im 
Vergleich zu BP38 zu beobachten, der durch den Rückgang der örtlichen Wärmeübergangs- 
koeffizienten auf der strömungsabgewandten Seite und auch in den äußeren Bereichen des 
Zahnfußes entsteht. Bei Erhöhung von Reg (BP19, s. Abb. 6.8d) auf 240- 10* und P = 66,48° 
setzt sich diese Tendenz fort und der Wärmeübergang auf der stromungsabgewandten Seite trägt 
nur noch marginal zum gesamten Wärmeübergang bei. Bei dieser Drehzahl ist der maximale 
Wärmeübergang, wie bei anderen Betriebspunkten mit ähnlichem Aufprallwinkel, im Bereich 
des strömungszugewandten Zahnkopfes konzentriert. Es ist anzunehmen, dass diese Tendenz 
erneut durch das flachere Eindringen des Ölstrahls und die geringere Ölmenge, die bei großen 
Aufprallwinkeln die strömungsabgewandte Zahnflanke erreicht, hervorgerufen wird. 
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Abbildung 6.8: Einfluss der Drehzahl Reo auf den flächengemittelten Wärmeübergang Nu* bei 
Reg = 3954 und den ortsaufgelösten Wärmeübergang Nu 


88 Ergebnisse 


Die im Vorangegangenen diskutierten Abhängigkeiten des Öl-Zahnrad-Wärmeübergangs zeigen 
eindeutige charakteristische Abhängigkeiten von der Strahlreynoldszahl Re, und der Drehzahl 
Rew. In Abbildung 6.9 sind die Verläufe des mittleren Wärmeübergangs Nu* für die diskutierten 
Ölvolumenströme mit Reg = 1757, 2639 und 3954 sowie Ergebnisse der dazwischen liegenden 
Werte von Re, = 2251 und 3295 zusammengefasst dargestellt. Der Düsenwinkel beträgt stets 
a = 20°. Es ist ersichtlich, dass der Verlauf von Nu* über Reg bei Rey = 2251 und 3295 
qualitativ mit den bereits diskutierten Verläufen übereinstimmt. Die ortsaufgelösten Messungen 
des Wärmeübergangs, die hier nicht detailliert betrachtet werden sollen, zeigen die gleichen 
charakteristischen Muster und Vorgänge mit steigenden Drehzahlen, die zunächst zu einem 
Anstieg und anschließend zu einem Abfall des mittleren Wärmeübergangs führen. Steigende 
Werte für Re, führen dabei zu steigenden Werten von Reg, bei denen der Abfall von Nu* eintritt. 
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Abbildung 6.9: Verlauf von Nu* in Abhängigkeit von Reo und Reg und für & = 20° 


Fazit 


In der vorangegangenen Analyse der Ergebnisse wurde erstmalig der Wärmeübergang an einem 
einspritzgekühlten Zahnrad isoliert von anderen Wärmeströmen und ortsaufgelöst beschrieben. 
Im betrachteten Betriebsbereich sind der eingespritzte Olvolumenstrom, beschrieben durch Reg, 
und der mittlere Warmeiibergang Nu* stets positiv korreliert (siehe Abbildung 6.9). Zwischen 
der Drehzahl und dem Wärmeübergang besteht eine komplexe, nicht-monotone Abhängigkeit, 
die von einem deutlichen Rückgang von Nu* ab einer kritischen Drehzahl und der Verschiebung 
des Maximums des Wärmeübergangs vom Zahnfuß zum Zahnkopf gekennzeichnet ist. Durch die 
vereinfachte Betrachtung des Vorgangs als kinematisches Modell kann diese Verschiebung erklärt 
werden. Bei geringen Drehzahlen ist der Aufprallwinkel, unter dem der Ölstrahl in die Zahnlücke 
eindringt, gering. Dadurch erreicht das Öl den Zahnfuß und bewirkt hier ein Maximum des 
Wärmeübergangs. Mit steigenden Drehzahlen wird der Aufprallwinkel größer. Der Zahnfuß 
wird ab einem gewissen Aufprallwinkel nicht mehr erreicht und das Maximum des Wärme- 
übergangs verschiebt sich in Richtung des Zahnkopfes. Die Deutung durch den Aufprallwinkel 
wird durch die Gegenlichtaufnahmen bei Re, = 2639 bestätigt. Der Rückgang von Nu* kann 
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durch die Gegenlichtaufnahmen ansatzweise erklärt werden. Bei Aufprallwinkeln von bis zu 
B ~ 55° prallt Öl an der strömungszugewandten Zahnflanke ab, trifft auf die gegenüberliegende, 
strömungsabgewandte Zahnflanke und führt dort zu signifikantem Wärmeübergang. Bei Auf- 
prallwinkeln von über B ~ 65° erreicht das abprallende Öl die stromungsabgewandte Zahnflanke 
nicht mehr, wodurch dort der Wärmeübergang zum Erliegen kommt. Bei den drei diskutierten 
Ölvolumenströmen tritt der Rückgang bei Werten zwischen B = 50° und 60° auf. Durch die 
Änderung der Strömungsform in der Düse für verschiedene Re, ist der kritische Winkel Ber, bei 
dem der Rückgang von Nu* eintritt, in Fälle anliegender und abgelöster Strömung zu unterteilen. 
Diese Fallunterscheidung wird in Abschnitt 6.3 genauer erläutert. 


6.2.2 Einfluss des Düsenwinkels 


Untersuchungen von Keller et al. (2017) zeigen, dass der Düsenwinkel einen deutlichen Einfluss 
auf den Aufprallvorgang hat. Die Messungen des Wärmeübergangs weisen eine deutliche 
Abhängigkeit vom Düsenwinkel auf, wie im Folgenden gezeigt werden soll. 


Zunächst soll der Verlauf des Wärmeübergangs für diskrete Düsenwinkel von & = 20°, wie 
in Abschnitt 6.2.1, und & = 45° bei Reg = 2639 und in einem Drehzahlbereich zwischen 
Re» = 48: 10* und 190 - 10* diskutiert werden. Anschließend wird der Einfluss der Variation 
des Düsenwinkels bei ansonsten gleichbleibenden Bedingungen diskutiert. 


In Abbildung 6.10 ist der Verlauf des mittleren Wärmeübergangs bei Re, = 2639 dargestellt. 
Insgesamt liegt der mittlere Wärmeübergang für beide Düsenwinkel auf einem ähnlichen Niveau. 
Im Gegensatz zu den bisher betrachteten mittleren Wärmeübergangsverläufen bei & = 20° wird 
bei & = 45° eine monotone Abnahme mit zunehmender Drehzahl beobachtet. Mit steigender 
Drehzahl nimmt Nu* stetig ab. 


Für Re, > 119: 10* sind die Verläufe von Nu* für beide Düsenwinkel etwa parallel, der mittlere 
Wärmeübergang ist für æ = 45° stets geringer als für & = 20°. In diesem Bereich liegen mit 
xX < 0,6 geringe Geschwindigkeitsverhältnisse vor. Dadurch ist der Einfluss des Düsenwinkels 
auf den Aufprallwinkel B gering (siehe Abbildung 3.7 in Abschnitt 3.3). 


Zur Erklärung des unterschiedlichen Wärmeübergangsniveaus kann die Darstellung des orts- 
aufgelösten Wärmeübergangs in Abbildung 6.11 herangezogen werden. Die Wärmeübergangs- 
verteilungen für Rey = 119- 10* bei & = 20° (BP68) und & = 45° (BP103) sind in den Abbil- 
dungen 6.11a und 6.11b dargestellt. In beiden Fällen sind qualitativ ähnliche Verteilungen auf 
der strömungszugewandten Flanke zu erkennen. Das Maximum des Wärmeübergangs wird im 
Bereich der axialen Mitte der strömungszugewandten Zahnflanke bei etwa st = 0,2 gemessen. 
Bei & = 45° ist der Bereich hohen Wärmeübergangs stärker in der Mitte konzentriert, während 
er bei & = 20° etwas gleichmäßiger verteilt ist. Auf der strömungsabgewandten Zahnflanke 
werden ebenfalls qualitativ ähnliche Wärmeübergangsverteilungen gemessen. Jedoch ist das 
Niveau bei & = 45° deutlich geringer als bei & = 20°. Dieses geringere Niveau ist ursächlich für 
den geringeren mittleren Wärmeübergang. 


Eine mögliche Ursache für diesen Niveauunterschied kann im unterschiedlichen Aufprall des 
Ölstrahls auf die Zahnflanke gesehen werden. Mit steigenden Düsenwinkeln sinkt die Neigung 
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Abbildung 6.10: Einfluss des Diisenwinkels & auf den flächengemittelten Wärmeübergang Nu* 
bei Reg = 2639 


des Ölstrahls, in Richtung der strömungsabgewandten Zahnflanke abzuprallen. Somit ist die 
Kühlwirkung an der strömungsabgewandten Zahnflanke bei zunehmendem œ geringer. Dies führ- 
te bereits bei der Betrachtung des drehzahlabhängigen Wärmeübergangs bei hohen Drehzahlen 
zu einem Rückgang von Nu*. 


Bei einer auf Reo = 167: 10* gesteigerten Drehzahl wird ein ähnlicher Effekt beobachtet. Die 
Wärmeübergangsverteilungen bei dieser Drehzahl sind in den Abbildungen 6.11c und 6.11d 
dargestellt. Bedingt durch den hohen Aufprallwinkel von etwa B = 68° ist der Wärmeübergang 
wie in anderen betrachteten Fällen um den Zahnkopf konzentriert. Erneut ist die Verteilung von 
Nu bei & = 20° stark auf der strömungszugewandten Zahnflanke konzentriert und dort etwas 
gleichmäßiger verteilt. Der Bereich hohen Wärmeübergangs erstreckt sich bis zum Zahnfuß, 
während er bei œ = 45° dort sehr gering ist. Bei & = 20° wird auf der strömungsabgewandten 
Zahnflanke noch ein geringfügiger Wärmeübergang gemessen, während er bei & = 45° dort 
nahezu zum Erliegen kommt. Wie bei Reg = 119 - 10* kann das geänderte Abprallverhalten des 
Öls für die verschiedenen Düsenwinkel als Ursache für diesen Unterschied gesehen werden. Das 
Öl erreicht bei & = 45° die strömungsabgewandte Zahnflanke nicht. 


Im niedrigeren Bereich von Rew zwischen 48- 10* und 94 - 10* werden große qualitative Un- 
terschiede des mittleren Wärmeübergangs für & = 45° und & = 20° gemessen, wie aus Ab- 
bildung 6.10 ersichtlich ist. Während bei & = 20° steigende Drehzahlen gleichbedeutend mit 
steigenden Werten von Nu* sind, sinkt der Wärmeübergang bei & = 45° in diesem Bereich 
kontinuierlich. 


Die Verteilungen des Wärmeübergangs bei & = 20° und & = 45° unterscheiden sich sowohl 


6.2 Charakterisierung des Öl-Zahnrad-Wärmeübergangs 91 


(a) BP68, Reo = 119: 10*, œ = 20°, B = 57,95° (b) BP103, Reo = 119: 104, œ = 45°, B = 58,33° 


(c) BP70, Rew = 167: 10*, œ = 20°, B = 67,22° (d) BP105, Rew = 167: 104, a = 45°, B = 69,52° 


Abbildung 6.11: Ortsaufgelöste Darstellung des Wärmeübergangs Nu bei Reg = 2639 und hohen 
Drehzahlen 


für Reg = 48 - 104 als auch für Reg = 71 10% stark, wie in Abbildung 6.12 deutlich wird. Bei 
a = 45° ist bei beiden Drehzahlen ein breiter Bereich in axialer Richtung in der Nähe des 
Zahnfußes auf der strömungszugewandten Flanke zu erkennen, der hohe Wärmeübergangsko- 
effizienten aufweist. Die Konzentration des Maximums in der Nähe des Zahnfußes entspricht 
den Erkenntnissen des vorangegangen Abschnitts für kleine Aufprallwinkel. Das Maximum in 
diesem Bereich ist jeweils geringfügig höher als das bei @ = 20°. Auch der Wärmeübergang auf 
der strömungsabgewandten Flanke ist erhöht, stellenweise um den Faktor zwei. Das geringfügig 
höhere Wärmeübergangsniveau wird durch den Aufprallwinkel erklärt, der für & = 45° bei glei- 
cher Drehzahl Reo stets kleiner ist als für & = 20°. Im Fall von BP108 in Abbildung 6.12b mit 
Rew = 48 - 104 und & = 45° beträgt der Aufprallwinkel B = 3,06°, was einem nahezu senkrechten 
Eindringen des Ölstrahls in die Zahnlücke entspricht. Dabei wird, bei rein kinematischer Be- 
trachtung des Vorgangs, auch die streömungsabgewandte Flanke direkt durch den Strahl getroffen. 
Dies kann als mögliche Ursache des erhöhten Wärmeübergangs an der strömungsabgewandten 
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Flanke gesehen werden. Die große axiale Ausdehnung des Bereichs hohen Wärmeübergangs 
kann ein Hinweis auf ein axiales Ausbreiten des Öls am Zahnfuß gewertet werden, das durch das 
senkrechte Eindringen des Öls hervorgerufen wird. Bei & = 20° und ansonsten gleichbleibenden 
Bedingungen (BP, s. Abb. 6.12a) beträgt der Aufprallwinkel B = 22,49° und ist somit deutlich 
flacher. In diesem Fall ist der Wärmeübergang an der strömungsabgewandten Flanke und seine 
axiale Ausbreitung am Zahnfuß geringer. 


Bei einer auf Re» = 71 - 10* erhöhten Drehzahl wird für & = 20° (BP34) die in Abbildung 6.12c 
dargestellte Wärmeübergangsverteilung gemessen. Hier ist im Vergleich zu Reg = 48 104 
(BP33) eine Verschiebung des Wärmeübergangs in Richtung des Zahnkopfes erkennbar, was 
erneut durch den erhöhten Aufprallwinkel von B = 39,08° erklärt wird. Bei Rew = 71: 104 
und & = 45° (BP107, s. Abb. 6.12d) beträgt der Aufprallwinkel ß = 29,97°, was zu einem 
vorwiegend am Zahnfuß konzentrierten Wärmeübergang führt. Das hohe lokale Maximum von 
Nu = 25 auf der strömungsabgewandten Flanke, das auch bei Reo = 48 - 104 gemessen wird, 
kann jedoch nicht durch den Aufprallwinkel erklärt werden. Möglicherweise ist die tatsächliche 
Strahlgeschwindigkeit us aufgrund eines abgelösten Ölstrahls höher als die berechnete, wodurch 
die strömungsabgewandte Zahnflanke auch bei Reo = 71: 10* und & = 45° vom Öl erreicht 
wird. 


Die Verteilungen des Wärmeübergangs bei & = 45° und Reg = 48 10* und 71 - 10* in Abbil- 
dung 6.12b und 6.12d weisen auf der strömungsabgewandten Flanke ein lokales Minimum bei 
yt = 0,375 und etwa s* = 0,83 auf. Auf Grundlage der verfügbaren Hochgeschwindigkeits- 
aufnahmen und aus Simulationen geht nicht hervor, wodurch dieser Effekt entsteht. Die Linie 
s" = 0,83 liegt jedoch im konkaven Bereich des ZahnfuBes und wird möglicherweise nicht von 
an der strömungsabgewandten Zahnflanke abprallendem Öl getroffen. Für eine Überprüfung 


dieser Beobachtung und dem unerwartet hohen Wärmeübergang auf der strömungsabgewandten 
Zahnflanke bei Reo = 71: 10* und & = 45° sind zusätzliche Informationen aus Experimenten 
oder numerischen Simulationen notwendig. 
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(a) BP33, Re, = 48: 104, & = 20°, B = 22,49° 
0,5 


(c) BP34, Re, = 71: 104+, a = 20°, B = 39,08° (d) BP107, Reg = 71: 104, œ = 45°, B = 29,97° 


Abbildung 6.12: Ortsaufgelöste Darstellung des Wärmeübergangs Nu bei Re, = 2639 und 
niedrigen Drehzahlen 


Um den Einfluss des Düsenwinkels eingehender zu betrachten, ist der Verlauf von Nu bei 
Reg = 2639 und Reg = 143: 10* im Bereich von & von —30° bis 45° in Abbildung 6.13 
dargestellt. Da Reg, Rew und somit È für die betrachteten Betriebspunkte konstant sind, ist der 
Aufprallwinkel B eine Funktion von & mit einem Minimum bei & ~ 25°. Im Bereich zwischen 
a = —30° und 20° ist der mittlere Wärmeübergang nahezu konstant. Bei einer Vergrößerung 
des Düsenwinkels darüber hinaus in Richtung positiver Werte ist ein Absinken des mittleren 
Wärmeübergangs zu beobachten. 


Durch die Analyse des ortsaufgelösten Wärmeübergangs kann auf die Ursachen des Verlaufs von 
Nu* in Abhängigkeit von & geschlossen werden. In Abbildung 6.14 sind die Wärmeübergangs- 
verteilungen der Betriebspunkte 47, 44, 43 und 104 dargestellt. Bei einem Düsenwinkel von 
a = —30° und B = 72,51° für BP47 in Abbildung 6.14a wird eine Wärmeübergangsverteilung 
gemessen, deren Maximum im Bereich des Zahnkopfes liegt. Dies entspricht der Erwartung aus 
den vorangegangenen Betrachtungen für ähnlich flache Aufprallwinkel. Bei einer Erhöhung des 
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Abbildung 6.13: Einfluss des Diisenwinkels & auf den flächengemittelten Wärmeübergang Nu* 
bei Re, = 2639 und Rew = 143 - 104 


Düsenwinkels auf & = 0° (BP44) ist nur eine geringfügige Erhöhung des mittleren Wärmeüber- 
gangs zu beobachten. Die Verteilung des Wärmeübergangs zeigt jedoch deutliche Unterschiede 
(siehe Abbildung 6.14b). Der Bereich hohen Wärmeübergangs erstreckt sich über einen größeren 
Bereich der strömungszugewandten Zahnflanke, wodurch der mittlere Wärmeübergang erhöht 
wird. Diese Verschiebung ist auf den steileren Aufprallwinkel und das damit verbundene tiefere 
Eindringen des Öls zurückzuführen. Auf der strömungsabgewandten Zahnflanke ist jedoch 
ein geringfügiger Rückgang des Wärmeübergangs zu verzeichnen, was einen gegenläufigen 
Effekt darstellt. Dieser Rückgang kann als Hinweis auf geändertes Abprallen des Ölstrahls 
und reduziertes Auftreffen auf der strömungsabgewandten Flanke gewertet werden. Die beiden 
gegenläufigen Effekte gleichen sich im Bereich von & = —30° bis 10° gegenseitig aus, wodurch 
Nu* etwa konstant bleibt. 


Der Effekt der Umverteilung des Wärmeübergangs setzt sich für größere Düsenwinkel fort. Aus 
der in Abbildung 6.14c dargestellten Wärmeübergangsverteilung für & = 20° (BP43) ist eine 
wesentlich gleichmäßigere Verteilung auf der stromungszugewandten Zahnflanke erkennbar. Der 
Rückgang des Wärmeübergangs an der stromungsabgewandten Zahnflanke überwiegt, es kommt 
zu einem beginnenden Rückgang von Nu*. Bei weiterer Drehung der Ölbrause in Richtung 
der Zahnradrotation zu & = 45° (BP104) wird eine Wärmeübergangsverteilung gemessen, die 
auf der gesamten Höhe der strömungszugewandten Zahnflanke hohe Werte aufweist und etwa 
die Form eines Rechtecks annimmt. Der Wärmeübergang auf der strömungsabgewandten Seite 
ist im Vergleich zum Wärmeübergang bei kleineren Düsenwinkeln deutlich reduziert. Der 
Rückgang des dortigen Wärmeübergangs verursacht maßgeblich die Reduktion des mittleren 
Wärmeübergangs Nu*. 
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(c) BP43, & = 20°, B = 63,80° (d) BP104, œ = 45°, B = 65,07° 


Abbildung 6.14: Einfluss des Düsenwinkels @ auf die Wärmeübergangsverteilung Nu bei 
Reg = 2639 und Reo = 143 - 104 


Fazit 


Die vorangegangene Betrachtung liefert die Erkenntnis, dass der Düsenwinkel sowohl den 
mittleren Wärmeübergang als auch die Verteilung des Wärmeübergangs beeinflusst. Größere 
Düsenwinkel führen dabei zu einer gleichmäßigeren Verteilung des Wärmeübergangs auf der 
strömungszugewandten Zahnflanke und tendenziell höheren Wärmeübergängen in diesem Be- 
reich. Auf der gegenüberliegenden, strömungsabgewandten Zahnflanke werden abnehmende 
Wärmeübergänge beobachtet. Im betrachteten Fall überwiegt dieser Rückgang ab einem Düsen- 
winkel von & = 20°, was zu geringeren Werten von Nu* führt. Die Vor- und Nachteile eines auf 
die strömungszugewandte Zahnflanke konzentrierten Wärmeübergangs sind vom Anwendungs- 
fall abhängig und können im Rahmen dieser Arbeit nicht bewertet werden. Die diskutierten 
Messungen zeigen jedoch, dass die Wärmeübergangsverteilung durch den Düsenwinkel gezielt 
beeinflusst werden kann. Insgesamt ist die Abhängigkeit des mittleren Wärmeübergangs vom 
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Diisenwinkel wesentlich geringer als seine Abhängigkeit von der Drehzahl und dem Ölvolumen- 
strom. 


6.2.3 Einfluss der Temperatur 


Um den Einfluss der Öltemperatur auf den Wärmeübergang zu quantifizieren, wird im Folgenden 
der Verlauf von Nu* in Abhängigkeit der Drehzahl bei & = 20° und einem Ölvolumenstrom 
von V = 4,81/min bei 3, = 60°C mit dem Verlauf von Nu* bei 9, = 80°C und ansonsten 
identischen Parametern verglichen. Die Messungen bei V, = 80°C wurden im vorangegangenen 
Abschnitt diskutiert. In diesem Vergleich sind © und ß für beide Temperaturen gleich, da 
der Ölvolumenstrom in beiden Fällen gleich ist. Die Strahlreynoldszahlen sind aufgrund der 
temperaturabhängigen Viskosität des verwendeten Öls verschieden und betragen Rey = 1619 
bei 60°C und Re, = 2639 bei 80°C. Der Verlauf von Nu* in Abhängigkeit der Drehzahl für 
ÙV, = 60°C ist in Abbildung 6.15a dargestellt. 


Über den betrachteten Drehzahlbereich ist ein kontinuierlicher Rückgang des Wärmeübergangs 
mit Erhöhung der Drehzahl bis auf Rew = 84 - 10* zu beobachten. Während bei Reo = 28 - 10* 
mittlere Nußeltzahlen von etwa 0,8 gemessen werden, sinkt dieser Wert auf ein Minimum von 
etwa 0,46 bei Rew = 84-104. Aus den Darstellungen des ortsaufgelösten Wärmeübergangs 
in Abbildungen 6.15b bis 6.15e wird deutlich, dass wie in allen bisher beobachteten Fällen 
das Maximum des Wärmeübergangs mit steigender Drehzahl vom Zahnfuß zum Zahnkopf 
verschoben wird. Der Wärmeübergang an der strömungsabgewandten Zahnflanke nimmt, wie 
auch bei den Messungen bei & = 80°C, kontinuierlich ab und führt zu einer Abnahme des 
mittleren Wärmeübergangs bei erhöhten Drehzahlen. 


Die Verschiebung des Wärmeübergangsmaximums auf der strömungszugewandten Zahnflanke 
zum Zahnkopf beginnt im Vergleich zu den Messungen bei V, = 80 °C und Reg = 2639 (Abb. 6.5) 
bei wesentlich niedrigeren Drehzahlen bzw. Aufprallwinkeln. Bereits bei Reg = 70 10* und 
P = 58,73° (BP123 Abbildung 6.15d) wird im Bereich des Zahnkopfes ein erhöhter Wärmeüber- 
gang gemessen. Diese Verschiebung wurde bei ö = 80°C erst bei B ~ 67° beobachtet. Eine 
mögliche Ursache ist in der geringeren Strahlreynoldszahl zu sehen, bei der noch keine Ablösung 
des Strahls zu erwarten ist. Hierdurch sind die Strahlgeschwindigkeit u, und der tatsächliche 
Aufprallwinkel stets geringer als für abgelöste Strahlen, was zu einer geringeren Eindringtiefe 
führt. 


Der mittlere Wärmeübergang ist für BP123 im Vergleich zu niedrigeren Drehzahlen deutlich 
reduziert. Der Wärmeübergang an der streömungsabgewandten Zahnflanke ist bereits erheblich 
zurückgegangen (siehe Abb. 6.15d). In Betriebspunkten bei 80 °C mit ähnlichen Aufprallwinkeln, 
bei denen beginnender kopflastiger Wärmeübergang gemessen wird, ist der Wärmeübergang auf 
der strömungsabgewandten Flanke noch wesentlich stärker ausgeprägt (vgl. BP29, Abb. 6.7c 
und BP43, Abb. 6.5d). Dies kann mit der durch geringere Temperaturen erhöhten Viskosität des 
Öls zusammenhängen, die das Abprallen des Öls durch Dissipationsvorgänge hemmt. Für eine 
Überprüfung dieses Effektes sind jedoch numerische Strömungsuntersuchungen erforderlich. 
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(e) BP126, Reg = 112-10*, B = 70,64° 


(d) BP123, Rew = 70: 10*, B = 58,73° 


Abbildung 6.15: Einfluss der Drehzahl Rew auf den flächengemittelten Wärmeübergang Nu* 
bei Rey = 1619 und den ortsaufgelösten Warmeiibergang Nu bei reduzierter Öltemperatur 
Ù, = 60°C 
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6.3 Ableitung einer Korrelation aus den gewonnenen Daten 


Aus der Diskussion der Ergebnisse in Abschnitt 6.2 wird deutlich, dass alle der untersuchten 
Parameter 


e Drehzahl n 
e Olvolumenstrom V 
¢ Düsenwinkel & 


e Oltemperatur % 


den Wärmeübergang am Zahnrad beeinflussen. Jedoch zeigen die Verläufe des mittleren Wär- 
meübergangs Nu ein ähnliches Verhalten, wenn sie in Abhängigkeit des Aufprallwinkels ß 
betrachtet werden, der aus V, n und & berechnet wird. Für die Entwicklung von Höchstleistungs- 
getrieben ist es von Interesse, die im Rahmen dieser Arbeit gewonnenen Daten in vereinfachter 
Form abzubilden. Dies legt es nahe, die gemessenen mittleren Wärmeübergangskoeffizien- 
ten durch eine Korrelation auszudrücken, die die beobachteten charakteristischen Einflüsse 
berücksichtigt. Eine solche Korrelation soll im Folgenden beschrieben werden. 


Die Messung des ortsaufgelösten Wärmeübergangs ergibt qualitativ ähnliche Verteilungen bei 
gleichen Werten für ß. Außerdem wird durch die Analyse des Stands der Wissenschaft zur 
Prallkiihlung durch Flüssigkeitsstrahlen deutlich, dass der Wärmeübergang an unbewegten 
ebenen Oberflächen stets eine Proportionalität zu Re? und somit zur Wurzel des Impulses des 
Strahls besitzt. Wird im Fall der Zahnradkühlung die Relativgeschwindigkeit u, herangezogen, 
kann der Ansatz 


h=w”:f(ß) (6.5) 


motiviert werden. Der Wärmeübergang wird nach Gleichung 6.5 als Produkt aus ue? und einer 
Funktion f, die nur vom Aufprallwinkel B abhängt, beschrieben. 


Der Verlauf von h/(u,)° in Abhängigkeit von ß ist in Abbildung 6.16 für alle gemessenen 
Betriebspunkte bei 0; = 80°C dargestellt. Der bezogene Wärmeübergang weist eine starke 
Abhängigkeit von ß auf. Dabei wird zwischen anliegenden und abgelösten Düsenströmungen 
unterschieden. Als anliegende Strömungen werden solche mit Re, < 2251 angesehen, bei denen 
in den Messungen keine abgelöste Düsenströmung beobachtet wird. Bei diesen anliegenden 
Strömungen konnte festgestellt werden, dass der Aufprallwinkel, ab dem der Rückgang des 
Wärmeübergangs einsetzt, kleiner als der bei abgelösten Strömungen ist. Der Aufprallwinkel, ab 
dem dieser Rückgang auftritt, wird im Folgenden als kritischer Aufprallwinkel, Ber, bezeichnet. 


Sowohl für anliegende als auch für abgelöste Düsenströmungen wird ein ansatzweise linearer 
Verlauf von h/(u,) bis zu Ber beobachtet. Bei Überschreiten dieses Wertes ist eine stärkere 
Abnahme von h/(u,)° festzustellen, die durch ein Polynom zweiten Grades beschrieben werden 
kann. Dieser charakteristische Verlauf ist auf die in Kapitel 6.2.1 beschriebenen Vorgänge 
zurückzuführen. Die kritischen Aufprallwinkel für anliegende und abgelöste Düsenströmungen, 
Ber und Born, unterscheiden sich. 
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Abbildung 6.16: Werte des Ausdrucks h/(u,)° in Abhängigkeit von ß für alle gemessenen 
Betriebspunkte bei 9% = 80°C 


Unter Berücksichtigung der Form der Düsenströmung und Unterteilung in unter- und über- 
kritische Aufprallwinkel können alle Datenpunkte durch den Ansatz aus Gleichung 6.5 gut 
wiedergegeben werden. Dazu wird für beide Formen der Düsenströmung jeweils ein unter- und 
überkritischer Bereich definiert, deren Grenzen durch die jeweiligen kritischen Aufprallwinkel 
definiert sind, mit 

Ber = 46° 


für anliegende Düsenströmungen und 


Bern = 57° 


für abgeléste Düsenströmungen. Der unterkritische Bereich B < Per wird dabei durch ein 
Polynom ersten Grades, der überkritische Bereich durch ein Polynom zweiten Grades abgebildet: 


h = u” -(pıß?+p2ß + p3). (6.6) 


Durch eine Fehlerquadratminimierung können die Polynomkoeffizienten zu Gleichung 6.6 
gefunden werden, die in Tabelle 6.2 aufgeführt sind. 


Aus dem Vergleich der berechneten Wärmeübergangskoeffizienten Ap mit den gemessenen 
Werten h in Abbildung 6.17 geht hervor, dass die Beschreibung des mittleren Wärmeübergangs 
durch Gleichung 6.6 und die Koeffizienten aus Tabelle 6.2 im Allgemeinen eine sehr gute 
Übereinstimmung mit den gemessenen Werten aufweist. Die Standardabweichung des Fehlers 
beträgt dabei o = 10,7%. 


Der durch Gleichung 6.6 beschriebene Zusammenhang unterliegt der Einschränkung, dass er 
nur für die untersuchte Düsenkonfiguration gültig ist. Eine Übertragung des Wärmeübergangs 


100 Ergebnisse 


Tabelle 6.2: Polynomkoeffizienten des stiickweisen Polynomansatzes nach Gleichung 6.6 


Düsenströmung B< Per B > Pa 
pı=0 pı = 0,06391 

anliegend p2 = —0,8798 p= —9,829 
p3 = 123,5 p3 = 403,2 
pı=0 pı = 0,09888 

abgelöst p2 = —0,7135 p2 =-—15,78 
p3 = 125 p3 = 663,8 


im Bereich der Stauströmung für geänderte Diisendurchmesser ist denkbar. Die Auswirkungen 
von Veränderungen des Diisendurchmessers oder der Abstände zwischen den Öldüsen werden 
jedoch nicht durch die im Stand der Wissenschaft wiedergegebenen Korrelationen erfasst. Diese 
Einschränkung kann aufgehoben werden, sobald experimentelle oder auch numerische Daten 
mit entsprechend geänderten Düsenkonfigurationen verfügbar sind. 


Amax 


Abbildung 6.17: Vergleich der gemessenen und der durch die Korrelation berechneten mittleren 
Wärmeübergangskoeffizienten, jeweils bezogen auf den maximalen gemessenen mittleren Wär- 
meübergangskoeffizienten Ras 


Die getroffene Aufteilung der Düsenströmung in anliegende und abgelöste Bereiche ist not- 
wendig, da die Berechnung von % und ß bei abgelösten Strömungen durch die geänderte 
Strahlgeschwindigkeit fehlerhaft ist. Die tatsächliche Strahlgeschwindigkeit wird im Rahmen 
dieser Arbeit nicht gemessen. Unter Verwendung gemessener Strahlgeschwindigkeiten und 
dadurch korrigierter Berechnungen besteht die Möglichkeit, dass die Verläufe beider Bereiche 
zusammenfallen und somit nur eine Aufteilung in Bereiche unter- und überkritischer Aufprall- 
winkel getroffen werden muss. Eine Berechnung der Strömungsgeschwindigkeit ausschließlich 
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anhand der vorliegenden Daten zur Düse und der Strömung ist nicht möglich. 
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7 Zusammenfassung 


Die Verbesserung von Triebwerken hinsichtlich ihres Treibstoffverbrauchs ist durch jahrzehnte- 
lange Forschung weit vorangeschritten. Klassische Optimierungspotentiale sind dadurch weitest- 
gehend ausgeschöpft. Weitere, inkrementelle Wirkungsgradsteigerungen in diesen Bereichen, 
etwa durch verbesserte Materialien oder Kühlmethoden, sind daher enorm kostenintensiv. Das 
Konzept des Getriebefantriebwerks hat das Potential, einen disruptiven Entwicklungsfortschritt 
auszulösen. Jedoch stellt das benötigte Höchstleistungsgetriebe eine neue und hochkomplexe 
Komponente dar, deren Entwicklung eine technologische Herausforderung ist. Das zur Umset- 
zung des Konzepts notwendige Getriebe muss hohen und sehr gegensätzlichen Anforderungen 
genügen. Die Leistungsdichte muss durch die Forderung nach geringem Gewicht und kleinem 
Bauraum extrem hoch sein. Als Konsequenz ist die Wärmestromdichte, die durch Reibungs- 
vorgänge im Zahnradeingriff entsteht, hoch. Eine wesentliche Herausforderung besteht in der 
Kühlung der hochbelasteten Zahnräder. 


Das Ziel dieser Arbeit war, ein tiefgreifendes Verständnis der Zahnradkühlung durch Öleinsprit- 
zung zu gewinnen und so einen Beitrag zur Entwicklung von Höchstleistungsgetrieben zu leisten. 
Dazu wurde der Einfluss der Betriebsparameter Drehzahl, Ölvolumenstrom, Düsenwinkel und 
Öltemperatur auf den ortsaufgelösten Wärmeübergang am Zahnrad analysiert. Einerseits wurden 
dadurch die relevanten Zusammenhänge aufgedeckt, die den Wärmeübergang beeinflussen und 
durch deren Kenntnis eine effiziente Auslegung des Kühlsystems erst möglich wird. Andererseits 
wurde so eine umfangreiche Datenbasis mit exakt definierten Randbedingungen geschaffen, 
durch die in Zukunft numerische Modelle validiert werden können. Erst durch diese Datenbasis 
sind in Zukunft kostensparende CFD-Simulationen des Öl-Zahnrad-Wärmeübergangs möglich. 


Bisherige Arbeiten und Prüfstände liefern keine ortsaufgelösten Messdaten des Wärmeübergangs. 
Zusätzlich ist durch das gemeinsame Auftreten von Reibung und Kühlung die Bestimmung von 
Wärmeübergangskoeffizienten in diesen Arbeiten nur eingeschränkt möglich. Daher wurde in 
dieser Arbeit eine neue Methode entwickelt und vorgestellt, mit der der Wärmeübergang an 
einem Zahnrad erstmals getrennt von Reibwärmeeinflüssen und ortsaufgelöst gemessen werden 
konnte. Durch die Analyse der Methode durch Monte-Carlo-Simulationen konnte deren hohe 
Messgenauigkeit nachgewiesen werden. Durch die Daten der Monte-Carlo-Simulationen konnte 
zusätzlich die Robustheit der ortsaufgelösten Messung bewiesen werden. 


Mit einem neu entwickelten und am Institut für Thermische Strömungsmaschinen gefertigten 
Prüfstand wurde der Wärmeübergang am Zahnrad für eine Konfiguration der Einspritzkühlung 
und eine Vielzahl von Betriebspunkten untersucht. Dabei wurde ein komplexer Zusammenhang 
zwischen den Betriebsparametern und dem Wärmeübergang beobachtet, der durch bisherige 
Arbeiten nicht abgebildet wird. Es zeigte sich die Existenz einer kritischen Drehzahl, bei deren 
Überschreitung ein deutlicher Rückgang des Wärmeübergangs zu verzeichnen ist. Durch Verwen- 
dung eines kinematischen Modells des Ölaufpralls konnten die beobachteten Zusammenhänge 
erklärt werden. Der Aufprallwinkel, unter dem der Ölstrahl auf die Zahnflanke auftrifft, hat 
dabei eine herausragende Bedeutung für den zu erwartenden mittleren Wärmeübergang. Durch 
die identifizierten Zusammenhänge wurde eine Korrelation abgeleitet, die den mittleren Wär- 
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meübergang im bisher erfassten Parameterbereich der Einspritzkühlung durch einen einfachen 
mathematischen Ansatz in sehr guter Näherung abbildet. 


Durch die Thermalmodellierung des Zahnrads unter Berücksichtigung der Reibung kann gezeigt 
werden, dass nicht nur der mittlere Wärmeübergang, sondern auch seine örtliche Verteilung 
zur Einhaltung von Temperaturgrenzen von Bedeutung ist (Bauer (2018)). Daten örtlich auf- 
gelöster Messungen des Wärmeübergangs sind durch die vorliegende Arbeit erstmalig und für 
einen weiten Parameterbereich verfügbar. Dies versetzt den Entwickler in die Lage, eine aus- 
reichende Kühlung der Zahnflanken sicherzustellen, ohne dabei Wirkungsgradverluste durch 
zu hohe eingesetzte Ölvolumenströme zu riskieren. Dadurch, dass somit die Versorgungs- und 
Übertragungssysteme im Ölkreislauf des Triebwerks kleiner dimensioniert werden können, 
stellt dies einen wertvollen Baustein für die Entwicklung eines Höchstleistungsgetriebes für das 
Getriebefantriebwerk dar. 


Um eine Basis für eine allgemeingültige Korrelation für den Wärmeübergang zu schaffen, muss 
der Wärmeübergang bei geänderten Geometrien untersucht werden. Durch die in dieser Arbeit 
entstandene Methode sind die dafür notwendigen Grundlagen geschaffen worden. Schrägverzah- 
nung, die Zähnezahl und der Modul sind Einflussgrößen, deren Wirkung bisher nicht bekannt ist 
und in zukünftigen Untersuchungen bestimmt werden muss. 


Das grundlagenorientierte Wissen zur instationären Prallkühlung bewegter Oberflächen ist sehr 
begrenzt. Die Beschreibung des Öl-Zahnrad-Wärmeübergangs als Modell aus Teilprozessen ist 
bisher nicht möglich, da diese Teilprozesse in der Literatur nur unvollständig beschrieben sind. 
Wenn durch zukünftige Arbeiten das Grundlagenwissen erweitert wird, kann die Modellbeschrei- 
bung des Öl-Zahnrad-Wärmeübergangs möglich werden. Durch die umfangreiche Datenbasis, 
die mit dieser Arbeit vorliegt, können solche Modelle validiert werden. Die Übertragung der 
Ergebnisse dieser Arbeit auf Zahnräder mit anderer Geometrie, aber auch weitere technische 
Fragestellungen im Kontext der Prallkühlung, kann dadurch ermöglicht werden. 
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A.1 Gemessene Betriebspunkte und Wärmeübergangsverteilungen 


A.1.1 Tabellarische Darstellung der Parameter 


Tabelle A.1: Parameter und Nu* der gemessenen Betriebspunkte 


BP n a Vs V Nu* 
1/min | ° °C l/min | - 
BPI 1988 | 10,78 | 76,46 | 1,56 | 0,140 
BP2 1988 10,76 | 78,04 | 2,45 0,219 
BP3 1988 | 10,77 | 79,29 | 3,28 | 0,300 
BP4 1988 | 23,49 | 79,57 | 3,28 | 0,298 
BP5 1988 | 30,06 | 79,71 | 3,28 | 0,295 
BP6 1988 | -0,07 | 79,85 | 3,28 | 0,321 
BP7 2965 | 11,27 | 80,53 | 3,26 | 0,269 
BP8 2804 | 12,72 | 68,82 | 3,9 0,286 
BP9 2804 11,74 | 71,02 | 4,38 0,330 
BP10 2804 11,79 | 71,92 | 4,59 0,441 
BP11 2804 11,77 | 72,6 5,22 0,620 
BP12 | 2804 | 11,8 73,34 | 6,86 | 0,823 
BP13 3099 11,77 | 75,12 | 4,58 0,423 
BP14 1995 11,79 | 76,27 | 3,32 0,318 
BP15 1998 20,33 | 77,51 | 2,21 0,205 
BP16 | 2805 | 20,21 | 782 | 683 | 0,852 
BP17 3101 20,21 | 79,08 | 8,36 0,943 
BP18 | 3785 | 20,29 | 80,05 | 7,14 | 0,904 
BP19 | 5030 | 19,53 | 81,43 | 7,09 | 0,677 
BP20 | 896 19,93 | 83,02 | 5,17 | 0,732 
BP21 897 20,12 | 82,47 | 3,71 0,646 
BP22 | 896 20,17 | 80,98 | 2,54 | 0,536 
BP23 896 19,97 | 78,46 | 1,14 0,170 
BP24 | 1568 | 30,98 | 78,4 | 2,85 | 0,301 
BP25 | 3783 | 19,57 | 79,55 | 10,33 | 1,096 
BP26 | 5054 | 19,29 | 80,74 | 10,28 | 1,265 
BP27 1863 19,22 | 82,14 | 4,26 0,568 
BP28 | 1998 | 19,51 | 80,05 | 3,19 | 0,302 
BP29 1499 19,48 | 80,31 | 3,21 0,402 
BP30 | 1002 | 19,49 | 80,47 | 3,23 | 0,564 
BP31 | 2981 | 19,39 | 80,67 | 3,19 | 0,289 
BP32 | 1996 | 19,4 81,9 | 4,84 | 0,738 
BP33 | 991 19,32 | 82,13 | 4,84 | 0,674 
BP34 1508 19,38 | 82,29 | 4,84 0,724 
BP35 2994 19,27 | 82,23 | 4,84 0,582 
BP36 | 2993 | 19,17 | 83,87 | 723 | 1,001 
BP37 1996 19,02 | 84,21 | 7,21 0,955 
BP38 | 1500 | 19,2 83,67 | 7,23 | 0,931 
BP39 | 1003 | 19,21 | 84,04 | 7,23 | 0,895 
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BP n a Vs V Nu* 
1/min | ° °C l/min | - 

BP40 1997 19,24 | 81,42 | 3,16 0,354 
BP41 | 2993 | 30,5 81,51 | 4,69 | 0,512 
BP42 2993 44,49 | 81,87 | 4,68 0,476 
BP43 | 2993 | 19,94 | 81,98 | 4,68 | 0,545 
BP44 | 2993 | 0,27 82,15 | 4,68 | 0,585 
BP45 | 3000 | -10,48 | 82,27 | 4,68 | 0,616 
BP46 | 2992 | -20 82,34 | 4,68 | 0,609 
BP47 | 2992 | -29,91 | 82,35 | 4,68 | 0,576 
BP48 | 2993 | 10,63 | 82,59 | 4,67 | 0,605 
BP49 | 3985 | 20 82,62 | 4,68 | 0,484 
BP50 | 3986 | 19,68 | 82,02 | 3,15 | 0,338 
BP51 2992 19,43 | 85,51 | 9,51 1,225 
BP52 1500 20,07 | 79,96 | 3,24 0,414 
BP53 | 1002 | 20,08 | 80,48 | 3,24 | 0,583 
BP54 | 1998 | 20,05 | 80,68 | 3,23 | 0,328 
BP55 | 2499 | 20,07 | 80,83 | 3,24 | 0,303 
BP56 | 3006 | 20,05 | 80,19 | 3,24 | 0,301 
BP57 | 3501 | 20,07 | 80,01 | 3,23 | 0,301 
BP58 | 4002 | 20,05 | 80,28 | 3,23 | 0,309 
BP59 | 748 20,08 | 80,62 | 3,23 | 0,595 
BP60 | 1998 | 20,05 | 81,01 | 4,11 | 0,482 
BP61 2496 20,02 | 81,1 4,11 0,396 
BP62 2977 20,05 | 81,11 | 4,11 0,373 
BP63 1499 20,06 | 81,08 | 4,11 0,717 
BP64 989 20,11 | 81,14 | 4,12 0,724 
BP65 | 3503 | 20,03 | 81,07 | 4,11 | 0,364 
BP66 | 4001 | 20,03 | 80,93 | 4,11 | 0,361 
BP67 1997 20,07 | 81,45 | 4,77 0,794 
BP68 | 2496 | 20,05 | 81,43 | 4,77 | 0,769 
BP69 | 2750 | 20,04 | 81,5 | 4,77 | 0,680 
BP70 | 3493 | 20,06 | 81,53 | 4,77 | 0,513 
BP71 | 1488 | 20,07 | 80,94 | 5,96 | 0,829 
BP72 1004 20,02 | 81,14 | 5,95 0,831 
BP73 | 1998 | 20,04 | 81,33 | 5,95 | 0,872 
BP74 | 2491 | 20,06 | 81,35 | 5,95 | 0,892 
BP75 | 2993 | 20,03 | 81,32 | 5,95 | 0,930 
BP76 | 3491 | 20,09 | 81,48 | 5,94 | 0,790 
BP77 | 3988 | 20,03 | 81,59 | 5,94 | 0,656 
BP78 | 1998 | 19,59 | 79,34 | 3,24 | 0,304 
BP79 | 1998 | -0,28 | 79,69 | 3,25 | 0,325 
BP80 | 1493 | -0,31 | 79,87 | 326 | 0,451 
BP81 | 1002 | -0,29 | 80,02 | 3,25 | 0,586 
BP82 | 747 -0,34 | 80,1 | 3,25 | 0,592 
BP83 | 2501 | -0,35 | 80,21 | 3,25 | 0,317 
BP84 | 3002 | -0,3 80,22 | 3,23 | 0,307 
BP85 | 3500 | -0,3 80,24 | 3,23 | 0,305 
BP86 | 4001 | -0,27 | 80,25 | 3,23 | 0,309 
BP87 1499 -0,19 80,13 | 4,74 0,787 
BP88 | 1997 | -0,18 | 80,29 | 4,74 | 0,808 
BP89 | 2495 | -0,22 | 80,37 | 4,73 | 0,720 
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BP n a Vs V Nu* 
l/min | ° °C l/min | - 

BP90 | 2752 | -0,21 | 80,52 | 4,73 | 0,659 
BP91 | 3003 | -0,23 | 80,67 | 4,73 | 0,594 
BP92 | 3499 | -0,22 | 80,66 | 4,72 | 0,514 
BP93 | 4000 | -0,18 | 80,63 | 4,73 | 0,486 
BP94 1001 -0,2 80,84 | 4,73 0,748 
BP95 | 1996 | -0,23 | 81,52 | 5,99 | 0,989 
BP96 | 2490 | -0,19 | 81,59 | 5,99 | 1,002 
BP97 | 3000 | -0,2 81,54 | 5,99 | 0,971 
BP98 | 3501 | -0,22 | 81,65 | 5,99 | 0,844 
BP99 | 3994 | -0,2 81,72 | 5,99 | 0,720 
BP100 | 4321 | -0,22 | 81,82 | 5,99 | 0,666 
BP101 | 1001 | -0,22 | 81,96 | 5,98 | 0,840 
BP102 | 1997 44,96 | 81,23 | 4,75 0,665 
BP103 | 2496 | 44,98 | 81,05 | 4,75 | 0,623 
BP104 | 3002 44,91 | 80,94 | 4,75 0,481 
BP105 | 3503 | 44,94 | 80,84 | 4,74 | 0,430 
BP106 | 3991 | 44,93 | 80,7 | 4,75 | 0,406 
BP107 | 1500 44,87 | 80,91 | 4,75 0,758 
BP108 | 1002 | 44,91 | 80,85 | 4,75 | 0,820 
BP109 | 1999 | 44,83 | 80,46 | 7,18 | 0,953 
BP110 | 2497 44,82 | 80,64 | 7,17 0,868 
BP111 | 2995 44,78 | 80,86 | 7,18 0,792 
BP112 | 3505 44,87 | 80,87 | 7,18 0,819 
BP113 | 4002 44,5 81,07 | 7,17 0,750 
BP114 | 1501 44,6 81,1 7,17 1,030 
BP115 | 1002 44,56 | 81,1 7,18 0,973 
BP116 | 1999 -0,2 81,13 | 7,17 1,081 
BP117 | 2495 -0,2 81,14 | 7,17 1,185 
BP118 | 2992 -0,19 80,99 | 7,18 1,179 
BP119 | 1001 | 19,56 | 61,34 | 4,66 | 0,798 
BP120 | 1500 | 19,55 | 60,19 | 4,67 | 0,698 
BP121 | 1998 | 19,58 | 60,59 | 4,67 | 0,628 
BP122 | 1986 19,56 | 61,4 | 4,65 0,717 
BP123 | 2496 | 19,61 | 60,89 | 4,66 | 0,501 
BP124 | 2993 19,57 | 61,12 | 4,66 0,461 
BP125 | 3491 19,62 | 61,16 | 4,65 0,467 
BP126 | 3986 | 19,54 | 61,13 | 4,66 | 0,495 
BP127 | 3994 | -0,7 78,38 | 6,19 | 0,781 
BP128 | 3684 -0,21 79,63 | 7,54 1,265 
BP129 | 3998 | 0,24 79,86 | 7,07 | 0,980 
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Abbildung A.1: Ortsaufgelöste Wärmeübergangsmessungen für BP1 bis BP6 
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(e) BP11 (f) BP12 


Abbildung A.2: Ortsaufgelöste Wärmeübergangsmessungen für BP7 bis BP12 
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Abbildung A.3: Ortsaufgelöste Wärmeübergangsmessungen für BP13 bis BP18 
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Abbildung A.4: Ortsaufgelöste Wärmeübergangsmessungen für BP19 bis BP24 
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Abbildung A.5: Ortsaufgelöste Wärmeübergangsmessungen fiir BP25 bis BP30 
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Abbildung A.6: Ortsaufgelöste Wärmeübergangsmessungen fiir BP31 bis BP36 
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Abbildung A.7: Ortsaufgelöste Wärmeübergangsmessungen für BP37 bis BP42 
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Abbildung A.8: Ortsaufgelöste Wärmeübergangsmessungen für BP43 bis BP48 
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Abbildung A.9: Ortsaufgelöste Wärmeübergangsmessungen für BP49 bis BP54 
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Abbildung A.10: Ortsaufgelöste Wärmeübergangsmessungen für BP55 bis BP60 
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Abbildung A.11: Ortsaufgelöste Wärmeübergangsmessungen für BP61 bis BP66 
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Abbildung A.12: Ortsaufgelöste Wärmeübergangsmessungen für BP67 bis BP72 
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Abbildung A.13: Ortsaufgelöste Wärmeübergangsmessungen für BP73 bis BP78 
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Abbildung A.14: Ortsaufgelöste Wärmeübergangsmessungen für BP79 bis BP84 
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Abbildung A.15: Ortsaufgelöste Wärmeübergangsmessungen fiir BP85 bis BP90 
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Abbildung A.16: Ortsaufgelöste Wärmeübergangsmessungen für BP91 bis BP96 
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Abbildung A.17: Ortsaufgelöste Wärmeübergangsmessungen für BP97 bis BP102 
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Abbildung A.18: Ortsaufgelöste Wärmeübergangsmessungen für BP103 bis BP108 
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Abbildung A.19: Ortsaufgelöste Wärmeübergangsmessungen für BP109 bis BP114 
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Abbildung A.20: Ortsaufgelöste Wärmeübergangsmessungen für BP115 bis BP120 
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Abbildung A.21: Ortsaufgelöste Wärmeübergangsmessungen für BP121 bis BP126 
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Abbildung A.22: Ortsaufgelöste Wärmeübergangsmessungen für BP127 bis BP129 
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A.2 Stoffwertkorrelationen des verwendeten Öls 


Im Rahmen dieser Arbeit werden die von Glahn (1995) angegebenen Stoffwertkorelationen für 
Mobil Jet Oil II verwendet. In den angegebenen Zahlenwertgleichungen ist Us; in K einzusetzen. 


Kinematische Viskosität v in 10~° m?/s: 


loglog(v + 0,8) = 8,962 182 — 3,527033 log %,. (Al) 


Dichte p in kg/m?: 
p = 1329,44 — 1,164 71%. (A.2) 


Warmeleitfahigkeit A in W/(mK): 


A = 0,15799 — 0,375: 1074 %,. (A.3) 


Spezifische Wärme c, in J/(kgK): 


Cp = 1054,86 — 2,88029 ¥,. (A.4) 


A.3 Herleitung des spezifischen Relativimpulses 


Die Relativgeschwindigkeit u,, mit der das Öl auf die Zahnflanke auftrifft, kann durch ihren 
radialen Anteil u, und tangentialen Anteil ur į berechnet werden: 


Uy = 4/ U2, +u2.. (A.5) 


Die Anteile sind jeweils die Summe der radialen und tangentialen Geschwindigkeitsanteile 
des Ölstrahls und des Zahns. Die Eigengeschwindigkeit des Zahns u; besitzt keine radiale 
Komponente. Somit wird Gleichung A.5 zu: 


ur = Vey + (uz — ust). (A.6) 


Die radialen und tangentialen Anteile von us können durch den Düsenwinkel & ausgedrückt 
werden: 


Üp = ye “cos? + (uz — us- sina)”. (A.7) 


Durch Umformen ergibt sich: 


Ur = y? — 2uzussin& + u?sin?@& + u2cos” a 
aaa 


=us 


= y — 2uzus + u2 + 2uzus (1 — sina) 


= V (uz — us)? + 2uzus (1 — sing). (A.8) 
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Durch Ersetzen durch us = Zu, und Umformen wird Gleichung A.8 zu: 


ur = yea — £)? + 2¥u2 (1 — sina) 


= wl (I =>) 99 901 — sina) 


= uz / 1 +32 —2 (1 — sine) 


und somit: 


“ = \/1+22-22(1-sine). (A.9) 
uz 
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